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STUDY OF CRYSTAL STRUCTURES AND PHASE TRANSITIONS OF 
SOLfD DEUTERIUM IODIDE 



Abstract 

r 

The crystal structure of solid deuterium Iodide was studied by 
thermal neutron and X-ray diffraction at tiiree temperatures. At ^.2K the 
low phase structure is ordered, monocMnlc space group C2/c, with 
16 molecules in the unit cell. At 8^K a two-fold disordered orthorhombic 
structure, space group Bbcm, is proposed with four molecules in the unit 
cell, each deuterium atom oscillating between two equivalent positions 180° 
apart. This middle phase is isomorphous with deuterium bromide. 

At 137K a twelve-fold disordered cubic structure Is proposed, space 
group Fm3m, with four molecules In the face-centred, unit cell. The 
deuterium atom occupies In turn twelve statistically equivalent positions. 
This structure is isomorphous with the high phases of deuterium chloride and 
bromide. 

Phase transitions between the monocUnic and orthorhombic phases at 
77K and the orthorhombic and cubic phases at 127K were studied by neutron 
diffraction in the temperature Intervals ^5-80K and I10-130K, Order 
parameters for these transitions have been calculated from the neutron 
diffraction data and agreement obtained with those calculated from calori- 
metrle data for the lower transition. 



ACKNOWLEDGEMENTS 

It is a pleasure to thank Dr E Sa'ndor, my supervisor, for his guidance 
and encouragement during the course of this worl<.' 

I would also like to thank: C Hotston, P Cruze and R Putney of 
Queen Mary College for technical assistance, Professor G Jones-, 
Professor A Ashmore and Professor D H Martin for the research facilities at 
Queen Mary College and Mr D H C Harris and his staff for assistance with the 
neutron diffraction experiments at A.E.R.E., Harwell. 

Finally I thank the Science Research Council for a research student- 
ship and the Neutron Beam Research Committee for providing neutron 
faci 1 i ties. 



TABLE OF CONTENTS 

Page 
Abstract 2 

Acknowledgements 3 

Table of contents k 

List of figures 7 

List of tables 15 

CHAPTER O N E - Review of previous work on HI and PI 18 

1.1 Introduction 18 

1.2 Calorimetric experiments 19 

1.3 X-ray experiments 21 
].h Optical experiments 23 
1-5 Dielectric experiments 2k 

(a) Static dielectric constant measurements 2A 

(b) Complex dielectric constant measurements 27 

1.6 Nuclear magnetic resonance experiments (NMR) 35 

1.7 Infra-red spectroscopy kS 

1.8 Theories of phase transitions In molecular crystals ^9 

CHAPTER TV/0 - E xperimental work 52 

2.1 Chemical preparations 52 

(a) Indirect preparation of HI/DI 53 

(b) Direct preparation of Dl 55 

2.2 X-ray diffraction study of Dl 57 

2.3 Neutron diffraction study of Dl 65 

(a) Dido 'Curran' powder di ff ractometer 66 

(b) Dido 'Badger' pavder d i ff ractometer 69 

(c) Pluto 'Panda' powder d i f f rfiClometer 69 
2.h CryostaLs used in neutron diffraction experiments 7^ 

(a) T.B.T. cryostat 7^1 

(b) B.O.C. cryostat 7^ 



Page 

CHAPTER THREE - X-ray povjde r diffraction study of the iodine 79 

sublattice in the three solid phases of PI 

3-1 Introductory theory of diffraction 79 

3.2 X-ray diffraction experiments 83 

(a) The low phase of Dl 83 

(b) The middle and high phases of DI 87 

CHAPTER FOUR - Neutron powder diffraction study of the crystal 92 

structure of the low phase of Dl at ^.2K 

k.\ Summary of the neutron diffraction experiments 92 

k.2 The unit cell of the low phase 9'* 

4.3 Models based on orthorhombic space groups 99 
k.k Models based on monoclinic space groups 102 
^.5 Refinement of the nxjdel based on space group C2/c (No. 15) 119 
4.6 Analysis of the intensity transform of the neutron powder 13^1 

pattern of the low phase 

^1.7 Conclusion ]^'^ 

CHAPTER FIVE - Neutron powder diffraction study of the crystal ] 47 
structure of the middle phase of Dl 

CHAPTER SIX - N'eutron powder diffraction study of the crystal 155 
structure of the high phase of Dl 

6.1 Zinc-sulphlde-type nwdel space group Fm3m (No. 225) 157 

6.2 Sodium-chlorlde-type model space group Fm3m (No. 225) 157 

6.3 Four-fold disordered model space group FT3m (No. 216) 158 
6A Six-fold disordered model space group Fm3m (No. 225) 158 

6.5 Eight-fold disordered model space group Fm3 (No. 202) 1 5B 

6.6 Twelve-fold disordered model space group Fm3m (No. 225) 158 

6.7 Freely-rotating or randomly disordered model 153 

6.8 Comparison of observed structure factors with structure 1 60 
factors calculated from disordered models 



Page 

CHAPTER SEVEN - Neutron powder diffraction study of the tvjo 1 70 

phase transitions in solid PI 

7. J The lower phase transition 170 

7.2 The upper phase transition 1 80 

7.3 Conclusion 192 

CHAPTER EIGHT - Summary 194 

8.1 Conclusions of the earlier experimental studies concerning 19^ 
the crystal structure of Dl 

8.2 The crystal structure of Dl as determined by X-ray and 196 
neutron diffraction 

8.3 The nature of the two phase transitions in solid Dl 197 

8.4 Suggested further experiments on HI/DI 199 

(a) Neutron powder diffraction 199 

(b) Neutron single-crystal diffraction , 200 

(c) Inelastic neutron scattering 200 

APPENDIX ONE - Computer programs 202 

Al . 1 Least-squares refinement of crystal structures 202 

Al.2 Least-squares refinement of cell dimensions 203 

AI.3 Calculation of powder peak positions 204 

A1.4 Calculation of random-orientation structure factors 204 

Al . 5 LPSG, Low phase symmetry generator 204 

A1.6 The RIetveld program 205 

APPENDIX TWO - Temperature measurement and control circuit 207 
of the variable temperature cryostat 

List of References 212 



List of figures 
Chapter o ne 

1.1 Specific heats of HI/DI at constant pressure as a function of 
temperature, Clusius and V/olf (19^7)- 

1.2 Static dielectric constants of HCl, HBr/DBr and H!/DI as a function 
of temperature, Cole and Havrilial< (1357), Brown and Cole (1953) 
and Swenson and Cole (ISS'*)- 

1.3 Imaginary versus real dielectric constant of HI/DI below the lower 
phase transition point. Cole and Havriliak (1957). 

].k Imaginary versus real dielectric constantof HI/DI above the lower 
phase transition point. Cole and Havriliak (1957)- 

1.5 Absorption conductance of D! and DBr as a function of frequency, 
Cole and Havri 1 iak (1957). 

1.6 Relaxation times for hydrogen hal ides as a function of temperature 
determined by dielectric absorption methods, Cole and Havriliak (1957) 

1.7 Imaginary versus real dielectric constant of HI as a function of 
temperature, Groenewegen and Cole (1966). 

1.8 Reduced locus of the imaginary dielectric constant versus the real 
part, Groenewegen and Cole (I966). 

1.9 Relaxation frequency as a function of temperature for HBr and HI, 
Groenewegen and Cole (I966). 

1.10 NHR line width versus temperature for HI, Alpert {]SkS) . 

1.1) Spin-lattice relaxation times versus temperature for DC! , DBr and Dl, 
Gcnln et al (i960). 



8. 

1.12 Spin- latt ice, rotating frame and spin-spin relaxation times, and 
second moment data for HCl , HBr and HI as a function of temperature, 
Genin et al (1968) . 

1.13 Relaxation times as a function of temperature derived from NMR and 
dielectric loss measurements on HI, Genin et al (1968) and 
Groenewegen and Cole (I966). 

l.l'* NMR dispersion signal in polycrystal 1 ine Dl at 96K and 7 MHz, Genin 
et al (1968). 

1.15 NQR coupl ing constants of polycrystal 11 ne DC 1 , DBr and Dl , as a 
function of temperature, Genin et al (I968). 

1.16 Infra-red transmission spectra of HI as a function of temperature on 
cooling through the lower phase transition at 70K, Hornig and 
Osberg (1955)- 

1.17 Infra-red transmission spectra of HI as a function of temperature, 
Hiebert and Hornig (1957). 

Chapter two 

2.1 HI and Dl sample preparation system used for the indirect method §2. la. 

2.2 Dl sample preparation system used for the direct method §2. lb. 

2.3 Dl sample purification system used for the direct method §2. lb. 
2. A Low temperature X-ray goniometer used at Queen Mary College. 

2.5 Liquid helium cryostat used on the Queen Mary College X-ray powder 
di ff ractomcter. 

2.6 The configuration of the Dido 'Curran' neutron powder d i f f ractometer. 

2.7 Photograph of the Dido 'Curran' neutron powder di f-^-actometer. 



2.8 The configuration of the Pluto 'Panda' neutron powder di f f ractometer. 

2.9 Photograph of the Pluto 'Panda' neutron powder d i ff ractometer. 

2.10 The variable temperature cryostat used for studying the lower phase 
transition of DI by neutron diffraction. 

Chapter three 

3.1 Low-angle X-ray powder photograph of DI at 84K and 50K. 

3.2 X-ray powder pattern of DI at 4.2K (low phase) chart recorder output. 

3.3 X-ray powder photograph of DI at 117K (middle phase). 
3.'» X-ray powder photograph of DI at I3OK (high phase). 

3.5 Lattice parameters of HI and DI as a function of temperature. 

Chapter four 

k.] Neutron powder diffraction pattern of DI at 4.2K (low phase) recorded 
wi th Di do 'Curran ' . 

4.2 Neutron powder diffraction pattern of DI at 4.2K (low phase) recorded 
with Pluto 'Panda' . 

4.3 Proposed structure for DI low phase space-group Fdd2 (No. 43), model No.l. 

4.4 Proposed structure for DI low phase space-group Fdd2 (No. 43), model No. 2. 

4.5 Proposed structure for DI low phase space-group Fdd2 (No. 43), model No. 3. 

4.6 The six methods of transforming the large orthorhombic unit cell to 
the smaller monocllnic unit cell. 

4.7 Middle phase structure of DI drawn In a monocllnic unit cell, type No,]. 



I»,8 Middle phase structure of Dl drawn in a monoclinic unit cell, type No.? 

A. 9 Middle phase structure of Dl drawn in a monoclinic unit cell, type No. 3. 

it. 10 Middle phase structure of Dl drawn in a monoclinic unit cell, type No. A. 

i».ll Middle phase structure of Dl drawn in a monoclinic unit cell, type No. 5- 

'».12 Middle phase structure of Dl drawn in a nranoclinic unit cell, type No. 6. 

A. 13 Middle phase structure of Dl drawn in a monoclinic unit cell, type No. 3 
projected onto the (OOl) monoclinic plane. 

k.\k Middle phase structure of Dl drawn in a monoclinic unit cell, type No. 4 
projected onto the (OOl) monoclinic plane. 

^.15 Proposed structure for Dl low phase space group Cc (No. 9) model No.l. 

k.](> Proposed structure for Dl low phase space group Cc (No. 9) model No. 2. 

i*.!? Proposed structure for Dl low phase space group Cc (No. 9) model No. 3- 

k.\B Proposed structure for Dl low phase space group Cc (No. 9) model No.^. 

k.\S Proposed structure for Dl low phase space group Cc (No. 9) model No. 5. 

4.20 Proposed structure for Dl low phase space group Cc (No. 9) model No. 6. 

4.21 Calculated (continuous line) and observed (crosses) intensities for Dl 
low phase, 'Curran' 4.2K model No. 3. 

4.22 Calculated (continuous line) and observed (crosses) intensities for Dl 
low phase, 'Curran' 4.2K model No. 4. 

4.23 Calculated (continuous line) and observed (crosses) intensities for Dl 
low phase, 'Curran' 4.2K model No. 5. 

4.24 Calculated (continuous line) and observed (crosses) Intensities for Dl 
low phase, 'Curran' 4.2K model No. 6. 



n 



4.25 Hidcjle phase structure of Dl drawn in an orthorhombic unit cell showing 
the 8 molecules of the basis chosen to produce 128 models. 

k.26 Middle phase structure of Dl drawn in a monoclinic unit cell type No. 6 
showing the 6 molecules of the basis chosen to produce 32 models with 
space group Cm (No. 8). 

^.27 Proposed structure for Dl low phase space group C2/c, model No. 5 in 
unit cell No.^. 

'».28 (a) (b) Observed (crosses) and calculated (solid line) intensities based 
on model No. 5 for Dl with Dido 'Curran' at '♦.ZK. 

A. 29 (a) (b) Observed (crosses) and calculated (solid line) intensities based 
on model No. 5 for Dl with Pluto 'Panda' at k.2K. 

4.30 Observed (crosses) and calculated (solid line) intensities based on 

model No. 5 for Dl with Pluto 'Panda' at 4.2K with the molecule tilted 

at 5° to the (OOl) plane. 

m 

'j.31 Resolution function of the intensity transform (using a 'window' 

function) of the neutron powder diffraction pattern of Dl recorded on 
the 'Panda' at ii.2K. 

k.32 Resolution function of the intensity transform (using no 'window' 

function) of the neutron powder diffraction pattern of Dl recorded on 
the 'Panda' at 't.ZK. 

A. 33 Intensity transform of the neutron powder diffraction pattern of Dl 
recorded on the 'Panda' at ^.2K. 

^.3^ Intensity transform of a neutron powder diffraction pattern of Dl 

calculated on the basis of the model shown tn figure k.]S. 

A. 35 Calculated neutron pov/der diffraction pattern of the iodine sublattice. 



12. 

it. 36 Intensity transform of the calculated neutron powder diffraction pattern 
of the iodine sublattlce showing the first six Iodine atom nearest 
neighbour distances calculated from the lattice parameters of table ^.2. 

^.37 Neutron powder diffraction pattern of Dl recorded on 'Panda' at A.2K 
minus the calculated neutron powder diffraction pattern of the iodine 
sublattlce. 

A. 38 Intensity transform of the neutron powder diffraction pattern of Dl 
recorded on 'Panda' at 't.2K minus the calculated neutron powder 
diffraction pattern of the iodine sublattlce. 

A.39 Neutron powder diffraction pattern of Dl calculated on the basis of 
the model shown in figure ^.19 minus the calculated neutron powder 
diffraction pattern of the iodine sublattlce. 

k.kO Intensity transform of the neutron powder diffraction pattern of Dl 
calculated on the basis of the model shown in llgure ^.19 minus the 
calculated neutron powder diffraction pattern of the iodine sublattlce. 



Ch3p_t^er five^ 

5.1 Neutron powder diffraction pattern of Dl at 8AK (middle phase) 
recorded with the Dido 'Curran' diff Tactometer. 

5.2 Two-fold disordered model for Dl middle phase, space group Bba2 (No.^1). 

5.3 Four-fold disordered model for Dl middle phase, space group Bbcm (No. 64) 

S.k Proposed two-fold disordered model for Dl middle phase, space group 
Bbcm (No. 64) . 

5.5 Observed intensities of neutron pov;der peaks for Dl middle phase 
recorded at 84K with Dido 'Curran'. 

5.6 Calculated intensities of neutron povjder peaks for tv/o-fold disordered 
model of Dl middle phase space group Bbcm (No. 64). 



13. 



Chapter six 

6.1 Neutron powder diffraction pattern of Dl at 137K (high phase) recorded 
with the Dido 'Curran' dt ff Tactometer. • 

6.2 Zinc-sulphide model for Dl high phase, space group Fm3m (No. 225). 

6.3 Sodium-chloride model for Dl high phase, space group Fm3m (No. 225). 

6.'* Four-fold disordered model for Dl high phase, space group F¥3m (No. 216). 

6.5 Six-fold disordered model for Dl high phase, space group Fm3m (No. 225). 

6.6 Eight-fold disordered model for Dl high phase, space group Fm3 (No. 202) 

6.7 Twelve-fold disordered model for Dl high phase, space group Fm3m (No.225)| 

6.8 Freely rotating model for Dl high phase. 

6.9 Calculated structure factors for four-fold disordered model of Dl 
high phase, space group F^Bm (No. 216). 

6.10 Calculated structure factors for six-fold disordered model of Dl high 
phase, space group Fm3m (No. 225). 

6.11 Calculated structure factors for eight-fold disordered model of Dl 
high phase, space group, Fm3 (No. 202). 

6.12 Calculated structure factors for freely-rotating model of Dl high phase. 

6.13 Calculated structure factors for twelve-fold disordered model of Dl high 
phase, space group Fm3m (No. 225). 

6.1A Observed and normalized structure factors for Dl high phase recorded 
on Dido 'Curran' at 137K. 



]h. 



Chapter seve n 



7.1 Temperature dependence of the intensity ratio R(T) = ^— ; ; — =r- 

^ l(220) + l(22D 

measured over the lower phase transition of Dl with 
increasing (X) and decreasing (O) temperatures. 

7.2 Order parameter q{T) derived from neutron diffraction experiments (+) 
and specific heat measurements (solid line) over the lower phase 
transition region of Dl. 

7.3 Temperature dependence of the intensity ratio I (002)/l (1 1 1) measured 
over the upper phase transition of Dl wth increasing (x) and 
decreasing (O) temperatures. 

7.'» Order parameters q'(T) derived from neutron diffraction experiments 
(+) and specific heat measurements (solid and dotted lines) over the 
upper phase transition of Dl . 



Appendix two 

A2.1 Temperature measurement and control circuit for the variable tempera- 
ture cryostat used in neutron diffraction experiments. 

A2.2 Temperature gradient across sample, as a function of temperature, for 
the sample tube used for the lower phase transition neutron 
diffraction experiment on Dl. 



15. 



List of tables^ 

Chapter one^ 

1.1 Melting, boiling and phase transition temperatures for HI/DI. 

1.2 The lattice parameters of HI measured by X-ray diffraction. 

1.3 Activation energies for molecular re-orientation determined from 
dielectric constant measurements. 

].k Static dielectric constant and maximum absorption frequency vs 
temperature for HI. 

1.5 Rotational correlation times and activation energies for DCl, DBr and 
Dl as determined by NMR experiments. 

1.6 Correlation times for translat ional diffusion in HCl , HBr and HI cal- 
culated by random-walk theory. 

1.7 Observed and calculated second moments of NMR line shape for HCl, HBr 
and HI. 

1.8 Quadrupolar coupling constants measured by NQR experiments for DCl, 
DBr and Dl . 

1.9 Frequencies of infra-red absorption lines for HI (middle phase). 



Chapter three 

3.1 Lattice parameters of HI and Dl measured by neutron and X-ray powder 
diffraction. 



16. 



Chapter four 

J^.l Results of the lattice parameter refinement obtained with the program 
PARAM for the low phase of Dl at 'i,2K. The refinement is based upon 
the powder diffraction pattern recorded with the 'Panda' with a neutron 
wavelength of A = 1.886^. 

't.2 Miller indices of Bragg peaks observed in the neutron powder diffraction 
pattern of Dl in the orthorhombic and the six monoclinic unit cells. 

h.l Summary of the unit cells and space groups considered for the structure 
of the low phase of Dl. 

Chapter five 

5.1 Structural parameters of Dt at 84K (middle phase). 

Chapter six 

6.1 R-factors obtained in the least-squares refinement of four disordered 
models proposed for the high phase of Dl. 

6.2 Structural parameters of Dl at 137K (high phase). 



Chapter seven 

7.1 Intensities of neutron powder peaks measured over the lower phase 
transition of Dl with increasing and decreasing temperature and the 
R(T) and q(T) values derived from these measurements. 

7.2 Temperature dependence of the anomalous specific heat over the lower 
transition of Dl, The values listed were calculated from the data 
of Clcsius and Wolf (19^7). 



17. 



7-3 Intensities of neutron powder peaks measured over the upper phase 

transition of Dl with increasing and decreasing temperature and the 
ratio and q'(T) values derived from them. 



y.k Temperature dependence of the anomalous specific heat over the upper 
phase transition of Dl. The values in this table were calculated 
from the data of Clusius and Wolf (19^7)- 



18. 



CHAPTER ONE 



Review of previous work on HI and PI 



1. 1 Introduction 

This study of the crystal structures and phase transitions of the 
hydrogen and deuterium iodides was undertal<en as the last of a series of 
experiments on the hydrogen and deuterium halldes at Queen Mary College. 

Due to similarities in the calorimetric and dielectric properties of 
deuterium bromide and deuterium Iodide, It was expected to find that the 
two compounds had similar structures. It was therefore interesting that, 
although the middle and high phases of deuterium iodide and deuterium bro- 
mide were isomorphous, their low phases had quite different symmetries. 
The low phase of deuterium iodide was in fact monoclinic as opposed to the 
orthorhombic low phase of deuterium bromide. 

The hydrogen and deuterium iodides, hereafter referred to as HI and 
Dl, have similar physical properties which are summarized below In table 1.1 
Three of the hydrogen halides, HCl, HBr and HI, possess more than one solid 
phase which are stable within well-defined temperature ranges. This thesis 
describes the determination of the crystal structures of the three solid 
phases of Dl and the mechanisms of the two phase transitions between them. 



The experimental worl< performed on HI/DI prior to this thesis will now 
be described In detail. 



Compound 


Melting Point 


Boi 1 ing Point 


T (lower) 


T^(upper} 


HI 
Dl 


223. 2K 
221. 2K 


237. 6K 
237. 25K 


70. IK 
77. 3K 


125. 7K 
I28.3K 



Table 1.1 Melting, boiling and phase transition temperatures for 
HI and Dl 

T- is the temperature where the phase transition is complete, when observed 
with increasing temperature. 

1 .2 Calorimetric experiments 

The specific heat of HI at constant pressure C^was first measured 

f p 

by Eucken and Karwat (192^) over the temperature range of 58 to 237K. 
Subsequently as part of a series of very accurate experiments designed to 
test the validity of the third lav/ of thermodynamics, the specific heats 
of the hydrogen and deuterium iodides were remeasured first by Giauque and 
Wiebe (1929) (HI only) and then by Clusius and Wolf (I9'i7) (both HI and Dl). 
The results of Clusius and Wolf are shown in figure 1.1. The earlier 
results obtained by Eucken and Karwat were similar only the statistical 
scatter was greater. 



Both the C versus T graphs In figure 1.1 show two anomalies which are 
attributed to second-order phase transitions. A first-order phase transition 
is defined as a discontinuous change In a thermodynamic variable (such as 
C ) and Is usually associated with a change of structural symmetry. A 
second-order phase transition is associated with a change In a first- 
derivative of a thermodynamic variable. The two transitions in figure 1.1 
are of this type and are also referred to as lambda transitions due to the 
characteristic shape of the C versus T graph In the transition region. 
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col. 




200 



tCk) 

Figure^ 1,1 Specific heats of HI and DI at constant pressure 

as a function of temperature, Clusius and Wolf (l947). 
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The two X-trans It ions in HI and Dl are associated with the 
reorientation of the molecules in the crystal. By considering the changes 
In the molecular configuration when undergoing a transition between two 
solid phases, it is possible to calculate the anomalous entropy associated 
with the two second-order phase transitions. This will be discussed In a 
later chapter. 

The C versus T graphs of HI and Dl are similar to that of DBr but 
are unl il<e those of HF/DF, HCl/DCl and HBr. The different behaviour of 
HBr and DBr is an unusual isotope effect which was studied by Johnson (1970) 
The HCl/DCl transition is of the first order and was shown to be due to a 
change from an ordered orthorhombic structure to a disordered cubic phase 
by Farrow (1969), and Sa'ndor and Farrow (1369). HF/DF has no solid-solid 
phase transitions and its structure has been determined, Johnson and 
Sandor (1975). 

1 .3 X-ray experiments 

The X-ray powder diffraction pattern of solid HI was studied by 
Natta (1930 and 1931), by Ruhemann andSimon (1930 using photographic 
techniques and later by Mauer and Keffer using counter methods (196A). 
The experiments gave information only about the positions of the iodine 
atoms. 



(a) Natta In his first paper (1930) reported that HI was cubic at 103K 
(middle phase). For his second paper (1930 he obtained better 
quality photographs and was able to assign Miller indices to powder 
lines which he had previously considered as being due to possible 
impurities such as Ice or iodine (caused by dissociation of the HI 
samples during irradiation by X-rays). 



22. 



He reported the lattice parameters of a tetragonal face-centred 
lattice at lOOK and also described photographs taken at 188K, 112K, 
90K and 77K. These showed tetragonal symmetry but the lattice para- 
meters were not reported. He did not appear to have observed the 
cubic phase above 125K, nor the orthorhomblc nature of the middle 
phase, as was found In our experiments. Discrimination between the 
tetragonal and orthorhomblc symmetry may have been beyond the 
resolution of the X-ray photographs, which is not surprising as the 
orthorhomblc lattice parameters differ by only ]%. The absence of 
the cubic phase can only be Implied from his report as lattice para- 
meters are only quoted for a temperature of lOOK. 

(b) Ruhemann and Simon (1931) recorded the X-ray powder diffraction 
pattern of solid HI with a ^8mm diameter camera using CuKa radiation. 
Photographs were taken of the three solid phases at the following 
temperatures: 21K (low), 82K (middle) and 125K (high phase). All 
the powder patterns were Indexed on the basis of a face-centred 
tetragonal structure. The non-observation of the cubic high phase 
could probably be explained by a systematic error In the temperature 
measurement. 

(c) Mauer and Keffer (196^) studied X-ray diffraction patterns from a 

flat plate of HI powder by a counter di f f ractometer at six temperatures, 
They found that the structure of HI was face-centred tetragonal In 
the low and middle phases, and face-centred cubic in the high phase 
above 125K. 



A summary of the published X-ray work on Hi is given in Table 1.2. 
Some recent X-ray experiments on the middle and high phases of HI 
are described in chapter three (Simon 1970). 
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Temperature 



R 



c/a 



Phase 



Author 



103K 

77K 

90 K 
112K 
188K 
lOOK 

21K 

82K 

125K 

A.2K 
39K 
77K 
120K 

I3OK 
I8OK 



6.18 



t.io 

5.94±.01 
6.03±.01 
6.19±.01 

5.987 
6.000 
6, 0^12 
6.106 
6.2i*6 
6.295 



1.0 



1.08 

1.08 
1.08 
1.08 

1.082 

1.082 

1.080 

1.072 

1.0 

1.0 



fee 



fct 



fee 



Natta (1930) 



Natta (1931) 



Ruhemann and 
Simon (1931) 



Mauer and 
Keffer UsSk) 



Table 1.2 



The lattice parameters of HI measured by X-ray diffraction 



1 .k Optical experiments 

The optical isotropy of the hydrogen halldes was studied In 1938 by 
Kruis and Kalschew. They investigated the rotation of the plane of 
polarization of a beam of polarized light Incident on a 1mm thick sample of 
the crystalline hydrogen halide. The temperature of the sample was varied 
from hOK to 220K and the polarization of the transmitted beam observed. 



Their results indicated that HCl , HBr and HI were optically isotropic 
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In the high phase and anisotropic in the lower phases. The phase transition 
from the anisotropic lower phases to the isotropic high phase could be 
observed In each compound but transitions between the anisotropic lower 
phases of HBr and HI could not. These results are consistent with the high 
phase of each compound possessing cubic symmetry and the other phases lower 
symmetries. 

1 .5 Dielectric experiments 

(a) Static dielectric constant measurements 

The static dielectric constant of an isotropic medium e is defined 
by the equation: 



E 



E + 4ttP 



(1) 



where E_ Is the applied electric field, _D the displacement field and P_ is 
the polarization Induced Inside the material by the applied field. P is 
equal to the dipole moment per unit volume, the dipole moment of a sample 
of N atoms being defined as: 



p - I r q 

•'- *•, — n^n 
n=l 



(2) 



where r^ Is the position vector of the nth atom with charge q . The 

-1 8 

dipole moment is measured In Debye units (10 e.s.u.). This is approxP 

mately the charge of one electron multiplied by one AngstrUm. In S.I. 
units : 



ID = 3.336 X 10 
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The dielectric constant Is a dimenslonless number. The dipole moments 
of the hydrogen halide molecules are: 
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y(HF) = 

y(HCl) = 

y(HBr) = 

y{HI) = 



1.91 D 
1.12 D 
0.85 D 
0.^40 D 



These values were obtained from the temperature dependence of the 
static dielectric constant. The technique, using the Clausius-Mossott i 
and the Languin-Debye equations, is described by Kittel (1967a). 

The static dielectric constant of solid HI was first studied by Smyth 
and Hitchcocl< in 1933. A later series of experiments performed by Cole 
and Havrilial< in 1957 covered a wider temperature range and their results 
are reproduced in figure 1.2. The low phase section of figure 1.2 for 
HI/DI shows a low value for the dielectric constant Indicating that the 
structure is ordered. As the temperature rises, an increase In the 
dielectric constant is observed and can be attributed to the greater 
orientatlonal freedom of the molecules. The temperature range over which 
the transition occurs is approximately 20" K wide suggesting that the onset 
of molecular re-orlentatlon Is gradual. 



The subsequent decrease of the dielectric constant with increasing 
temperature suggests a disordered structure in the middle phase. The 
behaviour of the dielectric constants of HBr and DBr is very similar to 
HI/DI. HBr has an additional solid phase, mentioned in section 1.2, which 
is stable just below the upper phase transition point. The static 
dielectric constant of HBr at llSK can be seen to differ from DBr and have 
an additional minimum value in figure 1.2. In HCl the single phase 
transition is very sharp, and Is accentuated by a discontinuous rise in the 
static dielectric constant (see figure 1.2). 
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Figure 1,2 Static dielectric constants of HCl,HDr,DBr, HI and 
DI ae a function of temperature. Cole and Havrlliak (l957). 
Brown and Cole (1953) and Swenson and Cole (1954), 




Figure 1.3 Imaginary versus real dielectric constant of HI and DI 
below the lower phase transition point. Cole and Havriliak (1957) 
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Figure 1.4 Imaginary versus real dielectric consta-.t of III end DI 
above the lower phase transition point, Cole and Havriliak (1957), 
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The decrease of the static dielectric constant with temperature is 
similar to that of a disordered structure. The resultant dipole moment is 
reduced by the greater thermal motion of the molecules. This Debye-behaviour 
is described in liquids by the Langevin function as discussed by Kittel 
(1967b). 



L(x) = ctnh(x) - 1/x 



(3) 



kT 



, p is the dipole moment, E the electric field, k the Boltzmann 



constant and T the absolute temperature. 

Above the upper phase transition point, a sharp decrease in the static 
dielectric constant occurs which can be attributed to a decrease in density 
and/or an increase in or ientational disorder, possibly from two to three 
dimensions. This Is followed by a further continuous decrease due to the 
larger thermal motion as the sample approaches its melting point. 

(b) Complex dielectric constant measurerrents 

An oscillating external electric field may be applied to the material 
being studied and the molecular dipoles will attempt to re-orient In 
phase with this field. However, the intermolecular forces, thermal motion 
and the inertia of the molecules retard this re-orientation process. 

The average time taken for a molecule to return to equilibrium from a 
polarized state via a re-orientation process is defined as the relaxation 
time. To describe this behaviour mathematically the dielectric constant 

JL 

Is treated as a complex quantity, e". This complex dielectric constant 
has been measured in HI and Dl by Cole and co-workers. 



The frequency -dependent behaviour of the dielectric constant is Lest 
represented by the empirical formula of Cole and Cole (19^): 
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e* = e' - U" = (e -e )/[l + (iwi )^""] 

O 00 



(A) 



This is a modification of the theoretical formula for polar liquids 
derived by Debye (1929). 



e -e =(e-e)/(l+iuT) 

TO o " 



(5) 



where e is the static dielectric constant. e is referred to as the 

oo 

infinite-frequency value and represents the electronic polarizabi 1 i ty since 
at high frequencies the ionic and orlentat ional polarizabi li ties are small 
due to the inertia of the molecules. x is a generalized relaxation time, 
a is an empirical constant varying between and 1 and cj is the frequency. 



The complex dielectric constant appears as the arc of a circle with 
its centre on or below the real x-axIs when drawn in the complex plane of 
an Argand diagram. This empirical formula of Cole and Cole, eqn.(4), has 
been found to describe adequately the behaviour of the complex dielectric 
constant of a number of substances. Such diagrams are referred to in the 
literature as 'Cole-Cole' plots after the original paper of Cole and Cole 
(19^1). Each point on the arc represents a different frequency and the 
absorption maximum is at the point of maximum e''. 



Figures 1.3 and 1.^ show the complex dielectric constant of HI/DI 
at temperatures below and above the lower phase transition point, Cole and 
Havriliak (1957). The arcs drawn by Cole for figure 1.A were stated to 
be exactly semi-circular although only frequencies up to 500 Hz had been 
recorded. The subsidiary arcs were suggested by Cole as continuations to 
higher frequencies. Above the lower phase transition point the centre of 
the arc lies on the c'-axis, while below the transition point the centre of 
the arc lies below the e'-axis. To verify the semi-circular shape of the 
arcs in figure 1.^, Cole plotted logarithms of the frequency against a 
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function called the absorption conductance which is proportional to oje'' and 
is measured in Siemens or reciprocal ohms and obtained the lines of figure 
1.5. When tijT <<1 only low frequencies are considered and Debye's equation 
can be approximated to give: 



we 



CO 



(6) 



Cole plotted log^Q(we'') versus log^Q(^) and obtained a straight line 
of slope 2. This confirmed that a, the empirical constant of Cole and 
Cole (1941), was zero and that the frequency dependence of c , the complex 
dielectric constant had a Debye form of relaxation. When performing these 
experiments in 1957 Cole and Havrlliak were unable to extend the frequency 
range of their measurements beyond 1 MHz and were therefore unable to 
observe any dispersion in the high phases of HI and D! . 



Since the maximum absorption occurs at the resonance frequency, the 

relaxation time t v/hich is characteristic of this absorption is the 


reciprocal of the frequency. Cole plotted the logarithm of this relaxation 
time versus the reciprocal of the temperature and obtained' the curves of 
figure 1.6. A frequency term A related to the probability of the dipole 
undergoing a relaxation process was defined by the expression: 



A = ;p exp(E^/2.3RT) 




(7) 



where E is an activation energy and R Is proportional to the Boltzmann 
constant. Cole applied the relation: 



'°9io ^0 



log^gA + Eg/2.3RT 



(8) 



to these data and obtained a series of activation energies tabulated in 
table 1.3. 
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Figure 1.5 Absorption conductance of DI and DBr as a function of 
frequency, Cole and Havrlliak (l957). 
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Figure 1.6 Relaxation times for hydrogen halides as a function of 
temperature determined by dielectric absorption methods, Cole 
and Havriliak (l957). 




Figure 1,7 Imaginary versus real dielectric constant of HI as a 
function of temperature, Groeuewegen and Cole (l966). 
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Compound 


Temperature K 


log,oA 


E^ J mol"' 


HI 


68-62 


-13.6 


9,200 




80-70 


-17.2 ■ 


15,000 


Dl 


75-63 


-13.3 


9,200 




100-78 


-15.6 


13,000 



Table 1.3 



Activation energies for molecular re-orientation determined 
from dielectric constant measurements. 



The high frequency dielectric dispersion of HI was measured by 
Groenewegen and Cole in 1966. The results were represented by an empirical 
expression similar to (h) referred to as a skewed-arc function which is 
shown for their data in figure 1.7. 



e" - e = (e„-e: )/{! + itot ) 

oo o <» 



B 



< 6 < 1 



(9) 



Table 1.^ summarizes the values of the static dielectric constants e and 
the frequency of maximum absorption loss u for HI. 









Temperature K 


e D 




(!) HHz 
m 


78 


17.^ 


A. 6 


80 


15.9 


8.4 


Bk 


13.^ 


21.5 


88 


11.8 


tih 


92 


10.6 


86 



Table l.A Static dielectric constant and maximum absorption frequency 
versus temperature for HI 
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A value for e was obtained by extrapolating the low temperature data 

to high frequencies giving e = 2.9D. Using this value a normalized 

function ( — ;^-^) was plotted and is shown in figure 1.8. The value of the 

o ™ 
empirical constant 6, eqn.{9), defined by Groenewegen and Cole (1966) was 

0.67. 



In the discussion of their results, Groenewegen and Cole compare the 
behaviour of e for HI with that of simple polar liquids such as liquid 
glycols studied by Davidson and Cole (1951, i960 and 1961). They remark 
upon the similarity between the two results and that such behaviour has not 
been observed in any other solid. At low frequencies the dielectric con- 
stant obeys the simple Debye type of relaxation. However, the width of 
the absorption maximum is broader than that predicted by a single relaxation 
time. This difference in width persisted in all their experiments at high 
frequencies. 

■if 
The differences between the loci of e above and below the lower 

phase transition point are more noticeable at lower frequencies and are 

attributed to a persistence of dipole correlations over a longer time, 



I .e. 



t -d-a) 
^(t) = ( — ) below the transition point (10) 



(fi(t) <t exp( — ) above the transition point (11) 



((>(t) Is the dipole correlation function representing the ratio of the dipole 
moment at a time t to the dipole moment at a time zero. An average Is 
tal<en over the whole sample for the case of no applied electric field. 
(f(t) Is discussed by Kubo (1957), Glarun (i960) and Cole (1965). 
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Figure 1,8 Reduced locus of the Imaginary dielectric constant 
versus the real part, Groenewegen and Cole (l9£56). 
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Figure 1.9 Relaxation frequency as a function of temperature 
for HBr and HI, Groenewegen and Cole (l96e). 




Figure 1.10 NMR line width versus temperature for HI, Alpert (l949). 
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Glarum and Cole discuss the possibility of describing the relaxation 
time as a distribution of relaxation times. This is qualitatively similar 
to a diffusion process. Two theories which describe the behaviour of the 
dipole correlation function by the use of sums of exponent ional decay 
functions are discussed. Zwanzig (I963) analyses the dipole-d ipole coupling 
by a series expansion in terms of (1/kT) where k and T are the Boltzmann 
constant and absolute temperature. Though their results are in qualitative 
agreement wl th the experimental data their dipole correlation functions, 
eqn.(lO), decay too quickly to describe the observed relaxation times. 
This discrepancy is thought to arise because the energies of the dipole- 
dipole interactions in HI are larger than kT. This means that Zwanzig's 
series expansion is of doubtful validity. 

A second theory due to Adam (1965) attributes the reduction in the 
relaxation time to the presence of defects In the solid. This theory also 
agrees qualitatively with the experimental data but requires an excessively 
high concentration of defects for quantitative agreement. In addition 
the theory has difficulty in explaining the behaviour of the relaxation 
constant at the phase transition. 

In figure 1.9 the relaxation frequencies are plotted as a function of 
temperature. Their behaviour is considered to provide evidence for a 
co-operative effect at the phase transition. The frequencies are fitted 
by Groenewegen and Cole to an Arrhenlus rate law of the form: 



w = A exp(-AH /RT) 

m 



(12) 



ah" Is the enthalpy of the barrier to free rotation and A is the frequency 
at which the dipole crosses the barrier. 



For HI in the middle phase: 
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A - 1.3 X lo'^ Hz and H = 12,600 J. 

This is considered a high value for A, as 'stiff molecular vibration 
frequencies are of the order of 10 Hz. Groenewegen and Cole conclude 
that the correct interpretation of the data requires more sophisticated 
theories and better understanding of the crystal structures of the various 
sol id phases. 

1 .6 Nuclear magnetic resonance experiments (NMR) 

NMR experiments with solid Hi were first performed by Alpert {19^9). 
The first MMR worl<. had been reported by Purcell , Torrey and Pound and by 
Bloch, Hansen and Packard (19^6). The theory of NMR has been covered 
extensively by Abragam (1961) and a short introduction only will be given 
here. 

The nuclei of all atoms acquire an induced magnetic moment when placed 
in an external magnetic field. If a radio-frequency magnetic field Is 
introduced In a plane perpendicular to the applied steady magnetic field, 
the magnetic moment of the nucleus will precess about the axis of the 
uniform field. if this radio-frequency field has the correct frequency, 
a resonance wi 1 1 occur at a critical Larmor frequency v , where: 



y H 



11^ 



(13) 



H is the static magnetic field, y is the nuclear magnetic moment, 1 the 
nuclear angular momentum quantum number and T\ the Planck constant divided 
by 2-n. 



For diamagnetic materials, i.e. those with no intrinsic nuclear 
magnetic moment, the ratio of v to H Is independent of the nature of the 
material containing the nucleus which is however coupled to the material and 
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can lose energy to it. There are two main methods of energy transfer in 
the sample causing a line-width in the resonance absorption curve. 



a) 



b) 



Spin-lattice coupling or relaxation is the exchange of energy between 
the nuclear magnetic moment and the thermal motion of the lattice. 
Spin-spin relaxation is the exchange of energy between different 
nuclei by means of their magnetic moments. 



These relaxations are described by characteristic times; T, the 
spin-lattice relaxation time and Tj the spin-spin relaxation time. In 
general one may say that T^ depends upon the local magnetic field at the 
nucleus. it is therefore strongly dependent upon inhomogenei ties in the 
applied field H^. These act as a limit to the accuracy of line-width 
measurements in cases where large magnetic pole faces and field strengths 
are used. NMR is a useful technique for studying second-order phase 
transitions such as those observed in the hydrogen and deuterium halides. 



At the time of Alpert's wori< there were two theories of the phase 
transitions which will be briefly discussed here. 



a) 



b) 



Pauling (1930) and Fowler (1935 and 1936) advanced the theory 
that below the transition point the molecules were librating while 
above the transition point they rotated freely. The molecular 
rotation was assumed to be a co-operative phenomenon and a lambda- 
point transition was deduced from a statistical analysis. 

Frenkel (19^6) and Landau (1937) proposed an order-disorder 
transition mechanism similar to that observed In alloys. According 
to this view, below the lower phase transition point the molecules 
have fixed orientations (are ordered) while above the transition point 
they librate between equilibrium orientations (are disordered). 
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Alpert analysed re-orlentat ion effects in liquid norma 1-propa no 1 and 
n solid ice-l and obtained order-of-magnl tude values for the characteristic 
lorrelation or 'flipping' time x . These are 1.8 x 10 s for ice and 

s for propanol. These two times characterize re-orientation between 
squilibrium positions and free rotation respectively." 

For an (a) -type phase transition mechanism, the width of the_ resonance 
ine is expected to decrease rapidly on entering the higher phase due to the 
large change in correlation time. The lower limit of the line-width would 
te determined by the inhomogenei ty of the magnetic field. For a (b)-type 
:ransition a slight discontinuity is expected at the phase transition point 
lue to a change in crystal symmetry or density. This would be followed 
»y a decrease in line-width due to a gradual Increase in the frequency of 
(scillatlon as the temperature increased. Again a lower limit would 
►xist due to the field inhomogenei ty. The curve shown In figure 1.10 shows 
itrong evidence for a {b)-type transition mechanism until the field 
nhomogeneity limit is reached at 200K. These data suggest that the 
;orrelatIon times for the molecular re-orientation are of the order of 
) X 10 s falling to 5 x 10 s at the higher temperatures. Alpert observed 
[a)-type transitions in CHr and CH-D at their upper phase transitions, indlca- 
:ing free-rotation. 



No further NMR experiments were performed on HI until 1968 when 
ienin and co-workers studied HI/DI over a temperature range of 50-200K and 
jbtained values for the four parameters: 

T. the spin-lattice relaxation time, 

T^ the spin-spin relaxation time, 

T, the rotating frame relaxation time, 
Ip 

jnd M^ the second-moment of the polycrysta 11 ine NMR curve profile. 
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Two new parameters are Introduced here; T. , the rotat i ng-f rar 
elaxation time, is a useful experimental quantity measured by an Interesting 
echnique. 



In the compound under study a pulsed magnetization is produced at 90° 
o the static magnetic field H^ of eqn.(I3) and Is maintained by a similar 
ransverse magnetic field which Is rotated at the Larmor frequency. The 
pin-spin coupling Is nullified by rotating all the nuclear magnetic moments 
n phase. However, the energies involved are much lower than for T, , or 
pin~Iattice, type relaxation processes and thus the frequency range of 
e-orlentations observed is In the KHz range instead of MHz. Essentially, 
5 the name implies, the spin-lattice decay time T, is measured In a new 
rame of reference which is rotating at the Larmor frequency. M- Is the 
Econd moment of the resonance curve profile f(oj) measured in a polycrystal 1 ine 
ample. In general the moments are defined as: 



(u-o) ) f(w)d(, 



Abragam (1961) 



(1^) 



It is a useful function indicating the width of the resonance curve 
lich can also be calculated from the Hamiltonlan of a simple spin system. 



Figure 1.11 shows the temperature dependence of T, in the high phases 

f Dl , DBr and DCl. Table 1.5 lists the rotational correlation times x 

c 

id the activation energies AE in the high phase. The tabulated values 
ire derived from the results plotted in figure 1.11 using the formulae: 



-^ { ^ } -i- 



(15) 



ind 






(16) 



1. Is the quadrupole moment and the other parameters and constants have their 
jsual meanings. 
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Figure 1.11 Spin-lattice relaxation times versus temperature for 
DCl,DBr and DI, Genln et al, (l968). 
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Figure 1,12 Spin-lattice O , rotating-f rame j^ and spin-spin Q 
relaxation times and second moment data Q for tlCl,HBr and 
HI as a function of temperature, Genin et al. (l968). 
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Compound 


T at upper transition 


AE 


00 


DCl 
DBr 
Dl 


2 X 10~'^s 

-12 
1.5 X 10 s 

1 .'i X 10 


■3560 J 
4560 J 
4060 J 


3 X 10"'^5 

2 X I0"^^s 

-14 

3 X 10 '^s 



Table 1.5 Rotational correlation times and activation energies for DC!, 
DBr and D) as determined by NHR experiments. 



Genfn attributes this rapid rate of re-orientat Ion , which is nearly 



.8 



10 MHz, to the 'twelve-fold disordered' structure of the high phase pro- 
posed by Sa'ndor and Farrow (1967 and 1969). This is statistically a face- 
centred cubic structure in which the Dl molecule points randomly towards each 
of its twelve nearest-neighbours. The orientation of the molecule changes 
by the deuterium atom moving to different equilibrium positions. The iodine 
atoms remain fixed in their face-centred positions. The model will be dis- 
cussed in a later chapter. 

The decrease in T] with increase in temperature, which occurs in the 
high phase of Dl and DBr near the melting point, was attributed by Genin to 
translatlonal diffusion of the molecules. An explanation of the inter- 
pretation of resonance line narrowing In terms of diffusion is given by 
Kitte! (1967c). Evidence of diffusion was also provided by the minima 
shown in T, in figure 1.12. Using the random-walk theory of diffusion and 
the method of Torrey (1953) and Resing and Torrey (1963), Genin was able to 
calculate correlation times for translatlonal diffusion which are listed in 
table 1.6. As can be seen, the correlation times for translatlonal 
diffusion ore of the same order of ragnitude as for rotational re-or i entdt icn . 
However the calculated activation energies are larger as would be expected 
In a crystal lattice. 
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Compound 


Energy 


00 


HCl 

HBr 
HI 


16300 J 

25900 J 
25500 J 


6 X 10~'^s 

1 X lO'^'^s 
8 X 10^^ 5 



Table 1,6. Correlation times for translational diffusion in 
HCl, HBr and HI calculated by random-walk theory. 

Genin also observed an anomalous curvature of the T, versus tempera- 

Ip ^ 

ture graph in figure 1.12. He attributed this curvature to quadrupolar 

interaction between the iodine nuclei. A minimum value of T, due to this 

Ip 

interaction was calculated assuming that slot; translational diffusion was at 

a minimum. A value of 1.9 x 10 s was obtained for T, which is in agree- 

Ip ^ 

ment with the data of figure 1.12. The anomalous curvature was not observed 
In HCl because the quadrupole moment was much smaller. 



Values of M„, the second moment, are listed in table 1.7. These were 
calculated by considering the hydrogen atom to be at the iodine site, 
I.e. at the time average of the twelve disordered positions. 



Compound 


M2 obs. 


M2 calc. 


HCl 


11 X 10"^s"^ 


9.2 X lO'^s'^ 


HBr 


8 X lO'^s"^ 


7.6 X 10"^5"2 


HI 


4 X 10"^s"2 


5.3 X lO'^s"^ 



Table 1.7. Observed and calculated second moments of NMR 
line shape for HCl, HBr and HI. 
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Another minimum occurs in T. in the upper region of the middle phase 
(see figure 1.12), Genin attributes this minimum to a 'three-plane' 
re-orientation process. He assumes a slow re-orientation between adjacent 
sites in the x-y plane, with fast re-orientations between sites in the x-z 
and the y-z plane. By 'fast' are meant re-orientation processes occurring 
at frequencies greater than uj, , the re-orientation frequency in the rotating 
frame system. The sites referred to are in the planes (101), (lOl), (loT) , 
(ToT) , etc. 



In the lower region of the middle phase of HI, a minimum is observed 
in T, near the phase transition point (see figure 1.12). Genin attributes 
this to a 'two-plane' re-orientation process in the x-z and y-z planes. 
In a 180° re-orientation the intramolecular dipole interaction is unchanged 
unlil<e the intermolecul ar interaction which is different. Thus in Dl the 
relaxation times for a 180*" re-ortentat ion would be different from those in 
HI. Expressions for the change In relaxation time have been derived by 
O'Reilly and Tsang (1967). For a 2:1 mixture of DI:HI, the values of T, 
and T. were only slightly shifted and it was concluded that 90° re-orientations 
alone are important. Genin calculated a value for T, for this process and 
claimed good agreement with his observed data. 



In the lov; phase of HI a single-plane re-orientation was proposed 
between the (IIO) and (llO) planes. T. was calculated for this process 
and agreement claimed with the observed data. The correlation times calcu- 
lated from these data are shown in figure 1.13 together with the correlation 
times obtained from dielectric constant measurements by Cole and Havriliak 
(1957) and Groenewegen and Cole (1966). As mentioned earlier, the low 
frequency effects are determined from the T. measurements and the high 
frequency effects from the T, and dielectric constant measurements. 
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Figure 1.13 Relaxation times as a function of temperature 
derived from NMU and dielectric loss measurements on HI, 
Genln et al . (1968) and Groenewegen and Cole (1966), 




Figure 1.14 Nr-Ol dispersion signal In polycrystalline DI at 
96K and 7MHz, Genln et al. (19G8). 
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Figure 1.15 NQR coupling constants of polycryatalllne DCl, 

DBr and DI as a function of tempc^rature, Genin et al.(l968), 



Genln also studied polycrystal 1 ine DI by a technique known as 
continuous-wave dispersion which has been described by O'Reilly (1958). A 
polycrystalline sample is placed In a constant magnetic field which Is 
modulated by an audio-frequency magnetic field. Nuclei which possess a 
nuclear quadrupole moment, where the magnetic quantum number is >_ h exhibit 
peaks in the modulated signal. These peaks result from the different 
energy levels available to the nucleus in the local field and are an average 
over all orientations of the powder sample. 

DI was studied in the low and middle phases at 7 MHz. Figure l.U 
shows an observed dispersion signal for DI at 96K. The quadrupolar splitting 
i and the coupling constant C are connected by the equation: 



^ 3 ^ h 



(17) 



where e is the electronic charge, q is the field gradient tensor and Q. is 
the nuclear quadrupole moment. 

Table 1.8 lists the quadrupolar coupling constants C derived from the 
powder NdR spectra of Genin at 50K. Figure 1.15 shows the variation of C 
with temperature. By using his one- and two-plane re-orientation processes. 
Genin predicted minima In the value of C at the same temperatures as the 
minima in T, in figure 1.12. These minima are not seen in figure 1.15 and ^ 
their absence is explained by Genin on the grounds that the 90° re-orientation 
sites do not possess equal occupational probabilities. 



Genin concludes by stating that the relaxation times derived from 
magnetic resonance experiments and dielectric-loss measurements agree to an 
order of magnitude. For free rotation they would be expected to differ by 
a factor of three. The results of the NHR measurements show cleariy that 
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the onset of the re-orientation processes coincides with the phase 
transitions. 



Compound 


h 


DC I 


190 


KHz 


DBr 


153 


KHz 


Dl 


113 


KHz 



Table 1.8. Quadrupolar coupling constants measured 
by NQR experiments for DCl , DBr and Dl . 



1.7 Infra-red spectroscopy 

Infra-red spectroscopy is an experimental technique which determines 
either: the complex dielectric constant e by reflecting infra-red 
radiation from a crystal or the absorption coefficient which is measured as 
a function of frequency by transmission through a thin film or plate. 

The vibrational motion of a molecule in a crystal may be described by 
a system of normal modes. The frequencies and intensities of the absorption 
lines are characteristic of these modes. By the use of group theory on 
the frequencies of the lines the data can be interpreted to provide a possible 
site or point group symmetry for the crystal. Review papers on this 
technique have been written by Bhagavantam and Venkataryuda (1939), Halford 
(19A6) and Winston and Halford (19^9). 

The infra-red spectrum of solid HI was first studied by Zunino (1936) 
and later by Hornig and Osberg (1955). The infra-red reflectivity for 
molecular crystals such as the hydrogen halides is generally weak so that the 
absorption coefficient is usually measured. 
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Hornig and Csbcrg recorded the infra-red spectrum of HI at 83K and then 
continued to cool their sample by evaporating liquid nitrogen at a reduced 
pressure. Figure I.I6 shows the transmission spectra after 30 minutes and 
1 hour, the final spectrum being recorded at 66K. ' The type of spectrum 
cliange observed by Hornig and Osberg depended on the thicl<ness of the film. 
No change was observed In the infra-red spectrum of the thickest films due 
probably to incomplete transition. Although the frequency of the main 
absorption line in HI was temperature independent, additional smaller sharp 
peal<s appeared In the spectrum below the lower phase transition point (see 
figure 1.16). Table 1.9 summarizes Hornig's results. 



Frequency cm 


THz 


Coinnents 


2120 
2130 
2110 
2075 
2180 


63.5 
63.8 
63.2 
62.2 
65.3 


Main absorption 1 ine 

) Subsidiaries present In 

) , 

) the low phase 

Weak band 

Very weak, observed only 
in thick f I Im band 



Table 1.9. Frequencies of infra-red absorption lines for HI 
(middle phase) 



The last two bands are probably sum and difference bands of a 50 cm 
lattice vibration frequency, coupled with the main absorption line. 
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Hornig discusses HI very briefly. He appears to have misinterpreted 
the data of Kruls and Kaischew (1938) as to the optical isotropy of HI. 
Apart from this he suggests two possible structures for the low phase. 
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Figure 1,16 Infra-red transmission spectra of HI as a function 
of temperature on cooling throuRh the lower phase transition 
at 70K, Mornie and Osberg (l955). 




Figure 1,17 

Infra-red transmission spectra of HI 
as a function of temperature^ Htebert 
and Hornlg (l957). 
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{a) Parallel or anti-parallel chains along the tetragonal axis, 
(b) Sets of four molecules arranged perpendicular to the tetragonal 
axis. (it is not made clear In what arrangement). 



The main component of the Infra-red absorption spectrum is associated wit^ 
the frequency of the Intramolecular stretching vibration. As part of a 
set of calculations of force constants for HC 1 and HBr, Hornig touches 
briefly upon HI . 
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is the potential energy change when the crystal is distorted by the i th and 

k I 

jth co-ordinates r. and r. In the kth and Jlth unit cells, then for solid 



HI; 



for HI as a gas ; 



f^^ = 266 Nra 
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where f., is the diagonal element of the force-constant matrix. The 
difference between the gas and solid phase values of f,, is taken by Hornig 
as an indication of considerable molecular interaction in solid HI. 

Hiebert and Hornig in 1957 recorded the Infra-red absorption spectra of 
H! using a liquid helium cooled infra-red spectrometer. The infra-red 
spectra are shown in figure 1.17. The following conclusions were drawn; 

(a) Between 25K and 70K, HI crystallizes in an ordered structure 
different from that of HCl and HBr. 

(b) Below 25K, HI undergoes a transition to a phase of lower 
symmetry with possibly more molecules in the unit cell. 
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(c) This lowest phase transition Is reversible, shows no hysteresis 
and has a negligible heat of transition. 

in 1964 Reeves, Robinson and Taimsalu published a conference abstract 
on the vibrational spectra of five phases of solid Hi. The lowest phase 
was reported to have a vibrational pattern similar to HCl and HBr suggesting 
Isomorphous structures. 



In a thesis abstract in 196A, Reeves reported that HI has a helical 
hydrogen-bonded structure with chains randomly directed parallel or anti- 
parallel to the tetragonal axis of the crystal. This phase was stable 
below lOK. 
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Savoie and Anderson (1966) attempted to clarify the position in a 
footnote to their paper on HCl/DCl and HBr/DBr. In a reference to Reeves 
et al they point out that the early experiments, described in sections 1.2 
to 1.6 of this chapter, and some unpublished Raman spectroscopy experiments 
on HI disagree with Reeves' results. Reeves' results are attributed to a 
metastable phase such as have been observed in solid argon and ammonia. 
The method of depositing thin films from the vapour phase onto a cold sub- 
strate was considered to favour the occurrence of such metastable phases. 
The three solid phases of Dl observed in the present neutron diffraction 
studies are in agreement with the results of Giauque, Cole, Alpert, Genin 
etc. who all studied bulk samples and observed three solid phases. 
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1.6 Theories of phase transitions in molecular crystals 

A detailed review of the different theories of phase transitions in 
molecular crystals has been given by Johnson (1970). The work of the 
following authors has been reviewed; Pauling (1930), Chang (1937), 
Kirkwood (19A0), TIsza (19^7). Bragg and Williams (193^), Powles (1952) and 
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Krleger and James (195^) • Further relevant papers have been published 
since this review; Hanamura (1970), KobayashI et al. (1969) and 
Hanamura (1970). 

Kobayashi et a! propose a two-dimensional mathematical model of a 
hydrogen hallde using two sub-lattices which are the low phase structures 
of Del and DBr proposed by Sdndor and Farrow (1967) and S^ndor and Johnson 
(1968). Order parameters are derived in terms of the dtpole moments of 
the sub-lattices using the consistency relations and Internal field 
approximation of Krieger and James (195^)- These equations are solved by an 
expansion method enabling the internal energy to be calculated as a 
function of temperature. The force constant used in the mathematical 
model is allowed to have three different ranges of value. The calculated 
internal energy predicts the orders of the lower phase transitions of each 
compound. HCl lies in one range giving a first-order transition and 
HBr/HI in another giving a second-order transition. (HI was assumed to 
be Isomorphous with HBr in the low phase). No numerical values of transition 
temperatures were computed. 
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Hanamura summarizes the published structural data on the hydrogen halldes |^ 

'■■i 
and concludes that 90° zig-zag bonding Is most likely. Dipolar and jji 

quadrupolar bonding energies are tabulated for each hydrogen halide. The 

upper phase transition In HBr and HI is said to be due to ordering of chains 

of molecules by quadrupolar interactions, the lower is due to dipolar 

Interactions. 



The Internal field approximation is used to calculate potentials betvjeen 
the chains. Three zig-zag models are considered and their dipolar energies 
calculated. For the low phase, the model corresponding to the proposed 
structure of Sender and Johnson (I968) has the lowest energy. If only 
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quadrupolar coupling is retained, the middle phase of HBr is generated as 
a random mixture of two of the three structural models considered by 
Hanamura. 
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CHAPTER TWO 



Experimental work 



The experimental work is divided into three sections; the preparation 
of samples of HI and 01 of high chemical and isotopic purity, the X-ray 
powder diffraction study of Hl/DI and the neutron powder diffraction study 
of Dl. The neutron diffraction work was carried out at the Atomic Energy 
Research Establishment, Harwell, using powder dl ff ractometers attached to 
the research reactors "Dido" and "Pluto". 

2. 1 Chemical preparations 

In both preparations Pyrex glass systems were used with ground-glass 
vacuum joints, hydrogen-free PTFE vacuum grease and greaseless valves made 
by Springams Ltd. The glassware was carefully cleaned with detergents, 
chromic acid and distilled water In turn. The dry components were then 
assembled, evacuated and maintained at a pressure of approximately 10 mm 
of Mercury for several days. During this period the glassware was heated 
by a hot-air blower to remove any residual internal film of water vapour, 
in all cases where chemicals were weighed, measured or placed into flasks, 
the operations were carried out in a glove-box filled with dry nitrogen. 

The reagents were purchased from the Koch Light Chemical Company with 
the following stated purities. 

Iodine 99.998^ 

Heavy Water 99.7^ 

Red phosphorus 99-999^ 



where ordinary water was used it was first doubly distilled in a conventional 
Liebig system. 
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In the preparation of Dl it was essential to keep the amount of HI 
impurity to a minimum as the incoherent neutron scattering cross-section of 
hydrogen is much larger than those of deuterium or iodine. For this reason 
Dl was at no time allowed to come into contact with any hydrogenous material 
with which isotopic exchange reactions could occur. 

(a) Indirect preparation of HI/DI 

This method used the reaction of red phosphorus and water on iodine; 



P + 51 + ^H^O = 5HI + H^PO^ 



(19) 



The iodine and red phosphorus were weighed, then ground together with 
a pestle and mortar in a glove box and the mixture was poured into flask A 
(figure 2.1). The preparation system shown in figure 2.1 was evacuated 
and sealed at V. Water or heavy water was allowed to run into flask A 
through a dropping funnel B. 



The Hl/DI evolving tn flask A was then passed through cold traps at 
C (cooled to -yO^C with a mixture of solid carbon dioxide and methanol) to 
remove any traces of iodine, water, phosphine (PH^) or phosphonium iodide 
(PHlI). The formation of the last two compounds was kept to a minimum by 
ensuring an excess of iodine in the initial reaction. The Hl/Dt was con- 
densed in flask D which v^as cooled with liquid nitrogen. The reaction, 
which proceeded steadily at room temperature, produced approximately 20ml 
of solid HI/Dl in 12 hours. The amount of water or heavy water running 
Into flask A from the funnel B was regulated in such a way that the pressure 
of the evolving gas was kept below 300mm of mercury, as recorded by the small 
aneroid vacuum gauge G. The rate of reaction could be slowed down by cooling 
flask A with liquid nitrogen. 



After 12 hours the flask D was sealed off and flask A replaced by an 




Figure 2.1 HI and DI sample preparation system used for the indirect method 2.1a 
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empty flask. The system was re-evacuafed and the HI/D1 transferred from 
to A through the cold traps several times and finally condensed in the 
1-metre storage tube E, which was subsequently stored in a 25-litre aluminium 
liquid-nitrogen dewar made by Spembley Ltd. Thewhole preparation was 
carried out In a photographic dark room with a weak "safe-light" because of 
the light sensitivity of Hl/DI. The sample produced was a pure white, 
polycrystall ine solid. 

(b) Pi rect preparation of PI 

To test whether the method of sample preparation had any influence on 
the neutron powder diffraction pattern, a sample of Dl was prepared by a 
direct method of catalytic combination. Iodine and gaseous deuterium 
were combined (at J400°C) in the presence of platinized asbestos as a catalyst 
In the apparatus shown in figure 2,2. This direct method used the reaction; 



D2 + l2 = 2DI (at ^00°C) 



(20) 



The rate of this reaction is very slow, less than 2ml of solid D! 
being produced per day. However the sample was expected to be of higher 
purity than the one prepared by the previous method due to the absence of 
any phosphorus compounds. The neutron diffraction patterns recorded for 
the two samples of 01 prepared by methods (a) and (b) were identical within 
statistical error. 

The electrolysis plant D of figure 2.2 was filled with heavy water and 
OjSO. in a dry glove-box. The 5-litre flask A and the catalyst were heated 
for a day at ^OCC to remove any grease or contamination. The Iodine was 
introduced into flask A In a dry glove-box. The system was evacuated and 
the 1-litre flask P was filled to a pressure of 1 atmosphere of deuterium 
at room temperature as recorded on the pressure gauge G. The deuteriufn was 
produced by electrolysis of the mixture of heavy v/ater and 02^*^1^ using a 
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Figure 2.2 DI sample preparation system used for the direct method 2.1b. 
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current of 5A, Ihe electrodes E and F being cooled by an external water 
supply W. The oxygen produced was removed by a rotary pump via a needle 
valve at V. 

The deuterium was dried by letting It pass through a condenser C at 
liquid nitrogen temperature. Any oxygen present was either condensed or 
converted to D„0 by heating with a hot air blower over some platinum cata- 
lyst in P. The 5-litre flask A, which was initially evacuated and sealed 
off by a valve, was opened to the remainder of the system so that the 
1-atmosphere of deuterium in the 1-litre f1asl< P now filled a volume of 
6 litres. The 5-litre flask containing iodine and deuterium at the 
resultant pressure of l/6th of an atmosphere was resealed by the valve and 
wrapped in an asbestos blanket and sealed completely by melting the glass 
at M. This was somewhat hazardous due to the risk of an explosion being 
caused by the glass-blowing torch. 



Flask A containing iodine and deuterium was heated for 12 hours at 

i<0O°C in an oven and after cooling was connected to the system shown in 
figure 2.3. The seal at B was broken by lifting the steel bar S with a 

magnet. The Dl produced in the flask was condensed in the evacuated 
system at D and purified in the same manner as described in section (a) 
above by cold traps at C. This latter operation was again performed in 
darkness. Both of these preparations are also described in less detail by 
Brauer (1963). 



2.2 X-ray diffraction study of Pi 

(a) X-ray diffraction photographs were taken of Dl at Queen Mary College 

with a Nonius and a Unicara Weissenberg goniometer using them as low- 

tcrnperature powder cameras. 
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Figure 2.3 DI sample purlf icatlon system used for the direct method 2.1b 




CO 



fftiM^m*M ,. ttut'ir-. 



- »^ * r ^— „- --:n^& j^-^ -:.^^. 



59. 



Samples of Dl were sealed into Linclemann glass capillaries of wall 
tliickness .Olmm and diameter .Imtn. The sealed capillaries were cooled 
by a continuous flow of cold nitrogen gas produced by boiling liquid nitro- 
gen in a pressurized glass dewar with an electrical lieater. The temperature 
of the sample was recorded by a copper-cons tantan thermocouple attached to 
the capillary and the thermoelectric emf was fed back to control a heater 
in the dewar arm as shown in figure 2.h. The temperature control system 
used in this experiment is described in detail by Farrow (I969) and 
Johnson (1970). The lowest temperature attainable with this equipment was 
approximately 83K. This limited the use of the cameras to the study of 
the middle and high phases of HI and Dl. 



The X-ray beam was monochromated by a point-focusing monochromator 
using a doubly-curved lithium fluoride single crystal plate, cleaved along 
the (001) planes. 



(b) X-ray powder diffraction patterns were recorded at ^.2K at Imperial 
College, London using a Siemens powder dl f f ractometer with a small liquid 
helium cryostat which retained liquid helium for two to three hours. 
The di ff ractometer and cryostat has been described by King (1967). A flat- 
plate reflection geometry was used which had two distinct disadvantages; 

(i) The flat Dl powder sample had to be supported by a transparent 
film mounted in front of it. A "Melinex" film of thickness 
0.001 in. was used for this purpose which gave a very Intense, 
broad diffraction peak which made it very difficult to observe 
the first dozen peaks in Dl. 

(ii) More inportant was the uncertainty of whether the sample was 
or was not poly-crystal 1 Ine. Since the powder was condensed 




SAMPLE- WITH 
ATTACHED THERMOCOJPLE 



AVER- LINE SCREENS 
FILM CASSETTE 



ON 

o 



Figure 2,4 Low temperature x-ray goniometer used at Queen Mary College, 
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directly from the vapour phase, the size and orientation of 
the crystallites could not be altered. Previous experiments, 
using the Nonius camera as described above, had shown that 
on average only lOS of the samples condensed were completely 
polycrystal 1 ine. The diffraction photographs contained 
powder lines with higher intensity spots due to larger 
crystallites within the sample. Because the Siemens diffracto- 
meter recorded an integrated intensity from the whole surface of 
the sample plate it was important to know If the sample were com- 
pletely polycrystal 1 Ine. A reflection photograph was unsatis- 
factory due to focussing problems and a transmission photograph 
was not possible as the sample and container were 5mm thick. 



(c) To overcome these difficulties a low- temperature X-ray powder 
di f f ractometer was assembled at Queen Mary College Incorporating the 
following features. 

(i) A large, well-tested cryostat capable of retaining liquid helium 
for 50 hours. 

(li) Cylindrical sample geometry, using thin-walled Lindemann glass 
capl 1 laries. 

(iti) An X-ray proportional detector either scanning continuously 
with a chart-recorder output or counting at fixed angular 
intervals and printing the recorded Intensities on a teleprinter 
enabling weak peaks to be observed. 

(d) These three features will be described In more detail below. 

( i ) The cryosta t 

A British Oxygen Company modular research cryostat was used to keep 
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the sample at 4.2K. The liquid helium capacity was 3 litres and it 
retained liquid helium for approximately 50 hours. The cryostat is shown 
in figure 2.5- The incident X-ray beam is shown at B passing through 
three cylinders surrounding the sample S (l<nown as 'tails'). The tails 
were made of solid beryllium ingots machined into tubes with a window 
thickness of 0.005 in. The three tails serve as: 

(a) an outer vacuum jacket, 

(b) a radiation shield maintained at the temperature of liquid 
nl trogen, 

(c) a radiation and vacuum shield at liquid helium temperature. 

These components were made at A.E.R.E., Harwell and A.W.R.E., Aldermaston. 

Due to the thinness of the tails the cryostat was supported in an 
aluminium alloy stand with a large section machined away leaving a minimum 
amount of material in the area of the diffracted beam. Aluminium was 
used to avoid X-ray fluorescence caused by the copper radiation. The 
sample space was thermally insulated from the liquid helium container LH by 
an evacuated exchange space ES . The sample of Dl was condensed as a solid 
into a Lindemann capillary tube using liquid nitrogen and an ancillary 
glass apparatus (Johnson 1970). The capillary was sealed by melting at 
each end and was then mounted on the tip of a long thin stainless steel 
tube of low thermal conductivity. Thermocouples were attached to the 
capillary with varnish and a low power heater (lOmW) of fine gauge wire 
was wound onto the capillary above the sample region. The capillary was 
dipped quickly into liquid nitrogen outside the cryostat to form a poly- 
crystalline sample of solid Dl. The long tube and sample were next 
transferred to the exchange space ES while liquid nitrogen was poured 
slowly over the capillary to prevent the sample melting and annealing into 
a powder with preferred orientation. 
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Flguru 2,5 Liquid helium cryostat used on the Queen Mary College 

x-ray powder dif f rnctometer. 
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The exchange space ES was then evacuated at V. The sample v/as aligned 
in the beam by means of two transverse slides T and three vertical tilting 
bolts H. The sample was first aligned in the X-ray beam optically using 
a fluorescent zinc sulphide screen. Final alignment was completed using 
short exposure X-ray transmission photographs. The position of the 
capillary could easily be detected due to absorption in the thermocouple 
leads and sample. Electrical leads to the sample passed through a neoprene 
vacuum seal at E (Batchelder and SIdey I969). 

(ii) To test the nature of the sample a short exposure X-ray diffraction 
photograph was taken at low angles to study the (ill) and (002) powder 
peal<s. After several attempts it was usually possible to form a good 
powder sample with continuous powder lines. The sample temperature was 
monitored by means of a gold-iron versus chromel thermocouple attached to 
the capillary with an Ice reference junction in a glass dewar. The 
temperature was controlled by means of an Oxfprd Instruments 'Harwell' 
temperature controller with the low power lOmW heater (described previously) 
attached to the capillary. 



(til) The electronic control system of the dl ff Tactometer was constructed 
of 'Harwell 2000 series' electronics. The X-rays were detected by a 
proportional counter filled with xenon. The amount of "white" radiation 
recorded by the counting system was reduced by the use of an amplifier and 
a pulse-height analyser on the signal from the proportional counter. 
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The instrument could either be used as a continuously scanning 
counter with a rate-meter and chart recorder output, or with a preset step 
length and count time in which case the output was printed on an ASR33 

teleprinter and also punched on a paper tape. The scattered Intensities 
were then drawn on a computer-controlled graph plotter. 



65. 



The X-ray powder di f f ractometer was made by Maclean Engineering Ltd., 
and was a less sophisticated version of the Harwell 'Apex' goriometer 
(Baker et al. 1966). 

2 . 3 Neutron diffraction study of DI 

The neutron diffraction experiments were carried out with three 
different powder di ff ractometers (Dido 'Curran', Dido 'Badger II' and Pluto 
'Panda') at the 25 MW thermal reactors Dido and Pluto (SRC, 1973). The 
three dl ff ractometers have the following features in common. 

- (i) The use of boron trifluoride proportional detectors. The boron 

is enriched in the isotope B which has a high neutron absorption 
cross-section and undergoes a nuclear capture reaction yielding 
a pulse-producing alpha particle according to the expression, 



B'°(n,a) Li^ 



(21) 



(ii) The counting time is controlled by a monitor counter. Counting 
for fixed times at small angular increments is unreliable due to 
fluctuations in reactor power, particularly in the last week of 
a four week reactor cycle when spent fuel elements are removed 
so changing the neutron flux at certain beam holes. To overcome 
this difficulty a small fission detector is placed in the incident 
monochromatic beam in front of the sample and acts as a monitor 
while removing less than 1^ of the incident beam. 

(ill) Automatic data printing. The use of this technique is necessary 
because the counting rate is so low that the recording of a powder 
pattern can take several days. The reactor operates for 2A days 
before ref ue 1 1 ing. 
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( i v) The use of 'Boral' (boron carbide compressed between aluminium 
sheets) shielding. Boral is used to define the incoming 
monochromatic neutron beam and the collimator is adjusted by 
taking neutron radiographs of the sample. with a polaroid camera 
fitted with a 1 i th i um- loaded zinc sulphide scintillator. Small, 
highly absorbing cadmium strips are fixed to the top and bottom 
of the sample container to define the sample position in the 
radiograph. 
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The next section will describe each di f f ractometer in more detail 



(^) Dido 'Curran' pov/der di ff ractometer 

The bulk of the neutron diffraction work was performed using this 
instrument of which the experimental layout is shown in figures 2,6 and 2.7- 
There are five counters CI to C5 fixed at 10* intervals relative to each 
other in a shielding block SH. The block can move on a track T and is 
supported on a cushion of compressed air, the total angular range covered 
by the first counter (Cl) being from -25** to 35° in 29. The sample was 
at S and could be rotated at half the angular velocity of the counter 
block. 



Monochromatic neutrons of wavelength 1.06A were provided by a lead 
monochromator crystal M which reflected neutrons at an angle of 22° from 
the (ill) lattice planes. By rotating the sample at half the angular 
velocity of the counter bank and using the full angular range available it 
was possible to record the diffraction pattern from d = y^S" to 17-5° wfth 
five counters from five sets of crystallites. Other regions were covered 
with fev/er counters, the whole procedure providing a useful check on 
preferred orientation in the powder sample and (after correction for 
detector efficiencies) a composite set of data from five counters. This 
will be shown in section A.5. 
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jre 2.6 The configuration of the Dido Curran neutron powder diff Tactometer, 
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(b) Dido 'Badger' pcv/zder di f f ractomete r 

Th)5 machine was a simpler version of the 'Curran' with a single 
detector and a monochromatic Incident wavelength of 1 . 1 1 3A reflected from 
the (Hi) lattice planes of a lead monochromator.' The neutron flux was 
less than half that of the 'Curran' and in the angular range of 9 = 5" to 15° 
the angular resolution was comparable. This machine was used for studying 
the upper phase transition in DI which involved scanning three powder 
peal<s at different temperatures. For this purpose the machine was quite 
adequate. 

(c) Pluto 'Panda' powder di ff ractometer 

In this machine the monochromator housing could be rotated and hence 
the wavelength and take-off angle could be varied, A 'squashed' 
germanium monochromator crystal was used in transmission at wavelengths of 

o o 

1.89A and 2J\A. The 'squashing' increased the mosaic spread of the 
crystal and hence the reflected intensity (Bacon 1955a). 

The energy distribution of the neutron beam emerging from the moderator 
in the reactor is approximately Maxwell ian with a characteristic tempera- 
ture of 60°C and a peak intensity at a wavelength of approximately K'^A. 
In contrast to the energy spectrum of an X-ray tube with its sharp 
characteristic peaks superimposed on a continuous background, the neutron 
spectrum has no characteristic wavelength. Moreover at neutron wavelengths 
greater than I.IA a monochromatic beam may become contaminated by shorter 
wavelengths. These are caused by second-and third-order Bragg reflections 
from the monochromator and occur at one half and one third of the 
fundamental or intended wavelength. 
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The configuration of the 'Panda' neutron powder di ff ractometer is 
shown in figures 2.8 and 2-9, The main beam from the core of the reactor 
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Figure 2»8 The configuration of the Pluto Panda neutron powder diffractometer. 
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enters Liirough a collimator at B and falls on the monochroniator at H. The 
monochromatic neutron beam is then Incident on the sample at S. A 
detector C moves on a circular track with its centre at S, The carriage A 
moves on another circular track centred at M enabling the take-off angle 
of the monochromator and the wavelength to be varied. The wavelength can 
be determined from the diffraction pattern of a cylinder of nickel powder, 
for which the lattice parameters are well known. The diffraction pattern 
of the nickel standard also reveals the presence of any short wavelength 
{higher order) contamination. 



To study diffraction patterns at long wavelengths It Is necessary to 
use the reflections from the (111) and (113) lattice planes of the 
germanium monochromator. Since germanium has the space group Fd3m (no. 227), 
the reflections (222) and (226) are absent preventing the occurrence of 
second-order contamination in the monochromatic Incident beam. Using the 
germanium (113) reflection and a wavelength of 1.89A the only contamination 
came from the third harmonic of the incident radiation, i.e. the (333) 

o 

germanium reflection for which A = 0.63A . This was negligible for three 
reasons : 

o 

(l) The contribution of the wavelength O.63A to the energy spectrum 
is smal 1 . 
(ii) The geometrical structure factor for the (333) germanium 
reflection Is smal 1 . 
(iii) The reflectivity of the sample increases with wavelength according 
to the formula. 

R « X3 (22) 

Thus the reflectivity for X = 0.63A is 27 times smaller than that of the 
fundi^menta 1 wavelength. 
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With a wavelength of 2.^A the third order contarninat ion from the 0.8A 
vvavelength was more serious due to Its greater contribution to the energy 
spectrum. This was eliminated by the use of a pyrolytlc graphite filter 
In which the hexagonal £-axis of the crystallites is oriented along the 
direction of the incide/it beam. The other crystal axes are oriented 
randomly in the (OOl) plane for each crystallite. Only neutrons with 
wavelengths which satisfy the Bragg equation for graphite are diffracted 
away from the direction of the incident beam Into Debye-Scherrer cones 
formed by the diffraction from crystallites having their a_-axes randomly 
oriented. These graphite powder diffraction neutron peal<s are absorbed In 
the cadmium shielding surrounding the filter. 



The neutron diffraction tn the filter reduces the intensity of the 
transmitted neutron beam at certain characteristic wavelengths. For 
pyrolytlc graphite aligned with Its C_-axis along the Incident monochromatic 
beam a low transmission occurs at 0.8A. The widths of the transmission 
minima caused by diffraction in the filter vary from one filter to the next 
according to the variation In the angular distribution of the direction of 
the £-axis in the filter. Accounts of these filters are given by Loopstra 
(1966) and Shiraneand Minkiewicz (1370). 

The advantages of the variable monochromator tal<e-off angle were; 



(i) By choosing suitably long wavelengths overlapping powder peaks 
could be partially or wholly resolved, 
(il) When the diffracted beam is parallel to the polychromatic primary 
neutron beam from the reactor, the di f f ractometer is In a 
parallel focusing geometry (Bacon 1955b). In this arrangement 
the angular resolution of the di ff ractometer is highest when 
the Bragg atigle of diffraction from the sample is equal to the 
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lasted at least ten days. As thermal hysteresis effects were expected to 
be present during the phase transition it was Important to keep the liquid 
helium bin full and the sample at a constant temperature. 

The sample tube used in this experiment was of a more complex design 
consisting of a narrow thin-walled vitreous silica tube approximately 10mm 
In diameter mounted co-axial ly inside a similar outer tube of nearly twice 
the diameter. The B.O.C. cryostat and the double-walled tube are shov^n in 
figure 2.10. The solid Dl sample S was in the inner tube. The annular 

exchange space between the tubes ES was kept at a pressure of less than 

-6 

10 mm of mercury by a vacuum pump at V. The whole assembly was immersed 

in the liquid helium bin LH surrounded by a radiation shield cooled by liquid 
nitrogen LN . 



I-';! 



The outer vacuum tail of the B.O.C. cryostat was made of aluminium 
with a thin window for the neutron beam B to pass through. Apart from a 
small PTFE spacing disc T of low thermal conductivity the sample in the 
inner tube was thermally isolated from its surroundings and so, in time it 
would have slowly cooled down to 4.2K by radiation to the liquid helium LH. 
By having two small heaters fitted above and below the sample S it was 
possible to vary its temperature above 4.2K. The minimum temperature 
attainable with this arrangement was limited (1) by the conduction of heat 
down the walls of the inner tube as well as along the heater and thermocouple 
leads H and TC and (ii) by thermal radiation from the top of the sample 
tube {vyhich was at room temperature). These effects were minimized by the 
use of very fine v/1 res and by platinum foil P in the sample tube to act 
as a radiation shield. A compromise had to be made however between the 
fineness of the electrical leads (to reduce thermal conduction) and the 
resultant increase of electrical resistance in the control circuits. The 
temperature measurement, thermocouple calibration and circuits used in 
temperature control will be described in appendix two. 
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Figure 2.10 The variable tomtx^rature ci-yostat used lor studylns the 
lower phase transition of Dl by neutron diffraction. 
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The sample v/as transferred from the storage tube to the specimen tube 
by vacuum sublimation. The specimen tube was cooled by liquid nitrogen 
in a small dewar to condense the Dl vapour coming from the storage tube 
which was kept at a slightly higher temperature. . During the sample trans- 
fer the annular exchange space ES surrounding the sample S was filled with 
helium gas at one atmosphere pressure. Convection of the helium gas in 
the exchange space enabled the latent heat of sublimation to be extracted 
from the Dl sample reasonably quickly and then transmitted to the liquid 
nitrogen. A sample transfer was completed in approximately h hours. 



I i 
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After the sample had been transferred to the double-walled specimen 
tube, the exchange space ES was evacuated and the specimen tube transferred 
to a helium cryostat. The sample then cooled down to ^OK by radiation to 
the liquid helium over a period of some hours. Continuous evacuation of 
the exchange space was necessary (i) to prevent condensation of oxygen or 
nitrogen impurities in the helium gas on the surface of the exchange space 
and (ii) to prevent the condensation of solid nitrogen from the atmosphere 
through the neoprene electrical seals onto the surface of the exchange 
space. Nitrogen has a high coherent neutron scattering cross-section and 
even small quantities can seriously change a powder diffraction pattern 
particularly at low 9-angles. 



The cryostat was refilled with liquid helium every 12 hours. The 
liquid helium transfer lines were precooled with liquid helium and the 
cryostat filled In approximately 5 mintues. The temperature of the outer 
vitreous silica tube and the sample varied negligibly during the transfer 
due to the Insulation of the evacuated exchange space. 



The bnck-l ransfer of the sample was a reversal of the transfer v;Ith 
the added problem of blockages in the narrow specimen tube. These were 
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always quickly cured by means of a hot-air blower and liquid nitrogen at 
suitable points. 
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CHAPTER THREE 



X-r ay powder diffraction study of the iodine sublattice in the 

three solid phases of PI' 



3. 1 Introductory theory of diffraction 

Thermal neutrons and X-rays when scattered by crystals give diffraction 
peaks (Bragg reflections) at angles which satisfy the Bragg equation: 



nA = 2dsine 



(23) 



where n is a positive integer, X is the wavelength of the incident radiation, 
d is the interplanar spacing and 6 is the diffraction angle (Bragg angle). 
The integrated intensity of a Bragg peak produced by a powder sample Is 
(Bacon I955ci) : 



'obs^^^^^ = 



C|F^(hkl)|A(e)je"^'^ 
s in6sin29 



(2^) 



where h, k and 1 are the Miller indices of the reflecting planes. C is a 
scale factor dependent on the flux of the incident .-KJ<t^''*i and the volume 
of the sample. F(hkl) is the structure factor and j is the multiplicity 

factor for the Bragg reflection. A(e) is the absorption correction for a 

-2W 
randomly polarized X-ray or neutron beam and e is the temperature 

correction factor. 

In equation (2A) the crystal structure is described by the structure 
factor F(hkl) which is defined by the expression: 



Na 



F(hkl) = I f^(6)exp [2TTi(hXj^ + ky^ + Iz^)] 
n=l 



(25) 



.A 
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where x , y and 2 are the fractional atomic co-ordinates of the n atom 
n n n 

in the unit cell. These are related to the absolute atomic co-ordinates 
x' , y' and z' by the expression: 



X = 



a 



y = 



xl 



and z = 



c 



(26) 



where a, b and c are the dimensions of the unit cell. ^p^^^ '^ ^^^ atomic 
scattering factor which will be discussed later. 

The temperature factor (Debye-Wal ler factor) W in equation (2A) 
describes the effects of thermal motion and can be represented in three 
forms : 

(a) An overall temperature factor B as in equation 2k: 



-2W r 2Bsin e, 
e = expL~ J 



(27) 



where 6 and X are the diffraction angle and wavelength defined in equation 
(23). 

(b) Each atom in the unit cell is assigned an isotropic vibration ampli- 
tude u which is related to the atomic temperature factor B^ by the 
express ion: 

(28) 



<u > 



B 
n 

16tt2 



The effect of the thermal motion of each atom on the observed intensity 
l(hkl) is now included in the expression for the structure factor given In 
equation (25) so that the observed intensity previously described by 
equation {2h) is given by: 



I . (hkl) = ^M^^j ^ |F'(hi<l)| 



obs 



s in0s in2e 



(29) 



, Na B sin e 

and F (hkl) - I f (6) exp [ZttI (hx^+ky^^+lz^) ] cxp [- — ] (30) 

n=l X^ 
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where the parameters are defined In equation (25). 

(c) In the general case a symmetric tensor of nine components {b?.} 

U 

is used to describe the thermal motion of each atom. The mean square 
amplitude of vibration of the n atom in the unit cell in the direction of 
the unit vector £ = (q,, <]2> ^7) 'S given by: 



3 3 . n 



i = l j=l ^ 2| *| I ''^1 

^v b X c 
where a_ = — — etc. are the reciprocal lattice vectors, 

a_. b^ X c 

A tensor of thermal motion {U..} can be defined as: 



and 



(U^.> 



%(a) 



{b?.} 

„LJL 



271^ I a' 



3 3 „ 
I I U?, q, q, 
i=l j=l 



I J ^1 ^J 



(31) 



(32) 



(33) 
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Since the tensor {b..} is symmetric, it has six independent terms. The 
number of independent terms may be further reduced by symmetry if the atoms 
occupy special positions in the unit cell, see Levy (1956). 

Using this general expression for the thermal vibration tensor the 
structure factor can be written as: 

Na 
F'(hl<l) = y f (e) exp[2Tri(hx + ky + Iz )] x 



n=l 



n ' n n 



exp[-(b"^h^ + b22k^ + b^^l^ + 2b'j2hk + 2b23kl + 2b'j'2h1)] (3A) 



m 



\i 



The preceding discussion has been pertinent to both X-ray and neutron 



82. 



scattering. There are however two points of difference. 

(a) The atomic scattering factor f (9). Since the nucleus is effectively 
a point source the angular distribution of the scattered neutrons is iso- 
tropic and the scattering factor f (e) is a constant, b . This is called 
the coherent scattering length and is a fundamental constant for each 
isotope of a particular element. X-rays are scattered by the electron 
core of an atom whose diameter is comparable to the wavelength of an X-ray. 
This X-ray scattering is described by an atomic scattering factor fpCs) which 
is proportional to the atomic number of the element and which decreases 
continuously with increasing scattering angle. 



(b) When X-rays are monochromat izcd by a crystal monochromator the observed 
intensity in equation {2k) is modified by a term p which corrects for 
polarization, Azaroff (1355)- 



2W 



(hkl) 



C r(hkl)[A(e)je" p 



sinesin2e 



(35) 



2 2 9 2 2 2 2 

where p = [(cos 2e^,cos p+sin p) cos 29+cos 2e^sin p+cos p] (36) 

(l+cos^2e^) 

n 



e^j is the Bragg angle of the monochromator and p is the angle between the 
normal to the plane of reflection of the monochromator and the normal to 
the plane of reflection of the sample. 



The following chapters will describe the determination of the crystal 
structure of the three solicj phases of Dl from the observed neutron and 
X-ray powder diffraction patterns. Standard algorithms have been derived 
and are widely used for refinement of crystal structures and lattice para- 
n^eters by iterative Iccts t-squares techniques. 
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The quality of the agreement between the intensities calculated for 
the proposed model and the observed intensities is assessed by a so-called 
R-factor which for the computer program used in this thesis (MPORFLS des- 
cribed in appendix one) is defined as: 




W d", -5 
n obs 



caic 



) 



(37) 



n=i 



I W I ^ 



n 2 



where fL. is the number of observations and I , and I , are the 
oos obs calc 

observed and calculated intensities respectively. S is the scale factor 
and tfie weight assigned to the n observation is W . 



3.2 X-ray diffraction experiments 

X-ray powder diffraction patterns were recorded in all three solid 
piiases of Dl. These served several purposes; to verify the purity of the 
samples in the initial period of the work, to measure the lattice parameters 
snd to provide values for the thermal parameters of the iodine atoms. 

(a) The low phase of Dj 

An X-ray powder diffraction pattern of the low phase of Dl was recorded 
at k.2K at Queen Mary College using the powder diff Tactometer and cryostat 
described in chapter two. It was very difficult to form good quality 
powder samples In the Lindemann glass capillaries used. Iodine has a 
iilgh X-ray absorption coefficient and for a cylindrical sample the "optimum" 
thickness producing maximum diffracted intensity is given by the expression: 
((nternational Tables for Crystallography 1959) 



nR = 1.5 



(38) 



where R is the radius of the sample and y Is the linear absorption coefficient 
defined by: 



I - \^ exp[-pt] 



(39) 
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where I is the transmitted intensity, 1 the incident Intensity and t the 
thickness of the sample. For Dl and CuK X-rays of wavelength l.S'tlZSA 
\i is 1015 cm and R (optimum) is 0.015 mm. The narrcwest available 
Lindemann glass capillary had a radius of COJiim, however capillaries of 
this size were always liable to become blocked while condensing the Dl gas 
into the sample space. To overcome this problem a 0.1mm radius capillary 
was used. The X-ray powder patterns obtained were of low intensity due to 

the absorption by the six layers of beryllium in the cryostat plus the fine 
detector colli mat ion necessary to eliminate the beryllium diffraction 
patterns. 

A suitable powder sample was grown and examined initially by 
Oebye-Scherrer photographs at 8i|.K and 50K shown in figure 3.1. (The first 
powder peak of beryllium occurs at a relatively large angle of 22° In 6). 
The first X-ray powder line recorded on the film occurred at about 9 = 10" 
for all of the samples Investigated. > 



An X-ray powder diffraction pattern of the low phase of Dl was recorded 
at 4,2K with relatively poor statistics using a count-rate meter and chart 
recorder. The chart recorder output is shown In figure 3.2. Because of 
the very low scattering power of X-rays for deuterium, the diffraction 
pattern obtained represents X-ray scattering from the Iodine atoms only. 



All the observed powder lines could be indexed on the basis of a face- 
centred orthorhombic lattice where the approximate lattice parameters were 
derived by using the program PARAM described in appendix one. The following 
values were obtained: 

a = 5-9^ ± .02 A 

b = 5.95 ± .02 A 

c = 6.22 + .02 A 

a = 6 *= y = 90' 



Pigur* 3.1 Low-mgle X-ray powd«r photogrmph of 
DI «t 84K uid 50K, ■ ^ 
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(b) The middle and high phases of PI 

X-ray powder diffraction patterns were recorded at temperatures of 
8^K, 93K, 103K. n7K, 128K, 129K, 130K and 133K in the middle and high 
phases using the Nonius goniometer described in chapter two. The patterns 
recorded indicated a face-centred orthorhombic unit cell for the iodine 
sublattice in the middle phase. The separation of the (020) . (200) doublet 
was only observed in a very few films, one of these is shown in figure 3-3. 
The iodine lattice in the high phase was found to be face-centred cubic. 
An X-ray photograph of the high phase of Dl tal<en at 130K is shown in 
figure 3.^- 

Good quality powder samples for the Nonius goniometer were not diffi- 
cult to prepare, approximately one in ten powder samples was satisfactory 
for diffraction studies. After some hours of X-ray exposure the sample 
decomposed giving rise to a series of low-angle diffraction spot due to 
free iodine. Ice crystals formed on the outside of the Lindemann glass 
capillary after approximately twelve hours. In one experiment a series of 
X-ray photographs were recorded on one film using a very narrow layer-line 
shield space. This technique is limited by the uniformity of the cassette 
circumference away from the central zone where the film is held in place 
by spr i ng cl i ps . 



The density profiles of the X-ray powder lines were recorded from these 
films using a Joyce-Loebl micro-dens it iometer. The values of the tempera- 
ture factor derived from different films varied considerably. This was 
probably due to the high correlation between the absorption correction and 
the temperature factor. An additional source of error was the uncertainty 
in the volume of the powder sample. Most samples were not solid cylinders 
but hollow cylinders of varying thicl^ness lining the inside wall of the 
capi 1 lary. 



Figure 3,4 X-ray powder photograph of DI at 130K (high phase). 
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[Figure 3,3 X-ray powder photograph o* DI at 117K (nlddle phase), 
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Lattice parameter refinements v;ere carried out on all recorded data 
and the results are shown in Table 3.1 and figure 3-5. The lattice para- 
meters of HI recorded by Mauer and Keffer (196A) are shown In figure 3.5 
as crosses In comparison with the lattice parameters of Dl (shown as dots) 
obtained by X-ray and neutron diffraction experiments. 

Some X-ray measurements on the middle and high phases of HI have been 
reported since the worl< of this thesis was carried out. Simon in 1970 
has measured the lattice parameters of HI at lOOK, I3OK and 210K using a 
photographic powder technique and has also found that the middle phase of HI 
is orthorhomblc. The complete results are given in table 3.1. 
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Flgur© 3,5 Lattice parameters of HI and DI as function of temperature. 
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C HAPTER FOUR 

Neutron powder diffraction study of the low phase of 

Dl at h.2K 



k. 1 Summary of the neutron diffraction experiments 

The crystal structure of the low phase of D! at 't.2K was first studied 
by neutron powder diffraction techniques using the Dido 'Curran' diffracto- 
meter (see §2. 3a) and a neutron wavelength of 1.06a. The sample was con- 
tained in a vitreous silica tube immersed in liquid helium Inside a vanadium- 
tall cryostat (see 52. A). The neutron powder diffraction pattern recorded 
in this experiment is shown in figure ^.1; the broad low hump in the 
background at low 29-angles is due to neutrons scattered from the liquid 
helium and the walls of the vitreous silica sample container. As mentioned 
earlier In §2.^, the neutron scattering of vanadium is almost totally 
Incoherent and isotropic. This experiment was subsequently repeated with 
the same di f f ractometer and the same neutron wavelength but using a Dl 
sample prepared by a different method (see 5§2.1a, b) . The neutron powder 
pattern recorded in this experiment was similar to the one shown. 



Neutron powder diffraction patterns of the low phase of Dl were also 
recorded with wavelengths of 2.66A and 1.886a using the Pluto 'Panda' 
dl ff ractometer under similar experimental conditions as described above. 
The former of these (X = 2. 668) was a trial experiment and its results were 
not used in the final analysis due to the suspected presence of X/3 contamina- 
tion in the monochromatic neutron beam. Nevertheless, the improved 
angular resolution of the diffraction pattern recorded with X = 2.66A showed 
that (Dony of iho. powder peaks, which appeared to be single in the diffraction 
pattern recorded with a wavelength of 1.06a using the 'Curran' di f f ractonieter, 
were In fact multiple peaks containing several components. 
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Finally, the neutron powder diffraction pattern of the lov; phase of Dl 
was recorded with the 'Panda' di f f ractometer using a wavelength of 1.886A. 
This diffraction pattern had good angular resolution over a wide angular 
range and also good statistical accuracy; it is shown in figure ^1.2. 



i),2 The u nit cell of the low phase 

The most interesting feature of the neutron powder diffract.ion of the 
low phase of D! was the presence of a distinct peal< at low angles not seen 
in the X-ray powder pattern (see the peak (llO) in figure ^.1). Using 
approximate lattice parameters derived from the X-ray powder diffraction . 
patterns of the low and middle phases of HI and Di (see §3-2) all neutron 
powder peaks observed in the low phase of Dl could be assigned Miller 
indices. A large orthorhombic unit cell was assumed with all three edges 
twice as large as the edges of the orthorhombic unit cell derived from the 
X-ray powder diffraction patterns of the low phase. This large orthorhombic 
unit cell contained 32 Dl molecules. 

Since it was t<nown from X-ray powder diffraction studies that the 
iodine atoms formed a face-centred orthorhombic lattice with four iodine 
atoms in the unit cell, the much larger unit cell revealed by the neutron 
powder diffraction had to be attributed to a more complex arrangement of 
the deuterium atoms. This could not be detected by X-ray diffraction 
techniques because the X-ray scattering of the deuterium atoms is negligible 
compared with the X-ray scattering of the Iodine atoms. Thus, X-ray 
diffraction of Dl reveals only the position of the iodine atoms. In con- 
trast, the neutron scattering cross-sections of deuterium and iodine are 
approximately equal, hence the neutron diffraction pattern reveals the location 
of the deuterium atoms as well. 
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A least-squares refinement of the lattice parameters was carried out 
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using the computer program PARAM {Appendix one). Since crystal- 
Igionochromated neutron beams give sharp peaks only in the so-called 
'parallel-focusing' geometry (see Bacon 1955b), the neutron powder diffrac- 
tion pattern was recorded only on one side of the incident beam. As a 
result the direction (zero-angle) of the incident beam had to be determined 
by calculation. It was varied via the computer program PAI^M used to 
refine the lattice parameters. Table k,\ lists the results of the lattice 
parameter refinement obtained with the program PARAM for the low phase of 
Dl at A.2K. The refinement Is based on the powder diffraction pattern 
recorded with the 'Panda' di ff Tactometer using a neutron wavelength of 



^^1 



\ = 1.886S. 



At this stage it will be useful to summarize briefly the structure of 
the middle phase of Dl which will be described in greater detail in Chapter 
five. The crystal structure of the middle phase has a face-centred 
orthorhombic unit cell containing four molecules. The iodine atoms occupy 
face-centred positions and each deuterium atom has two equilibrium positions 
symmetrical about the Iodine atoms. The two equilibrium positions are In the 
(001) plane 180° apart and the deuterium atoms 'flip' between them at random 
(two-fold disorder). The space group Is Bbcm (no. 64) and a projection of 
the time-average of the structure onto the (OOl) plane is shown in 
figure 5-4. 
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The calorimetric data described In §1.2 show a X-anomaly In the 
specific heat of Dl in the lower phase transition region associated with an 
anomalous increase in entropy. This 'excess' entropy can be interpreted 
as a change in the conf I gurat lonal entropy of the solid due to an increase 

In disorder on transition and will be discussed In detail In chapter seven. 
Here it is sufficient to state that the observed anomalous entropy cliange 
between the low and middle phases of Dl is; 
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Observed 29 

15.60 

30.50 

30.60 

31.18 

3it.12 

36-88 

37.^0 

39.0it 

iiO.itO 

it 1.80 

i*5.15 

i*7.55 

ii8.52 

itS.Ol 

51.30 

52.80 

53.56 

55.00 

56.00 

59.23 

61.80 

63.01 

63.58 

65.0it 

68.90 

69.58 

71.78 

72.68 

7it.it0 

75.ii8 

78.00 

80.10 

8it.20 

87.60 

95. Ai* 



Calculated 29 

15.'*3l88 
30.3550i4 
30.it922ii 
3l.l9AOit 
3i|. 11903 
37.0it98l 
37.2798it 
39.02716 
itO.'t6930 
i)1.822it6 
A5. 25101 
A7. 56985 
it8. 62210 

itS. 89929 

51.2i45it0 

52.79119 

53.57296 

55.03itn 

56.02260 

59.258it6 

6I.7836I 

63.15091 

63.'t6l99 

65.05976 

68.9151'! 

69.it567it 

71.85058 

72.6it738 

7^4.20910 

75.52ii23 

78. 08179 

80.21388 

8it.llit57 
87.53A6i+ 

95.'*3A35 



20 



lbs- 



■29 



ale 



0.1it8l2 
0.1it'*96 
0.10776 
■O.Oli+Qij 
0.00097 
■0.16981 
0.12016 
0.01 28it 
■0.06930 
•0.02246 
10101 

01985 
10210 

11071 
05it60 
00881 
01296 
03i+11 
02260 
028it6 
0.01639 
o.litogi 
. 1 1 80 1 

01976 

0151'+ 
12326 
07058 
03262 
19090 
Oitit23 
08179 
11388 
0.085^43 
0.06536 
0.00565 



-0 

-0 

-0 







-0 

-0 

-0 

-0 



-0 

-0 



-0 





-0 
-0 
-0 



Old 

11.801530 

11.872205 

12.857706 

90.00 

90.00 

90.00 



New 

11.80 1530 

11.872205 

12.857706 

90.00 

90.00 

90.00 



The error of fit is 2.015030819 



Delta 

0.00 
-0.00 
0.00 
0.00 
0.00 
0.00 



Estimated 
Std. Dev. 

0.007836 
0.007^^09 

0. 00^1890 

0.00 
0.00 
0.00 



Correlation Coefficients 



■0.696itit 



■0.10287 
■0.27791 



a 
0. 
0. 
0. 



e 
0. 
0. 
0. 

0. 



Y 
0. 
0. 
0. 

0, 



TABLE it.l. Results of the lattice parameter refinement obtained with the 
program PARAM for the low phase of Dl at ii.2K. The refinement is based j.pcr 
the powder diffraction pattern recorded with the 'Panda' with a neutron wav.- 
lenath of X = 1.866^. 
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If this change is due to an Increase in conf i gurat ional entropy it 
should be equal to: 



AS,_2 = 1.101 = k In^ 



i^O) 



where k is the Boltzmann constant and njn^ are the number of equilibrium 
positions in the low and middle phases respectively. For the lower phase 
transition in Dl ; 



mi ddle 
low 



= 1.7^ 



(^1) 



This value supports the assumption that the low phase Is ordered and 
the middle phase is two-fold disordered. The heavy iodine atoms are 
expected to occupy essentially the same ordered positions in both the low and 
middle phases, while the deuterium atoms are expected to occupy positions 
which are likely to coincide with one of their equilibrium positions In the 
middle phase. 



H 



On passing through the lower phase transition point the two-fold 
disordered positions of the middle phase are expected to be achieved by 
re-orientation of the Dl molecules through l80°. Though re-or ientatlons 
through 90" seem to be less probable, they cannot be excluded at present 
from consideration. 



Having found the true unit cell of the low phase the next step was to 
determine the space group from the systematic absences In the observed 
powder pattern and then propose suitable models for the structure of the 
low phase. 
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4.3 Models base d on orthorhonbic space groups 

The neutron powder diffraction pattern of the low phase was initially 
assigned Miller indices on the basts of a large orthorhombic unit cell con- 
taining 32 Dl molecules. The observed peal<5 obeyed the selection rules 
for a face-centred lattice; 



h + k, k + 1, (h + 1) = 2n 



(^2) 



There are five face-centred orthorhombic space groups; F222 (No. 22), 
Fmm2 (No. '♦a), Fdd2 (No. A3), Fmmm (No. 69) and Fddd (No. 70). 
Of these the space group Fddd can be excluded because it also requires that 



in 



hkO; h + k = An (h,k = 2n) 



(A3) 



hence it would not allow the observed (A20) powder peak. Furthermore, the 
face-centred orthorhombic unit cell containing 32 Dl molecules should be 
compatible with the iodine sublattlce derived from the X-ray powder pattern 
of the low phase. This requires that the large orthorhombic unit cell 
should accommodate eight small face-centred orthorhombic unit cells 
containing four iodine atoms each, forming a 'super-face-centred' iodine 
sublattlce. In addition, the symmetry operations of the face-centred 
orthorhombic space group (when applied to the deuterium atoms) should provide 
an ordered subgroup of the two-fold disordered D! middle phase described In 
the previous section. 

It was found that, out of the four remaining orthorhombic space groups, 
F222 (No. 22) and Fdd2 (No. A3) could produce the 'super-face-centred' iodine 
"subiattice. The 32 iodine positions in space group F222 (No. 22) can 
occupy four 8-fold special positions, described in Wyckoff notation, as; 



1 



Se (x,0,0) , (x,0,0) 



(AA) 



100. 



8f (O.y.O), (0,7.0) 
8g (0,0,z), {0,0»D 
8i {l,yM, ilA-yA) 



('15) 
(A6) 

iky) 



Choosing x, y and z separately as J leads to an iodine sublattice of 
the required type. However, it is not possible to choose co-ordinates for 
the deuterium atoms which provide an ordered subgroup of the middle phase. 



In space group Fdd2 (No. 43) the 32 Iodine atoms can occupy two 16-fold 
special positions, described In Wyckoff notation as; 



I6b (x,y,z), (x,y,z), (^x, ^+y. l+z) , (i+x, ^y, l+z) 
(0,0,0; 0,l,li 1,0,1; hho) + 



(48) 



Choosing (x,y,z) as (1,1.0) and Ci,|,0) respectively leads to an Iodine 
sublattice of the required l<ind and a deuterium lattice which Is a subgroup 
of the mi ddle phase. 



Structure models of the low phase based on the orthorhombic space 
group Fdd2 (No. 43) were generated in the following way. The deuterium 
atoms were constrained to remain in the (OOl) plane and occupy one of four 
positions. Two of these coincided with positions they occupy in the middle 
phase and the others were derived from these by a 90° rotation of the 
molecule about the orthorhombic c-axis passing through the iodine atom. 
Since there were two molecules In the asymmetric unit of this structure 
there were 2^ = 16 possible rrradels (of which many were equivalent). There 
were only three different models among the sixteen possible and these are 
shown in figures 4.3 to. 4.5- 

The first model could be easily discarded since it could be described 
by a snail unit cell containing four molecules. The second model was 
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2=1/2 





Figure 4,3 Proposed structure for DI low phase space-group Kdd2 (No.43), 

model No* 1. 
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Figure 4,4 Proi)Osed structure for DI low phaee space-group Fdd2 (No. 43), 

model No. 2. 
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Fif^ure 4,5 Proposed structure for DI low phaoe space-group Fdd2 (No. 43), 

r idol Ko. 3. 
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theoretically possible but physically unlikely because the electrostatic 
repulsion is expected to prevent the deuterium atoms of neighbouring mole- 
cules facing one another. The third model was slightly more realistic 
since it contained ten adjacent deuterium atoms. Eventually all three 
models had to be discarded because the diffraction patterns calculated for 
them were quite unlike the observed pattern. 

k.k Mode ls based o n mono cl in jc space groups 

As an alternative to the large face-centred orthorhombic unit cell 
containing 32 Dl molecules the observed neutron powder peaks of the low 
phase could also be indexed on the basis of a C-centred monoclinic cell 
containing 16 molecules. The edges of the face-centred large orthorhombic 
unit cell can be transformed to the edges of the smaller monoclinic cell in 
six different ways. These transformations are described by the following 
matrices and are illustrated in figure ^.6. 



No.l 



fo 








a_' 




1 





b_' 
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1 





c' 




_"* 


i 




monoc 


1 in i c 





a_ 
b_ 
c_ 
orthorhombi c 



a, = 



b, = 



c, = 



a 
b 
/a +c 



11.802^ 

11.872S (kS) 

8.7278 



No. 2 



a_' 
b' 
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= 
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-i 
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nxjnoc 


1 inic 







b_ 
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orthorhombi c 



a^ = 



b., = 



■2~-2 



= /a +c = 



12.858S 

11.8728 (50) 
8.7278 



No. 3 





1 
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monocTi n i c 
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b_ 

£ 

orthorhombic 



b., = a 



= /b*' 



•2—2 



+c = 



12.858^ 
11.8028 (£0 
8.7508 



I 'i: 
ill 



4 



® 



® 



(D 
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c-a 



c.= 



bi = b 



A C,= - 



a + c 



c 



vQi = a 



C 



92= c 



® 



a3=c 


■^ h r* 




-^ U — ^^ 



'^b,= -g 



c 



vb3 = g 



c-b 
^ 2 



c 



^5=-l=2 



a+b\^'""'""-^ ^5=0 



c 



g5=g M' 




C 



6- 2 



rlBuro 4.6 Tl.e Blx r.ietl>cdE of transforming tho large orti,orliO.-nblc 
unit ci-,11 to the KSallcr mcmncllnic unit coli. 



b, = b 
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g6=b 



lOA. 
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a_ 
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= b 



~ a 



= 11.8728 



2 2 



« /b +c = 



11.8028 (52) 
8.7508 



b. = 



2 






■f 


a 
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n 


.8028 ' 


c 
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12 


.8588 (53) 


/^^+b^ 
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.3708 1 



No. 6 



a' 
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b' 
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i 
1 inic 


-i 






c 
L— J 

orthorhombic 



^6 



= b 



= v^2 



a +b = 



11.8728 
12.8588 (5M 
8.3708 



These matrices transform both the unit cell edges and the Miller indices 
of lattice planes. Table k.2 lists the Miller indices of the Bragg peaks 
observed in the negtron powder diffraction pattern'of the low phase of Dl 
In the orthorhombic and six monoclinic unit cells. Since the lattice 
parameters of all these monoclinic unit cells are simply related to those 
of the large orthorhombic cell, it was unnecessary to carry out least squares 
refinement of the lattice parameters in the monoclinic system. 

After the orthorhombic Miller indices had been transformed Into the six 
corresponding sets of monoclinic indices by means of the matrices (49) to 
(5^) they were examined for systematic absences. There are five C-centred 
monoclinic space groups; 

C2 (No. 5), Cm (No. 8), Cc (No.9) . C2/m (No. 12) and C2/c (No. 15). 



Of these C2, Cm and C2/m have no selection rules other than those 
characterizing a C-centred lattice, i.e. 

in hk\\ h+k = 2n, hOl; (h=2n) , OkO; (k=2n) (55) 
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TABLE 4.2. Miller indices of Bragg peaks observed in the neutron powder 
diffraction pattern of Dl in the orthorhombic and the six 
monoclinic unit cells. 
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On the other hand, the space groups Cc and C2/c have the additional 
requi rements that; 



in hOl ; 1 = 2n , (h«2n) 



(56) 



The types of orthorhombi c Miller Indices necessary to give a monocl I n i c 
Bragg peak of type (hOl)^, where 1 = 2n+l, are; 



Unit cells Nos. 1 and 2; (hOO^: 1^-h^ = 4n-2 
e.g. (002)^, (204)^. ... 



(57) 



Unit cells Nos. 3 and A; (01^1)^; l^-l<^ = ^n-2 
e.g. (002)^, (20A) ... 



(58) 



Unit cells Nos. 5 and 6; (hkO)^; h^-k^ = ^n-2 
e.g. (^20)^ which was observed. 



(59) 



Of the orthorhombic peaks listed In table 4.1 only the (420)^ peak was 
observed in the low phase of Dl , which excluded space groups Cc and C2/c 
for the monocl inic unit cells nos. 5 and 6. 



As discussed in 5^.2 the mfddfe phase of Of was expected to be 
related to the low phase by a 180° rotation of the Dl molecules In the (OOl) 
plane about an axis passing through the iodine atom. In addition the mono- 
cffnic space group had to accommdate the face-centred iodine sub~iattice 
with four atoms in the unit cell. In the space groups mentioned above the 
only symmetry operations besides the C-centring are £-glides, mirror planes 
and two-fold diad rotations about the b^-axis. Figures k.l to A. 12 show 
the middle phase of Dl referred to the six C-centred monocl inic unit cells 
to see whether these arrangements can be described by the above-mentioned 
symmetry operations. 
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Igures 4.10 - 4.12 Middle phase structure of DI drawn In a monoclinlc 

unit cell type No.s 4-6, 



no. 



Figures ^.7 to ^.12 show cltnographic projections of the middle phase 
arrangement onto the orthorhombic (OOl) plane. In cells No. 1 and 2 no 
subgroups of the atomic arrangement can have a £-glide or a mirror plane 
perpendicular to the b_-axis or a two-fold diad rotation parallel to the 
b-axis. The arrangements are pseudo-monocl inic and should be described 
by primitive trtclinlc cells consisting of eight molecules each. In cells 
No. 3 and No-'+ subgroups can have £- and n-glldes perpendicular to and 
2 and 2, rotations parallel to the b;-axis. Figures ^.13 and k.]k show 
clinographic projections of cells No. 3 and No.^ drawn onto the (OOl) mono- 
clinic plane to show the action of the monoclinic symmetry operations of 
space group Cc. It can be seen that in cells No. 5 and No. 6 £-g]ides 
and diad rotations are excluded. However, mirror planes perpendicular to 
the b-axis are allowed, thus these unit cells can accommodate atomic 
arrangements which are subgroups of the middle phase and have monoclinic space 
group symmetry Cm. 

Thus the following monoclinic space groups are compatible with the six 
possi ble un i t eel 1 s ; 



Cells No.l and No. 2 
Cells No. 3 and No.^t 
Cel Is No. 5 and No. 6 



None 

Cc and C2/c 

Cm 



If in the monoclinic unit cells No. 3 and No.^ each deuterium atom Is 
allowed to occupy one of the two disordered positions, it is possible to 
obtain eight different atomic arrangements for each unit cell, a total of 
sixteen. There Is a basis of two molecules In space group C2/c and four 
molecules (labelled 1 to A in figures ^.13 and k.]k) in the lower symmetry 
space group Cc, so allowing 2 (16) arrangements per cell. This Is reduced 
to eight per cell as the choice of the original direction is arbitrary. 
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igure 4,13 Middle phase structure of DI drawn In a monoclinic unit 

cell, type No. 3, projected onto the (OOI) monoclinic plane. 



igure 4.14 Middle phas» structure of DI drawn in a monoclinic unit 

cell, type No, 4, projected onto the (ooi) monoclinic plane. 
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The eight models generated with cell No. 3 are identical to the eight 
models of cell No.'*. Of these eight models three pairs are equivalent to 
each other so that there are five essentially different models which are shown 
in cell No.^ in figures 4.15 to A. 19. Model No. 1 corresponds to the 
ordered low phase structures of DCl, DBr and DF. The first two compounds 
are described by space group Bb2.m (No. 36) and the last by Bm2.b (No. 36) 
due to a change in the orientation of the major axis. ^k)del No. 2 is very 
similar to model No.l except that In model No. 2 the zig-zag chains are anti- 
parallel while in model No.l they are parallel. However, both models No.l 
and No. 2 can be referred to a small orthorhomblc unit cell containing four 
nx)lecul es . 

If each molecule In the low phase is allowed to re-orientate through 90° 
as compared with its position In the two-fold disordered middle phase, then 
this would lead to 128 further models (many of which would be equivalent 
to each other). Moreover several of these models are expected to be non- 
physical (as In the orthorhomblc Fdd2 models) because of unlikely configura- 
tions of adjacent dipoles. A computer program LPSG (Low Phase Symmetry 
Generator) described In Appendix one, was written to calculate the 
Intensities of the first thirteen powder diffraction peaks for each of these 
128 models. The 128 sets of calculated intensities were compared with the 
observed intensity pattern of the low phase of Dl , particularly in the 
region between the first low phase peak (lIO) and the subsequent Intense 
multiple peak (130-220) which remains practically unchanged In the middle 
phase diffraction pattern (figures A. 2 and 5.1). In this region there are 
no observed peaks although they are not excluded by any space group symmetry. 
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Eight of these models were found to give rise to calculated intensity 
distributions similar to the observed powder pattern. These were two sets 
of four equivalent structures of which one was model No. 5 (figure ^.19) and 
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4,17 Proposed structure for DI low phase space-group 
Co (Ko.9) models No. 1-3, 
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the other which is of a lower symmetry is shown in figure ^».20 as model 
No. 6. 



The diffraction patterns calculated for the four models with large unit 
cells, i.e. Nos. 3, ^, 5 and 6 shown in figures 4.17 to A, 20, are shown in 
figures k.2\ to A. 24 as continuous lines. The crosses represent the powder 
diffraction pattern of Dl at 4.2K recorded by the 'Curran' di ff Tactometer 
using a wavelength of X = 1.06A, Models No. 3 and No. 4 were very similar 
to each other (see figures 4.17 and 4.18) and the calculated diffraction 
pattern based upon model Mo. 5 was found to give the most reasonable fit, 
particularly at low angles (see figure 4.23). Models No. 3, 4 and 6 could 
be described by space group Cc, while model No.5 had an additional centre of 
symmetry characteristic of space group C2/c. 

Although model No.5 based on space group C2/c gave the best agreement 
with the observed data, other possible monoclinic space groups were also 
investigated for completeness. In order to examine all of these possible 
models the large face-centred orthorhombic unit cell was again considered. 

A basis of eight molecules was chosen (shown as 1-8 in figure 4.25) and 
the deuterium atoms were allowed to occupy one of the two disordered 
positions of the middle phase. The only permitted symmetry operation was 
that of face-centring so that the total number of models generated 
was 2' (128). Of these 128 many were equivalent to one another so that 
there were 17 different structures. 5 of these had small unit cells and 
of the remaining 12, 3 had been previously considered in the cells Nos. 3 and 
4 and are shown In figures 4.17 to 4.19. 

Diffraction patterns were calculated for the remaining 9 structure but 
they were all found to disagree with the observed neutron powder diffraction 
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Figures 4.21 and 4«22 Calculated (continuous line) and observed (crosses) 

intensities for DI low phase, Curran 4,2K, 
models No. 3 and 4. 
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Figure 4.25 
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Figure 4,25 Middle phase structure of DI drawn in an orthorhombic unit cell 

showing the 8 molecules of the basis chosen to produce 128 models. 



Figure 4,26 Middle phase structure of DI drawn In a monocllnlc unit cell 
type No, 6 showing the 6 molecules of the basis chosen to 
produce 32 models with space-group Cm (No.S). 
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pa 



ttern of the low phase of Dl , particularly at low angles 



As discussed earlier cells No. 5 and No. 6 could contain structures with 
space group symmetry Cm, If the structure is assumed to form a subgroup 
of the middle phase (as shown in figures A. 11 and ^.12) it is necessary to 
consider 6 iodine atoms in cell No. 6. Of these atoms, h are in the two- 
fold Wyckoff special position '2a' and 2 are in the general position 'Ab', 



I .e. 



2a (x.O.z) 
Ab (x,y,z), (x,7,z) 
(0,0,0; 1,1,0) + 



(60) 
(61) 



Choosing (x,y,z) as (0,0,0), (0,0,^), (^,0,0), (1,0,1) and as (1,1,0), 
ililA) respectively leads to a structure which Is a subgroup of the middle 
phase. The unit cell No. 6 is shown In figure 4.26 with the six atoms 
labelled as 1-6. Allowing the deuterium atoms to occupy one of the two 
equilibrium positions of the middle phase leads to 2 (32) models for both 
cells No. 5 and No. 6. Examination of these models shows that they are all 
included in the 17 different structures generated within the large face- 
centred orthorhombic unit cell. Only the most general case remains of 8 
molecules in the orthorhombic unit cellwith 90° re-orientat Ion. This leads 
to 2x4'(32768) models, an excessively large number to examine. Table A. 3 
summarizes the unit cells and space groups considered. 
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k.S Refinement of the models based on space group C2/c (No. 15) 

Two sets of data were used in the refinement; a set of data recorded 
with X = 1.8868 using the 'Panda' dif fractometer and a combined set of data 
derived from the ten powder patterns recorded with the five counters of the 
'Curran' di f fractometer using X = 1.068. The ten 'Curran' diffraction 
patterns were corrected separately for background and detector efficiency 



1 '-'fii 

is 
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before combining them into a single composite data set. The ten corrected 
diffraction patterns showed only statistical variation between each other 
indicating a lacl< of preferred orientation. 



The Rietveld program described in Appendix one was used in an attempt 
to refine the model No. 5 shown in figures 4.19 and ij.Z?. It was necessary 
to use this program for two reasons. 

a) The superposition of monoclinic powder peaks 

Because the lattice of the low phase is pseudo-orthorhombic , the 

orthorhombic powder peak (^q^^qIq) which has a multiplicity of eight, 

becomes two superimposed monoclinic peaks, 



Cell No. 3 [lo.l<o.(-V^)l and [T^,k^^i-2^)] 



and Cell No. 4 [k^,h^.- (-2^)] and [ k^ , h^ , - (-2^) ] 



(62) 



(63) 



:■■::■: .t\ 



each with a multiplicity of four. 



On the other hand in the monoclinic space group C2/c, the geometrical 
structure factors obey the following identities; 



1-/1 



F(hkl) = F(hkD 
'm ^ 'm 



Thus for the two unit cells, 

Cell No. 3 l[l^,k^,(-^)] ^ |[T ,k^,-(_2_°)] 

and 

Cell No. 4 |[k^,h^.-(-2_E^)] ^ l[k^,h^,-(-2_°)] 



(6M 



F(hkl)^ = F(hkl)^ ^ F(hkl)^ = F(hkT)^ (W2". 1^=2") (65) 

F(hkl)^ =-F(hkl)^ ¥ F(hkl)^ =-F{hk-n^ W2"' ln,=2nH-l) (66) 



(67) 



(68) 
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Figure 4,27 Proposed structure for DI low phase space group C2/c 

model No. 5 In unit cell No, 4. 
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although for 
Cell No. 3 



'^Vo' (^^^ = '^W' -(^^1 



^'^^ k +1 k -1 

Cell Uo.k e[k^.h^.-(-^)] = e[k^,h^, -(-2^)] 



(69) 



(70) 



Thus in the powder diffraction pattern of the low phase, each peak 
is the superposition of two whose intensities are generally different. 
A similar situation occurs in a cubic powder pattern where the peaks (115), 
(333) and (006) , (22^) always overlap. 

b) The large number of unresolved peaks 

The powder intensity distributions of the low phase of Dl recorded by 
the 'Curran' and 'Panda' di f f ractometers at ^.2K, after correction for 
background are shown in figures A. 28 and k.lS (a,b) by crosses; the snraoth 
line is the intensity distribution calculated for the proposed model No. 5 
and the vertical lines represent the angular positions, of the powder peaks 
compatible with the space group. No powder data were recorded in the 
angular range of 17-28° with the 'Panda' dl f f ractometer (see figure A. 29a). 
The total number of powder peaks allowed in the angular range covered by the 
'Panda' experiment Is 250 (effectively 125 considering the effect of super- 
position). Very few of these v/ere recorded as 'single' peaks despite the 
relatively good angular resolution of the Instrument. In addition the high 
angle peaks of the 'Curran' diffraction pattern (see figure A. 28b) are broad. 
The version of the Rietveld program used for this work (Hewat 1970 was 
unable to cope with such broad peaks. 

The intensity distributions shown in figures ^.28 and ^1.29 are cal- 
culated for model No. 5 shown in figures A. 19 and ^.27. The intramolecular 
bond length chosen initially for the model was that of the middle phase and 
the molecules were constrained to point towards nearest neighbour Iodine 
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bonds. Apart from disagreements at higher angles, where the poorly defined 
background was expected to cause difficulties, the main differences between 
the observed and calculated diffraction patterns shown in figures '».28a and 
4.29a are in the peaks labelled 1 to h. These peaks have the following 

Mi 1 ler indi ces ; 



Peak No. 
(1) 
(2) 
(3) 



Cell No. 3 

(no)+(Tii) 

(3lT)+(312) 
(131)+(T32) 
(330)+(l33) 



Cell Uo.h 

(iTo)+(iTn 
(iTi)+(ilD 

(331) + (332) 
(3T0) + (33D 



■V 



For peaks (l) and (2) the calculated intensity is greater than the 
observed one and for peaks (3) and {h) the situation Is reversed; in unit 
cell No. 3 two peaks are of type (hhl) while in unit cell No, 4 all four peaks 
are of type (hhl). The peaks (5) and (6) in figures 4.28 and 4.29 are 
the second orders of peaks (1) and (2). The good agreement between the 
observed and calculated values of peaks (5) and (6) as well as the close 
similarity of the ten 'Curran' diffraction patterns, practically excluded 
preferred orientation effects as a possible cause for the discrepancies in 
peaks (1) to (4). 



To account for the observed discrepancies the deuterium atoms were 
allowed to rotate about the iodine atoms in the orthorhombic (OOl) plane. 
Appropriate pseudo-symmetry operations were used to connect the two mole- 
cules of the asymmetric unit (see molecules (l) and (2) in figure 4.27). 
The angular rotation covered a range of 0.5° to 7-5° in both directions and 
the intramolecular bond length was also allowed to vary separately from 
I.2R to 1.68. These distortions of the model produced either only slight 
improvements or, in some cases, gross disagreement between the observed and 
calculated diffraction patterns. 
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An examination of the geometrical structure for these peaks was more 
informative. The geometrical structure factors for the space group C2/c are 



A = 8cos2Tr(hx+lz)cos2Trky 
B « 



h+k=2n, l=2m 



(71) 



A = -8sin2TT(hx+lz)sin2TTky 
B = 



h+k=2n, l=2m+l 



(72) 



and A = B = 



h+k=2n+l 



(73) 



For the peaks labelled (1) and (2) in figures ^.28 and ^.29 based on 
model No. 5 (see figure A. 27) the geometrical structure factors are: 



Cell No. 3 

F^lT0) + F^Tll)=F^(l) 



F^(3lT) + F^Il2) = F^(2) 



2 

= 6^ { [cos2fx,. cos2iTy. .+cos2 TT x,2COs2TTy-„] + (7^) 

[sin2Tr(z,j-x,psIn2Try-^+sin27r(z-2-X32^^'"^^y'?2^ ^ 

= 6^ {[co52TT(2z-^-3x-^)cos2TTy j+cos2TT(2z.2-3x,2)cos2Try „) 
+ [sin2Tr(3x-^-z j)sin2Try j+sin2TT(3x.2-Z-2)sin2TTy ] } 

(75) 



Cell No.J^ 



F^(no)+F^(lTD = F^(l) = 6^ {[cos2Trx^,cos2iTy^^+cos2TTX2j2<=os2TTy^2^^ + ^76) 

[sin2Tr(xj^^-z^j)sm27ry^^+sin2TT(x2^2"2/,2^^'"^^^42^ ^ 



M 



r .~;i 



F (ni) + F^(lT2) = F^(2) = 6^ {[cos27r(x^j-2z^Pcos27ry^j+cos2T7(x^2'2z2,2)^°s2Try^2] 

+ [sln2Tr(xj^^+Z2^|)sIn2Try2j^+sin2Tr(X|^2'*"^Z,2^^'"^^"'^A2^ -^ 

(77) 

where the co-ordinates x. . etc. describe the j atom of the asymmetric unit 
{j = l ,2) and i is the unit cell chosen (l=3,^). 
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The unit cell edge-vectors (a b c ) and the Miller indices (h k 1 ) of 

— o—o— o o o o 

the large orthorhombic unit cell can be transformed to those of the mono- 
clinic cells No. 3 and No. A by the matrices (51) and (52). The fractional 

atomic co-ordinates (x ,y ,2 ) are transformed by. a different matrix (the 

000 

transpose of the inverse unit-cell transformation matrix). If 

— ~ i^' ■ ; '^ ^^^ '^'^' ^ ^^' ' matrix and ^ = 7 Rr,i ? 's the fractional 

co-ordinate matrix then R, , = (A )... 
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(51) 

(78) 



(79) 
(80) 



(52) 
(81) 



(82) 
(83) 



If the expressions for the geometrical structure factors of the peaks 
(1) and (2) given by equations (7^) - (77) are transformed by the matrices 
(79) - (83) into the orthorhombic unit cell it is seen that the geometrical 
structure factors of these peaks are tlie same in both cells No. 3 and No. 4 
and are given by: 



h (1)=6A {[cos2Tr(yQj+ZQ^)cos27rXQ^+cos2TT(yQ2+ZQ2^^<=s2TrXQ2^ "*■ 
[sin2TT(y -2 ,)sin2w +sin2iT(y---z._)sin2TiX-_] } 



01 01 



02 02' 



m) 



-•«^ 



I vfl 
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F^(2)=6i+ {[cos2Tr(y-,-32-,)cos2TTx ,+cos2TT(y -3z )cos2Trx ] + (85) 



'01 - 01 



^01 



02 ■'02' 



[sln2T:(y-,+3z„,)sin2TTX ,+sin27T(y +3z )sin2TTX ] } 



'01 ■" 01 



01 



02- 



Hence if the Dl molecules lie in the (001)^ plane It is found that for 

z,. 
Cell No.3 x^ = -|^ and Cell No. A Z/,, = 

and for the orthorhombic cell, z-. = 0. 

Having shown that the geometrical structure factors are the same 
Irrespective of whether calculated in cell No.3 or No. A, the remaining 
discussion will refer to cell No, 4. Substituting z^^ = Zj^2 = in 
eqns (76) and (77) shows that for a planar structure the geometrical structure 
factors of peaks (1) and (2) are equal. 

Cell No. A 

p2 (1) = F , {2)=6'4 {[cos2TTX^^cos2Try^^+cos2TTX^2'^°52TTy^2l ■•■ ^^^^ 

[sin2TTX|j^sln27ry^^+sin2TTX^2si"2Try^2l ^ 

The intensities of the observed powder diffraction peaks (l) and (2) 
can be seen in figure '4.28. It can be seen that the intensity of peak (2) 
is practically zero so that the observed powder diffraction peak intensities 
are described by the Inequality: 






Cell Uo.k 



I , (110 + iTD 

obs 




I , (111 + 111) 

ODS 



(87) 



_ - Bsin 9] 

j F^. (no + lll)e -pr- 

or obs '^ 



2 

-Bsin 62 

j F^^ (111 + 112)e 5-2— (88) 

■' obs ^ 



sin0, sin26 



sin92sin202 



w 



here 6, and O^ are the angular positions of peaks (1) and (2), B is the 



temperature factor, X the wavelength and j the multiplicity. 

The multiplicity j is ^ for all of these peaks while on increasing the value 

y V r / Bs in 0\ 1 
of e from peak (l) to peak (2) the temperature factor term lexpt tt — )i 



131. 



decreases by ^% and the Lorentz term (sinSsinZe) increases by a factor of A, 
Thus the observed structure factors obey the relation: 



^bs(^^ » 



F^ (2) 
obs 



(83) 



while from (86) 



r , (1) 



r , (2) 



(86) 



These conclusions are the same for both cells No. 3 and Uo.h. The only 
way of altering equation (86) is to allow non-zero values for z., and z.^ 
in equations (76) and (77) thereby allowing the Dl molecule to tilt out of 
the (001) plane in cell No. 4. This is unlil<ely because in the middle phase 
the planar arrangement of the Dl molecules is well established. 

Nevertheless calculations of the expected diffraction pattern were 
also carried out on the assumption that the arrangement of the Dl molecules 
in the low phase is not planar. The distortion of a planar arrangement can 
be carried out in the following way. If the space group Is C2/c there are 
two molecules in the asymmetric unit (see (1) and (2) in figure 4.27). 
Assuming that the two molecules do not lie in the (OOl) plane, their 
co-ordinates are related by the equations: 



Cell No. A 



x, = c + X, , 



Z2 = ± z, 



These relationships give four methods of distortion as follows 

1 ' ^ 1 ' ^ 1 ' ' ' ^2 * ^2 ' ^2 

(xj, y^ , Zj) , (x^, y2, -^2^ 

(xp y, , -z,), (x2, y2, ^2^ 

(x^ , y^ , -z,), (x2, y2, -^2^ 



(90) 



(91) 
(92) 
(93) 
(94) 



If the space group is Cc there are four molecules In the asymmetric unit and 
these are related by the following additional equations: 



,-^.. 

:.*■." 



I. .-s 



«l^ 
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X- = 



5-Xl , 73 



yi-5 



Z *= + Z . 

3 - 1 



^ = -><,. yi4 -- -yi 



± z 



(95) 
(96) 



These three expressions (90) , (95) and (96) give a total of sixteen methods 
of distortion as follows: 



(xpYpZj), (X2. 72,22^' (>^3. 73.23). (^ij.y^.z/^) 
(xpYpZ^, (x2,y2,Z2^' (x3.y3.23). (^l^*Yl^>~^^^) 



(xpYp-z^, (x2,y2.-22^' (X3, 73,-23), (x^^.y^.-Zi^) 



(97) 
(98) 



(99) 



Of these sixteen the number of non-equt valent methods of distortion was 
found to be ten. These ten modes were investigated for out-of-plane angles 
(0) of 2.5° and 5°. Only one mode of distortion showed any significant 
Improvement namely, 



u + Xp y2 = 5 - Yj 



^2 = '^] 



(100) 



in the space group C2/c. This model was investigated for 0-values in the 
range 0' to 15"; the best agreement, although only slightly better than 
the 0=0° case, was observed for = 7.5" and the powder pattern calculated 
for this case is shown in figure '♦.30. 



The results of this distortion on the calculated intensities of the 
powder peaks of model No. 5 in unit cell No. 4 may be summarized as follows; 
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Peak No. 



Mi 1 ler indices Cel 1 No. 4 



(1) 


(110) + (111) ; 


(2) 


(iTl) + (ITT) ; 


(3) 


(33T) + (332) ; 


w 


(330) + (333) ; 


(7) 


(202) + (200) ; 



(8) 



(002) + (ooT) 



(9) 



(^02) + (O^iO) 



Change in calculated intensity of 
powder peaf< 

little change, 

significantly reduced, 

1 i ttle change, 

little change, 

an unobserved peak which now has a 

finite intensity. 

reduced slightly, the presence of 

this peak and Its neighbour (9) is 

strong evidence of a planar 

structure. In the orthorhombic 

unit cell these are the (002) and 

(020+200) reflections. 

unchanged. 



li.6 Analysis of the intensity transform of the neutron powder pattern of 
the low phase 

Although only powder diffraction patterns were available for the low 
phase of Dl, it was decided to explore whether the application of Fourier 
transform techniques would provide any useful information on the distri- 
bution of the Dl molecules. For an Isotropic powder, amorphous solid or 
liquid a pair-distribution function g(r) may be defined as the probability 
per unit volume of finding a specific atom type at a distance r from another 
specific atom type. For Dl there are three such pair-distribution 
functions; 9i_i(''); probability of an iodine atom at distance r from 
another iodine atom, gr)_n(r); probability of a deuterium atom at distance 
r from another deuterium atom, g- .(r); probability of a deuterium atom at 
distance r from an iodine atom, where g, _(r) = grj_,(f*). 

Applying to Dl the relationships derived by Powles (1973) for compounds 
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containing two types of atoms in equal concentration we can write; 



'♦TTrp{- 



b^[g,_,(r)-l] + bJ[g-_-{r)-l] + 2b,b_ [g, - (r)-l ] 



I ^^l-l 



D'^D-D 



I D'^l-D 



^1 *2b,bjj-.b2 



-} = 



ri 



x2 



Tr/2 



I ., (6)sin[ r ]5 in2ede 

ui - ^ 



(101) 



where b. and b are the neutron scattering lengths of iodine and deuterium 
respectively, p is the number of atoms per unit volume and lfx,(9) is the 
intensity of the neutron powder pattern observed at the diffraction angle 6, 



If b_=0 the LHS of equation (101) becomes 



ATTrp[g,_,(r)-l] = d,_,(r) 



(102) 



giving the pair-distribution function of the iodine sub-lattice. 



Experimentally the powder diffraction pattern can only be observed up 

to a certain diffraction angle 6 . This restriction creates small 

^ max 

oscillations in the Fourier transform of the diffracted intensity v/hich can 
be reduced by the use of so-called window or modification functions M(e) 
(see Wright 197^). The use of these functions worsens slightly the 
resolution of the pair-distribution function g(r) in real space. 



For the neutron powder diffractioii pattern of the low phase of Dl 
recorded on the 'Panda' at X = 1.886a, the maximum value of 9 was kl" . 
To determine the resolution of the intensity transform of these data a 
very simple model of the pair-distribution function was considered, namely 
a 5-function at lA. The theoretical intensity distribution as a function 
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of scattering angle 6 for this simple model is that of a diatomic molecule 
namely; 



1 Yn\ • r^^r.slnei 
I (6) = s in[ J ^J 



^'"^ -J ^/rkvr^5'\ne-. 



(103) 



-r^: 



where r, 



= iS. 



The function 1(0) Is a sine curve whose amplitude slowly decays as G 
increases. Equation (lOl) can now be rewritten to include the window 
function on the RHS. 



6 
f max 



d(r) = 



l6lT 



l(e)M(e)sin[^^^!-!^|i^]sin2ede 



(104) 



and the window function chosen for this work, which decreases slowly in 
value from one at 9 = to zero at = 9 , Is that flue to Lorch (I969), 



M(9) 



s I 



n[ 



sin9 //Tisin9 \ 
max V ■ . } 



'sine 



(105) 



max 



Figure 4.31 and 4.32 show the intensity transforms of the function 

l(8)M(9) [equations (104) and (101)] and the function 1(9), with the maximum 

value of equal to that of the 'Panda' neutron powder diffraction pattern 

I.e. the Fourier transforms of the intensity distribution below an angle of 

47° of the model 6-function g(r) described above with and without the 

window function M(8) (105) respectively. The functions appear as broad 

6-functions centred about r = lA.with spurious small oscillations due to the 

finite value of 9 . For the window function of equation (105) the full 
max 

width of this 6-function at half Its maximum height (FV/HM) is given by the 
expression (Wright 1974) ; 
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^re 4.31 Resolution function of the intensity ^'^^^^^ ("«^^^,%*if °;;' '"^^"°"^ 
of the neutron, powder diffraction pattern of DI recorded on the Panda 

at 4.2K. ^ 

.arc 4.32 Jteaolutlon function of the intensity trans£onn(using ^^'"''^h/p^ndi""^ 
of the neutron powder diffraction pattern of DI recorded on the Panda 
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Ar(FWHM) = T ^—^ — 

max 



1.138 



(106) 



w 



here 9 for the 'Panda' data Is k7^ and X = 1.886^. 
max 



Figure ^.33 shows the Fourier transform of the observed neutron powder 
pattern of the low phase of Dl at A.2K. Figure k.3^ shows the transform 
of the powder pattern calculated for the model No. 5 (figure 4.27). It is 
seen that the two transforms are very similar. The first peak at ] .kA 
corresponds to the separation of the deuterium and iodine atoms within a 
molecule. The peaks at k, 8 and I2X correspond to first, second and 
third nearest Dl molecules In the crystal. 

Figure A. 35 shows the calculated neutron powder diffraction pattern 
of the iodine sublattice and figure 4.36 shows its Fourier transform. The 
vertical lines in figure 4.36 show the positions of the first six iodine 
nearest-neighbour distances calculated from the lattice parameters of 
table 4.2. There are no distances less than 4.18a [this is the separation 
between an atom at (0,0,0) and one at (c,c,0)] and the small oscillations 
below 3A in figure 4.36 are due to the low value of the maximum diffraction 
angle. 

Figure 4.37 shows the neutron powder diffraction pattern of the low 
phase of Dl recorded on 'Panda' at 4.2K [Dl / , x] minus the calculated powder 
diffraction pattern of the Iodine sublattice [1/ , \] . The Fourier 
transform of this function is shown in figure 4.38. Figures 4.39 and 4.40 
show the calculated powder diffraction pattern for the low phase of 
Dl [Dl, , x] based on model No. 5 minus the calculated iodine powder pattern 
[1/ , \] and the Fourier transform of this function. 



The two transforms (figures 4.38 and 4.4o) are the experimentally 
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Igure 4,33 Intensity transform of the neutron powder diffraction pattern of 
DI recorded on the Panda at 4»2K. 

igure 4,34 Intensity transform of a neutron powder dlffractloni pattern of 
DI calculated on the basis of the model showu in figure 4»ig, 
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i'lgure 4,35 Calculated neutron powder diffraction pattern of the iodine sublattice. 

'igure 4. 30 Intensity transform of the calculated neutron powder diffraction 
pattern of the iodine sublattlce showing the first six iodine 
neareat-nftlchbour distances calculated from the lattice parameters 
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Figure 4,37 Neutron powder diffraction pattern of DI recorded on Panda at 4,2K 
minus the calculated neutron powder diffraction pattern of the 
iodine sublattice* 

Figure 4,38 Intensity transform of the neutron powder diffraction pattern of 
DI recorded on Panda at 4,2K minus the calculated neutron powder 
diffraction pattern of the iodine sublattice. 
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figure 4»39 Neutron powder diffraction pattern of DI calculated on the hr.:;±3 
of the model shown in figure 4.19 minus the calculs,ted neutron 
powder diffraction pattern of the iodine sublattlce, 

figure 4.40 Intensity transform of the neutron po^^tJsr dtffrsction pattern of 

DI calculated on the basis of the model shown in li^uro 4, 19 minus; 
the calculated neutron pow<lor diffractioa pattern of the iodine 
' subiattice. 
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derived and the calculated (model No. 5) pair-distribution functions for 
D-l and D-D distances respectively. Equation (lOl) is the transform of Dl 
and equation (102) is the transform of the Iodine sublattice so that 
figures ^.38 and ^.^0 are the function; 



4.rp {bgfgp.pCr)-)] •^2b,bp[gp.,(r)-l]^ 



(107) 



b^ + 2b, b^ + b^ 



The two transforms of figure ^.38 and k.kO differ In the region of 
7^ but otherwise are very similar. The result is not surprising as this 
distance represents the diffraction peak (1) in figures A. 28 and A. 29. 
The difference between the iodine sub-lattice and the calculated Dl lattice 
transforms (figures 4.36 and A.34) can be seen as a broadening of the first 
peak centred at ^A, and represents the arrangement of deuterium atoms about 
the Iodine atoms. Due to the poor resolution in real space (> 1a) and 
because the transform Is only a one-dimensional representation of a three- 
dimensional structure, any direct interpretation of tlie transform in terms 
of a model is unfortunately not possible 



A. 7 Concl us ion 

Neutron powder diffraction patterns of the low phase of Dl at 4.2K 
were recorded with three wavelengths (l.06, 1.886 and 2.66A) and were 
assigned Miller indices on the basis of a large face-centred orthorhombic 
unit cell of 32 molecules. These observed data are reasonably consistent 
with each other and there is no evidence of preferred orientation in the 
samples. The criterion was established that the low phase structure is 
related to the middle phase of Dl by 180° (or possibly 90°) re-orlentat Ion 
of the deuterium atoms about the Iodine atoms. On this basis the space 
group Fdd2(No.'i3) was considered and 32 models were generated. Of these 
only 3 were separate non-equivalent structures giving calculated diffraction 
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patterns quite unlike the observed neutron powder diffraction pattern of 
Dl at k.lK (particularly at low angles). 

The large face-centred orthorhombic unit cell was transformed Into 
6 different C-centred monocllnic unit cells of 16 molecules each. The 
middle phase of Dl v;as considered in these cells and It was found that; 
2 cells supported no monocl i nic symmetry operations, 2 supported space 
groups Cc(No.9) and C2/c(No.l2) and 2 supported space group Cm(No.8). 
By allowing 180° re-orlentat ion of the deuterium atoms, the cells supporting 
space groups Cc and C2/c each produced the same 8 models. Of these only 
5 were separate non-equivalent structures and the calculated diffraction 
pattern of model No. 5 (figure ^.\S) gave reasonable agreement with the 
observed neutron powder diffraction pattern. Also In this unit cell the 
deuterium atoms were allowed to re-orientate through 90" intervals and a 
further 128 models were considered in space group Cc producing one new 
model whose calculated diffraction pattern was similar to the observed 
neutron pattern. 

At this stage the large face-centred orthorhombic unit cell was 
re-considered with 8 molecules in the asymmetric unit and by allowing 180° 
re-orientation of the deuterium atoms produced 128 models. Of these there 
were only 17 non-equivalent structures of which 12 had unit cells of 
32 molecules (these 12 included the 3 structures of space group Cc generated 
earlier). Difffraction patterns calculated on the basis of these 9 structures 
did not agree with the observed neutron powder diffraction pattern. 






The 2 nvonoclinic cells supporting space group Cm were considered and 
produced 32 nodels when allowing the 6 molecules of the asymmetric unit to 
rc--oricnl.nc through \iAj'\ All of thc5e32 rodels were equivalent to n;Gc__^els 
qenerated in the large orthorhombic unit cell and are included in the siii'^niary 
of r,t rucLurt V- considered shown in Table 4.3- 
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The model No. 5 (figure ^.19) was refined by the computer program of 
Rietveld (Appendix one). The analysis of the data is complicated by the 
large number of unresolved Bragg peaks. Several different modes of 
distortion of the structure were introduced but these produced no significant 
improvement in the agreement between the calculated and observed powder 
diffraction patterns which differed from each other in a few low angle peal<s. 

Fourier transformation of the observed intensity distribution is not 
very productive because of the poor resolution in real space and because the 
transformed intensity is only a one-dimensional average of the nuclear 
density function. 






To conclude tills chapter it can only be said that the structure of the 
low phase of Dl is probably very close to that described by the model shown 
in figures A. 19 and ^.27- The difference between the proposed model and 
the true structure is due to some small change in the symmetry, possibly 
involving a larger unit cell than the ones considered. Any further 
investigation of the structure will involve either recording powder patterns 
with much better resolution than is possible at present or else using 
diffraction data collected from a single crystal; if one can be grown. 
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4 molecules 

with 180° 
reorientation 

space group 

Cc 

8 possible 
models 

5 structures 



/i molecules 

with 90° 
reorientation 
space group 

Cc 

128 possible 
models 



6 molecules 

with 180° 
reorientation 

space group 

Cm 

32 possible 
models 



Tabic 4.3 3ui>L-nary ol! the unit cells and spaco-grou-ps considered for the 

structure of the low phase of DI . 
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CHAPTER FIVE 



Neut r on povjder diffraction study of the crystal structure of the 

middle phase of P I 



The neutron powder diffraction pattern of the middle phase of Dl was 
recorded with the Dido ' Curran' powder diff Tactometer described in 
Chapter 2.3(a) and shown in figures 2.6 and 2.7. The sample was contained 
in a thin-walled vitreous silica tube at a constant low temperature in a 
cryostat fitted with vanadium tails to provide isotropic background 
scattering. The outer container of the cryostat was filled with liquid 
nitrogen and the inner container (normally filled with liquid helium) held 
dry air. The equilibrium temperature of this cryostat arrangement wa^ 
about 8^K. At this temperature the heat conducted down the walls of the 
sample tube and the inner helium container from the exterior of the cryostat 
was balanced by the radiation loss to the liquid nitrogen container. The 
neutron powder diffraction pattern recorded at 8AK is shown in figure 5.1. 
The broad 'hump' in the background scattering over the 29 = lO^-lS'' range 
is due to the vitreous silica sample container. 

Using lattice parameters derived from the X-ray powder diffraction 
pattern recorded at 8AK (see table 3.1) the neutron powder peaks v/ere assigned 
Miller indices on the basis of an orthorhombic unit cell containing four 
molecules. The following systematic absences were observed; 
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Figure 5.1 Neutron powder diffraction patt-rn of DI at 84K (middle phase) 
recorded with the Dido Currcin dif f ractonietcr. 
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These conditions indicated a B_~ccntred lattice v;i th an a_- or £-gllde 
in the (OlO) plane. Of the twelve possible space groups satisfying the 
above condition of B-centring, only two were compatible with the full set 
of absences; these were Bba2 (No.M) and Bbcm (N0.6A). 

Since the iodine atoms were known from X-ray data to form a face-centred 
orthorhombic lattice, it was necessary to choose special positions for these 
atoms in the space groups Bba2 (No.Al) and Bbcm (No. 64) which would provide 
such a lattice ('super-special' positions). The analysis and tabulation 
of these positions had already been performed, Sandor (1968). 

If the special positions of these space groups are considered the 
following sites are obtained for the deuterium atoms. 

Bba2 (No. 41) general position 8b (Wyckoff notation) 

Bbcm (No. 64) special position 8f and general position l6g. 

These positions provide either eight or sixteen sites, Implying that 
the four deuterium atoms in the unit cell occupy either two- -or four-fold 
disordered positions. Since the recorded neutron powder diffraction 
pattern is a space average over all possible states of the structure it is 
necessary to apply weighting factors (l and I) to the scattering lengths of 
the deuterium atoms occupying two- or four-fold disordered positions 
respectively. Three models based upon the above atomic arrangements are 
shown in figures 5-2 to 5.4. 



Least-squares refinements using the program ORFLS (described in 
Appendix one) v/ere carried out on these models with the out-of-plane iodine- 

dcutcrium bond angle held constant. it was not possible to vary the 
z-co-ordinate of the deuterium atoms during the least-squares refinement of 
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Figure 5»2 Two-fold 

disordered model for DI 

middle phase, space group 
Bba2 (no. 41). 



Figure 5.3 Four-fold 
disordered model for DI 
middle phase, space group 
Bbcm (no. 64). 
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Figure 5.4 Proposed two- 
y^-^-S^y fold disordered model for 
^\^J DI middle phase, space group 

^ Bbcm (no. 64), 
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the models shown In figures 5.2 and 5.3 as the program returned very large 
values for the z-co~ord! nate (due probably to a singular solution of the 
set of norma? equations). A series of refinements was therefore carried 
out and the out-of-plane angle increased by 5° between each refinement. 

It was found that the nxjdel described by the space group Bbcm gave the 
lowest R-factor. The minimum R-factor occurred when the Dl molecules lay 
In the (00 1) plane in which case the models become equivalent. These 
refinements were carried out with unit weights. Subsequently a full 
least-squares refinement was performed using anisotropic temperature factors 
and the vjeighttng scheme proposed by Cruickshanl< (1953). The following 
parameters were varied by the program; 

(a) Scale factor, 

(b) For the deuterium atom; 

>^,y,^]],^22'^3i and bj2 U = 0, b^^ = b^^ = O) 

(c) For the Iodine atom; 

^11 •'^22 ^""^ "^33 tx=y=z=^0, b.2 = b , = bj, = O) 

The constraints applied to the parameters during the least-squares 
refinement were deduced from the theory of Levy (1956). They were due to 
the Iodine being in a special position at the origin and to the deuterium 
atom lying in a mirror plane. This latter condition constrains the thermal 
ellipsoid tensor to be symmetrical about the (OOl) plane. 



In addition to the atom multiplier of (l) applied to the deuterium 
atom scattering length to compensate for It being in a disordered position, 
a factor of (j;) was applied to the Iodine scattering length to correct for 
the effect of the Iodine atom being in a special position (^a of space group 
Bbcm (M0.6A)), 
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The results of the least-squares refinement are given in table 5.1- 
Figures 5-5 and 5.6 siiow in istograms of the observed intensities and those 
calculated. It was found that if the b^^ and h^^ parameters of iodine 
were allowed to vary simultaneously the fitted values of these two para- 
meters were correlated. Tv;o independent refinements were therefore 
performed to obtain values of these parameters. 
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Dl middle phase al 8i(K 

Lattice parameters 

a = 5.979 1 .OO7S 

b " S.OZh + .007^ 

c = 6.517 1 -009^ 

ct - e = Y - 90° 

Space group Bbcm {No. 64); number of molecufes in the unit cell is h. 

Co-ordinates of 8f special position in space group Bbcm (NcS^) {centro- 

symmetrlc) 



X+5 



g-X 



y 

y 



Y+l 



yH 





I 



+ [(0,0,0); (1,0,1)] 



Parameter/ 'i 

Iodine atoms , positional parameters: x=y=0 
thermal parameters: 



scattering length: 



b.p.025+.002 
b^j'^.OSS+.OOZ ^ 
b ,=,030+. 0007 

5.2 fm. 



deternined independently 
due to correlation 



Deuteri urn a toms^, positional parameters: x = -.171 + . 001 

y = .169+. 001 
z = 

thermal parameters: '^11 "^ .O56+. 001 

^22 = •039+.0015 

b,3 = .04ii+.001 

b,2 = .007+-002 

^23 = h3 = ° 
scattering length: 6.7 fm. 



Unit weight R-factor = 6.7% 
Weighted R-factor = 1.3% 

Refinement carried out over 5 cycles on intensities 

Bond length (uncorrected) = 1.44+ .01^ . , ^ , ,,^r, 

Bond length (corrected for librational effect) = 1.49+ .01)^ (Cruickshank 6;>6 

TABLE 5.1 Structural parameters of Dl at 84K (middle phase) 
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Figure 5.5 Observed intensities of neutron powder peaks for DI 
middle phase recorded at e4K with Dido Curran. 
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Figure 5,e CKlculated InirervKitics oi neutron powder peaks for two-fold 

disordered jfiOdel of IJI islddle phase, space group Bbcra (no.G4) 
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CflAPTER SIX 



Neut ro n powd&r diffraction study of the, crystal s tructure of tj-ie 

high phase of D! 



The neutron powder diffraction pattern of the high phase of .Dl was 
recorded with the Dido 'Curran' dt f f ractometer as described in section 2.3a' 
The sample was contained in a vanadium cryostat as in the middle phase 
experiments (chapter five) and maintained at a temperature of 137K by a low 
power electrical lieater. The ernf from a thermocouple mounted on the 
samp5e was monitored by a chart- recorder which controlled the power to the 
electrical hc.ijter by means of a microswitch. 

The neutron povjder diffraction pattern of the high phase of D! 
recorded at i37K is shown in figure 6.1. Using the lattice parameters 
derived ftcni the X-ray powder photographs (see table 3-l)t the neutron 
powder peaks v/ere assigned Hilier indices on the basis of a cubic unit celt 
containing four molecules. The systematic absences In the Miller indices 
we re ; 



hkl: h+l<, k+1, (1+h) = 2n 



(109) 



indicating an fee lattice. 



There are eleven face-centred cubic space groups of which five are 
compatible with these absences. 

F23 (Ho. 196), Fm.3 (K'o.202), F^32 (No. 209), F"^3m (No. 2)6) and 
Fm3m (No. 225). 



An ordered cubic high phase based on the space group Fm3m can be easily 
disregarded. Since tisere are four molecules per unit ct " 1 , In an ordered 
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Figure 6.1 Neutron ponder diffraction pattern of DI at 137K (high phase) 
recorded with the Dido Curran dj.f fractometer. 
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structure the deuterium atoms must occupy four-fold special posittonr. 
This condition leads to two possible models. 

6. 1 Zinc sulphide-type model space group Fm3m (No. 225) 
Iodine atoms in special position ha (Wycl<off notation), 

(0,0,0), (o.hl), iho,l). (^J,o) 

Deuterium atoms in special position he (Wycl<off notation). 
(ij5,n), (i,^,u), (ii.s.i;), (^,^,^) 
This gives a D-l bond length of 2.82^. 

6 . 2 Sodium chori de-type mo de] space gr oLi p. Fnjjn i (No. 225) 
Iodine atoms in specl^jt position ha (Wyckoff notation). 

(0,0.0), (o,.},i), ihoA). (^,^.0) 

Deuterium atoms in special position ^b (Wyckoff notation). 
ihhl), (1,0,0), (O.LO), (0,0, i) 
This gives a D- I bond length of 5-&5A. 



tn the gas phase tiie D-l bond length was found to be 1 . 609A from micro- 
v;ave nicasurenients by Cowan and Gordy (1958) and in the middle phase it v;as 
found to be 1 Ji^S from neutron diffraction measurements (Chopter five). 
Since these two models give unacceptable values for the D-l bond length they 
can be disregarded. Examination of the symmetry operations of the five 
space groups listed earlier shows that there are special positions of these 
space groups v^hich describe disordered structures. Four of these will be 
discussed and the special case of a freely-rotating or randomly disordered 
model. As some of these space groups contain the lower symmetry groups as 
sub-groups, the special positions of the disordered structures will be des- 
cribed in terms of the highest symmetry space groups; Fm3 (No. 202), 
F^r:-:^ (^Jo.2le) and Fr^;;-ii (No, 225). 



6 ■ 3 Four-fold disordere d modet spgc e group F"^3m (No. 216 ) 

Iodine atoms in special position ^ia (Wyckoff notation). 

(0,0,0), (0,^,^), iho,l). ilA.o) 

Deuterium atoms In special position ]6e (Wycl<:off notation). 
{x,x,x)» (x,x,x,), (x,x,x), {x,x",x) 
+ [(0,0,0); (0,^,i); (i,0,|); (^,^,0)] 

(>A SIx-fo _ 1d disordered model sp^CJ^group Fm3m (No. _ 225) 

Iodine atoms in special position ^a (V/yckoff notation). 

(0.0,0), (0,1,1), (l,0,i), (^,1,0) 

Deuterium atoms in special position 2he (Wyckoff notation). 

(x,0,0), (0,x,0), (0,0, x) 

(x,0,0), (0,x,0), (0,0, x) 

+ [(0.0,0); (0,^,^); (1,0,1); (1,1,0)] 

This model represents a kind of close-packed structure which would 
occur if the hydrogen bonds were very v:eak. 

6 . 5 Eight-f Q^!_d_dj^sorde red model space group Fm3 (No. 202) 

Iodine atoms in special position ha (Wyckoff notation). 

(cr,o,o), io.lA), (i,0,i), (^,^,0) 

Deuterium atoms in special position 32f (Wyckoff notation). 

(x,x,x), (x,x,x) , (x,x,x) , (x,x,x) 

(x x x) , {7 x x) , (x X x) , (x x x) 

+ [(0,0,0); (o,LD; (1,0,1); (l,Lo)) 



6.6 Twelve-fold di so rdered^mode_l_ s pace gro up Fm3m {No.2 23) 
Iodine atoms In special position ka (V^ycloff notation) . 

(0,0,0), (o.ij), (1,0,1), ihho) 
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Deuterium atoms in special position h&U (Wyckoff notation) 

{x,x,0), (x,0,x), (0,x,x) 

(x",x,0), (x,0,x), {0,x,x) 

(x,x,0), (x,0,x), iO,x,x) 

(x,x,0), (x,0,x), {0,x,x) 

+ [(0,0,0); (0,^J); il ,0 A) ; ilA,o)] 

This model represents the case when hydrogen-bonding between neigh- 
bouring iodine atoms Is strong. Since each iodine atom has twelve nearest- 
neighbours the deuterium atom can lie along each of the twelve nearest- 
neighbour bonds as It 'flips* between different disordered sites. 

6-7 Freely-rotat f ng or randomjy-d rsord ered model 

This model has been suggested by several workers although, as discussed 
in section 1.6, It is sha-;n to be unlikely. Since the neutron diffraction 
pattern records on average structure over space and time the two forms of 
the model, freely-rotating or randomly-disordered, cannot be d i st InguisHed 
(by neutron diffraction) and are described by the scime mathematical 
expressions. A set of structure factors was calculated by the following 
method. * 

The position of the deuterium atom relative to the iodine atom was 
described by a Gaussian radial distribution function, p(r); 



h"™^'| 



(r)dr = ^ [ -1^1 ]dr 



(110) 



where r Is the distance of the deuterium atom from the iodine atom, 1 Is the 
mean intramolecular bond length, a Is the width of the Gaussian and 
represents the mean amplitude of vibration along the bond and Z Is a 
normalising constant. 
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The mean bond length of the molecule t was taken to be 1 . ^tjA and the 
half -width of the Gaussian cf was taken as 0.08A, the mean vibrational 
amplitude of the deuterium atom calculated from spectroscopic data. Tlic 
normalizing constant Z was calculated from the equation; 



■J 

hi^y^'expi- ^^4^ ]dr 



cm) 



The structure factor was given by the expression: 



2 

F'(hkl) = /,b,exp("^-^'iJ^)+4b^ 



H-rrr p(r)sin( ^ ) ^^ (112) 



This expression was evaluated and integrated numerically using Simpson's 



rule. 



6.8 Comparison of observed structure factors with structure factors ca lcu- 
lated from dis ordered mode ls 

Projections of the seven models described above onto the (OO)) plane 
are siiov/n in figures 6.2 to 6.8. Solid circles Indicate atoms at z *= 
and open circles atoms at z = |, Histograms of calculated structure 
factors are shown in figures 6.3 to 6.13. The structure factors F' {hk}) 
contain the expression for thermal motion and are defined in chapter three 
equation (30}- l^ the histogram of the observed structure factors shown 
In figure (>.]h there are three places where the observed structure factor as 
a function of 2B deviates appreciably from a monotonical 1 y decreasing 
function. in particular these are: 

F' (?.22) > F' (113) 
F'(i33) => F'COO'4) 
K' (;-2'i) > F' (U2'i) 
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Figure G.2 Zinc sulphide 
model for DI high phase 
space eroup Fm^ni (no. 225) 




Figure 6.3 Sodium 
chloride model for 
DI high phase, space 
group Fm3m (no. 225) 




Figure 6.4 Four-fold 
disordered model for DI 
high phase, space group 
F'TSm (no, 216). 



^2, 




Figure 6.5 Six-fold 
disordered model for DI 
high phase, space Rroup 
Fm3m (no, 225), 




Figure 6.6 Eight-fold 
disordered model for DI 
high phase, space group 
Fiij3 (u6. 202). 
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Figure G.7 Twelve-fold 
disordered model for DI 
high phase, space {jroup 
Fm3m (no. 225), 





Figure 6,8 Freely-rotating 
model for DI high phase. 
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Figure 6.9 Calculated structure factors for four-fold disordered model 



of DI high phase, space group F43m (no. 21b). 
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Figure G.IO Calculated structure factors for i.lx-fold disordered model 
of DI high phasn, space group Fm:iT^ (no. 225 ), 
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Figure 6.11 Calculated structure factors for eight-fold disordered model 
of DI high phase, space group Fni3 (no* 202), 
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Figure 6.12 Calculated structure factors for freely-rotating model 

of DI high phase. 
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Figure 6.13 Calculated structure factors for twelve-fold disordered model 
of DI high phase, space group Fm3 in (no, 225). 



u 
o 

u 

C8 

o ^ 
u 

3 

u 

n 

CO 



(111) 

(002) 



DI 137K (HIGH PHASE) 



OBSERVED STRUCTURE FACTORS 



(022) 



(024)(224) 



^^4-'^ (133) 



(113) 



(004) 



(135) 
(044) 



10 



20 



30 40 

29 / DEGREES 



50 



60 



70 



Figure 6.14 Observed and normalised structure factors for DI high phase 

recorded on Dido Curran at 137K. 
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These relationships also appear in the histogram calculated on the 
basis of the twelve-fold disordered model (see figures 6.7 and 6.]k) but 
are absent in the histogram calculated on the basis of the randomly- 
disordered model (see figures 6.8 and 6.12) which shows an angular variation 
of intensity of a type normally associated with an anK>rphous material - 
Least-squares refinements were carried out on all but the last of the 
disordered models using peak intensities with unit weights as data. The 
results are shown in Table 6.1. 



Model 


Unit-weight R-factor (intensity) 


ii-fold 

6-fold 

8-fold 

12"fold 


15.9? 
36.8^ 
15.8? 



TABLE 6.1. R-factors obtained in the least-squares refinement 
of four disordered models proposed for the high 
phase of Dl . 



The twelve-fold disordered model gave the best fit and the results of 
a least-squares refinement of this model using anisotropic temperature 
factors are shown in table 6.2. Because the iodine and deuterium atoms 
were in special positions it was necessary to apply Levy's rule to derive 
relationships between the temperature factors (Levy 1956). The following 
parameters were varied during a least-squares refinement based on the 
diffracted intensity. 

(a) Scale factor, 

(b) For tlie deuterium atom; 

x(x=^y) , b|](b..=b ), b , and h^^ (^=0 . '^j 3'^'^23'"'^^ 
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(c) For the iodine atom; 



(x=--y=--z=0 , ^]2^^2Z^^\3^^^ 



The least-squares refinement of the 12-fold disordered model appeared 
insensitive to the value of the parameter b,^(D). The computer program was 
altered to increase the value of b--(D) from zero to .06 In steps of .005. 
This produced no change In any of the fitting or structural parameters. 
A refinement was also carried out using only the single peaks of the cubic 
pattern with the unmodified version of the ORFLS program v/hi ch refined upon 
the structure factor [p'j and not the intensity. This produced a worse fit 
because of the smaller number of data and also gave a singular solution when 
varying b ,(D). The weighting scheme used in these calculations was the 
one proposed by Cruickshank (1963). 
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PI hig h phase a t I37K 

La 1 1 i ce pa ramete r s 

a = b = c = 6.285 + .018 
a = 6 = Y = 90° 

Space group Fm3ni (No. 225); number of molecules in the unit cell is 4. 

Co-ordinates of 48h special position in space group Fm3m (No. 225). 

(x.x.O), (x,x,0), (x,x,0), (x,7,0) 
(x.O.x), (x,0,x), (x,0,x), (x,0,x) 
(0,x,x), (0,x,x), {0,x,x), (0,x,x) 

+ [{0.0,0); {0,1,1); (5,0, J); (^,^,0)] 



Param eters 

lodi ne atoms , positional parameters 
thermal parameters: 

scattering length: 



X = y = z 



^2 = ^23 
5.2 fm. 



33 



= .0177 + .0002 



=b|3=0 



Deuteri um atoms , positional parameters: 



scattering length: 

Unit weight R-factor = 6.3^ 
Weighted R-factor = 0.]% 
Bond length (uncorrected) = 1.42A 



X = y = . 1578 + .0002 
z = 



thermal parameters: b, , = b = .0^69 + .0001 



b - > .01 

b,2 = -.0289 + .0006 

b,3 = b23 = 

6.7 fm. 



TABLE 6.2. Structural parameters of Dl at 137K (high phase) 
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CHAPTER st:vEr: 

Neutron powder diffraction study of the tv;o phase 
transitions in solid Dl 



7- 1 The lower phase transition 

The analysis of the neutron powder diffraction experiments described 
in the preceding chapters revealed that the crystal structure of Dl is 
ordered monocltnic at ^,2K (low phase) and two-fold disordered orthorhombi c 
at ShK {middle phase). The specific heat measurements of Clusius and 
Wolf (19V) (see Chapter One) suggest that the transition from the low phase 
to the middle phase is a gradual one extending over a temperature range of 
A5-8OK. The progress of this transition can be followed by recording 
the gradual disappearance of the (110, 1 iT) powder peak of unit cell No.^ 
which appears only in the low phase. Hereafter all monoclinic Miller 
indices refer to unit cell No.A. 



To compensate for any variation In the measuring system with time 
(such as sample n-.ovemcnt while filling the cryostat with liquid helium or 
change in the neutron detection efficiency) the recording was extended to a 
nearby 'reference' powder peak, whose intensity did not change appreciably 
over the transition region. The most convenient choice for the reference 
peak appeared to be the (220, 222) powder peak, which at the time of the 
phase transition experiment was thought to be a single peak. Subsequent 
experiments witfi higher angular resolution showed that this peak also con- 
tained contributions from the (312", 31T) and (130, 1 3T) powder peaks which 
appear only in the low phase. This contribution was rather small to have 
iny i.pp rec I i^Ll i; viffcct on th.: results. The geoii^et r i cal structure foctcr 
for the proposed ntudel of the low phase (see figure h.]S) is; 



F{hkl) ■-- -8sin2jr{hx+l2)sin2irky (h+k=2n, l = 2n+l) 
F(hk]) = 8cos2TT(hx+lz)cos2Trky (h+k=2n, 1==2n) 
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(113) 
(11^) 



Since the observed neutron powder diffraction pattern Is a time average 
taken overthewhole sample, the molecules are "seen" by the neutrons in all 
permissible equilibrium orientations. If the fractional number of 
deuterium atoms in the 180° "flipped" position at temperature T is denoted 
by p(T), statistically averaged intensities of the low phase and reference 
peaks scanned over the lower transition region can be written in the form: 

F^(llT) + F^(llO) = 6/4 b^[l-2p(T)]2[5in22TT{x+y)+sin^2Tr{x-y)] (II5) 

F^(3lT) + F^(3lT) + F^(13T) + F^(130) = 

ek bjj[l-2p(T)] [sin^2TT(3x-y)+sin^2TT(3x+y)+2sin^2TT(x-3y)] (116) 



F^(220) + F^(222) = 51 2[b +b cosifTixcos^try) ]' 



(117) 



where bj and b^ are the iodine and deuterium thermal neutron coherent 
scattering lengths. 

To simplify the above expressions the origin of the lattice has been 
moved to an iodine atom so that the two equilibrium positions of the 
deuterium atoms In the middle phase are at {x,y,0) and (x,y,0). The 
deuterium temperature factors were assumed to be equal in both orientations 
of the molecule. 



The fraction p(T) increases from zero in the ordered low phase to O.5 
in the two-fold disordered middle phase. In analysing the results it will 
be convenient to use a new function q(T) defined as: 



q(T) = l-2p(T) - 




[F^dlO) + F^dlT)] , 

I 1 observed 



[F^(llO) + F^dlT)] . . 
1 T ordered 



(116) 
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Becouse of the gradual nature of the phase transition, it is not 

2 

possible to measure F_(hkl) , , directly over the transition region. 
'^ T ordered 

Nevertheless, an approximation to equation (118) can be made by assuming 
that the following ratio Is constant over the temperature range studied; 



f[fl{]\0} + F^dlT)] . , 
T T ordered 

ordered J [F^(220) + F^ ( 22"^ ]^^,^^^^ 



R(T) 



(119) 



This approximation is acceptable for the following reasons; 

a) The Lorentz-factor ratio R. (T) changes as the lattice 
parameters increase with temperature: 



Rl(t) = 



sine.j.(nT)stn2e^(llT) sine^(llO)sin20^(llO) 
sine^(222)sin2e^(227) s ine^(220)5 i n2e^(220) 



(120) 



This change is very slight since at all temperatures in the low phase: 

2sine^(llT) = 2sln6^(no) = sin9-|.(22T) = slne^(220) (121) 

so that from figure 3-5, \W = .389 and R|_(75) = -388. 

b) The ratio of the Isotropic temperature factors for the two peaks 
R^(T) wh< 



lere; 



R^-d) = 



exp [-Bs i n^e^ ( 1 iT) ] exp [-Bs I n^G^ (l 1 0) ] 



. 2 



. 2r 



exp [-Bs I n e^ (222) ] exp [-Bs i n 6^ (220) ] 



using equation (121) 



(122) 



R^(T) - exp[Bsin^G^(nT)] = exp [Bs i n^B^ (l lO) ] 



(123) 



changes very slightly over the temperature region investigated. At ^.2K 
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this ratio is 1.009 and at 84K Is 1.025 representing a 2l°/, change which is 
small compared with the standard deviations of the measured ratios. 

According to equations (115) and (116), the structure factors of the 
Bragg reflections; {110. llT) , (3l'2, 3lT) and (130, I3T) are proportional 
to [l-2p(T)]. The (220, 222) structure factor is an even function of x and 
y and is unchanged by the 'flipping' of the Dl molecules. The ratio of 
peak intensities, 



R,(T) = 



I Oio+iiT)_ 

I (220+222+3 1I+31T+1 30+ 13T) 



(I2i*) 



was recorded at different temperatures over the lower transition region 
using the Dido 'Curran' dl f f ractometer. 



At^.2K the structure was assumed to be ordered, i.e. p(T) was zero. 
At this temperature the relative peak intensi^ties of the I (l 10+1 iT) and 
I (312+31 1 + I3O+I3T) powder peaks were measured using the better angular 
resolution available with the Pluto 'Panda' di ff ractometer using a neutron 
wavelength of 2.66A. It was assumed that the ratio of the peak Intensities 
did not vary with temperature In the low phase, so that if; 



I (3I2+3II + I3O+I3T) 
r(^l) = 1(1 10+1 iT) 



(125) 



Is measured at ^.2Kwlth neutrons of 2.668(X.), then for the Dido 'Curran' 
diff ractometer at 1. 06^(^2) this ratio when corrected for the Lorentz-factor 
Is; 



rix^) 



r(x,) 



[Ad^dlO) - \]][hd^i3]2) - X^] 
[^d^llO) - A^][Ad^312) - \h 



(126) 



where d(llo) etc. are the lattice spacings for these reflections 
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The ratio; 



R(T) 



[i(no) + \{\]T)] 



observed 



(127) 



[1(220) + 1(222)] 



observed 



can be calculated from the observed intensity ratio Rj (T) equation (12^) and 
r(1.06) equation (126) as follows; 



R(T) 




R,(T) 
1 - R| (T) r(1.06) 



(128) 



Comparing equation (127) with (119) and (118) and using the approximations 
discussed above we can write; 



q(T) = l-2p(T) « R(T) 



(129) 



The anomalous X-peak observed in the specific heat measurements of 
Clusius and Wolf (19^7) over the lower phase transition (see chapter one, 
figure 1.1) can be attributed to the change in conf igurat ional entropy as a 
result of molecular re-orientation. 



From the definition of specific heat used in statistical mechanics; 

T 



t T 



Cp'"'(T)dT _ 
T 



(130) 



T is the temperature where the phase transition begins and C (T) is 
the 'anomalous' specific heat at temperature T. The 'anomalous' specific 
heat is defined as the amount in excess of that due to vibrational motion 
(shown as a dotted line in figure i .]) . 
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A1sc; 



S = k tn .Q 
c 



(131) 



where k Is the Boltzmann constant, S the molar conf i gurat ional entropy, 
^ the number of v/ays of distributing p(T)N reversed molecules on N lattice 
sites and N is the Avogadro constant. It can be shown that for a tv.-o-fold 
disordered system In which the two equilibrium orientations are equally 
probable; 



N! 



"^P^"^^^ ~ {Np(T))! {Nil-p(T)]n 



Using Stirling's approximation for large values of N; 



(132) 



In N! = NlnN-N 



(133) 



so that 



^^Sn^ = -kN{p(T)lnp(T) + [l-p(T)l]n[l-p(T)] 



p(T^)lnp(T^)-[l-p(T^)]ln[l-p(T^)]> 



{13M 



The ratio R. (T) and hence R(T) was determined experimentally at severa' 
temperatures over the lower phase transition of Dl from ^5K to 80K using the 
'Curran' neutron powder dl f f ractometer shown in figures 2.6 and 2.7. 



The majority of the measurements were taken at temperatures where 
R| (T) (the ratio of peak intensities) was expected to change appreciably. 
The sample was contained in a vitreous silica tube contained in a liquid 
helium cryostat as shown in figure 2,10. The temperature was controlled by 
the circuit described in Appendix two. Each measurement took approximately 
eight hours and after each temperature change the sample was allowed two 
hours to reach thermal equilibrium. Measurements were- first made with 
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increasing temperature until the sample reached the middle phase, then with 
decreasing temperature going back into the low phase and finally again with 
increasing temperature to close the cycle with some overlap and interpolating 
between the earlier measurements to check the temperature measuring system 
for dri ft. 

The phase transition appeared to be complete at 78K when recorded with 
increasing temperature, whereas the specific heat measurements Indicated a 
transition temperature of 77.3K, This discrepancy is probably due to the 
radial temperature gradient in the sample used in the neutron diffraction 
work. Differences in the isotopic purities of the samples used in the 
calorimetric and diffraction studies may also have contributed to it. The 
radial and axial gradients were measured later as described in Appendix tv/o 
and both were found to be approximately 0.5°K- A small thermal hysteresis 
of approximately 0.5°K was recorded. It is not possible to give a more 
accurate value for the width of the hysteresis due to the magnitude of the 
thermal gradients in the sample (^O-S^K). 

Figure 7.1 shov/s the temperature dependence of the ratio R(T) , 
measured with both Increasing and decreasing temperatures, and table 7.1 
lists the observed peak Intensities as well as the R(T) ratios and the 
order parameters q (T) derived from them. The expression for the configura- 
tional entropy equation (130) was evaluated graphically using the specific 
heat data of Clusius and Wolf (19^7). The anomalous specific heat and the 
order parameters derived from these data are listed in table 7.2. 



Clusius and Wolf measured the integrated heat input used in their 
measurement of the specific heat. This enabled values of the entropy change 
over the lov;er transition region to be calculated. Values of p(T) and 
hence q(T) were obtained from eqn.dS'^) by an interpolation technique. 
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Figure 7.1 Temperature dependence of the Intensity ratio 



u(t) = / I(li0)+I(lll) 
l(220)+l{222) 



measured over the lower phase transition of DI with 
increasing (x) and decreasing (o) temperatures , 



178. 



T^K 


Di rect ion 


1 (110+1 iT) 


1 (2 20+2 22) 


R(T) 


q(T) 


4.2 ± 1.0 


Increasing 


2383 


9935 


0.49 ± 0.03 


1 .0 ±0.06 


A7.8 ± 1.0 


• I 


5297 


22230 


0.49 ± 0.03 


1.0 ± 0.06 


61.8 ± 1.0 


1 1 


5368 


24066 


0.47 ± 0.03 


0.96 ± 0.06 


67.5 ± 1.0 


II 


6510 


32873 


0.45 ± 0.03 


0.92 ± 0.06 


68.5 ± 1.0 


11 


4675 


23535 


0.45 ± 0.03 


0'.92 ± 0.06 


71.8 + 0.5 


II 


3815 


22188 


0.42 ± 0.03 


0.86 ± 0.06 


7A.1 + 0.5 


II 


4725 


30175 


0.39 ± 0.03 


0.80 ± 0.06 


75.1 ± 0.5 


II 


4290 


32336 


0.36 ± 0.03 


0.735± 0,06 


76.2 + 0.5 


It 


3138 


31784 


0.31 ± 0.04 


0.63 ± 0.08 


76.2 ± 0.5 


II 


1047 


12077 


0.29 ± 0.04 


0.59 ± 0.08 


77.5 ± 0.5 


M 





32829 


0.0 ± 0.05 


0.0 ± 0.10 


76.9 ± 0.5 


Decreasing 


1240 


29615 


0.2 ± 0.04 




75.6 ± 0.5 


1 1 


3090 


28011 


0.33 ± 0.03 




7^.5 ± 0.5 


II 


3080 


24568 ' 


0.35 ± 0.03 




73.4 ± 0.5 


II 


4585 


31650 


0.38 ± 0.03 




69.6 ± 1.0 


II 


5170 


30920 


0.41 ± 0.03 




57.1 ± I.O 


" 


6648 


28900 


0.48 ± 0.03 




47.8 ± 1.0 


II 


6235 


26324 


0.49 ± 0.03 





TABLE 7.1. Intensities of neutron powder peaks measured over the lower 

phase transition of DI with increasing and decreasing tempera- 
ture and the R(T) and q(T) values derived from these 
measurements. 
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T°K 


r(T) 


p(T) 


q(T) 


VK 


Cp(T) 


p(T) 


q(T) 

1 


kk.O 


0.0 


.015 


.970 


66.5 


1.10 


.030 


.940 


55.0 


0.0 


.015 


.970 


67.0 


1.25 


,031 


.938 


56.0 


0.20 


.015 


.970 


67.5 


1.35 


.033 


.934 


56.5 


0.20 


.016 


.968 


68.0 


1.^5 


.034 


.932 


57.0 


0.25 


.016 


.968 


68.5 


1.65 


.036 


.928 


57.5 


0.30 


.016 


.968 


69.0 


1.90 


.038 


.924 


58.0 


0.30 


.017 


.966 


69.5 


2.20 


.041 


.918 


58.5 


0.^0 


.017 


.966 


70.0 


2./*5 


.044 


.912 


59.0 


0.^0 


.017 


.966 


70.5 


2.75 


.047 


.906 


59.5 


0.^0 


.018 


.964 


71.0 


3.20 


.051 


.898 


60.0 


0.50 


.018 


.96if 


71.5 


3.50 


.055 


.890 


60.5 


0.55 


.019 


.962 


72.0 


3.90 


.060 


.880 


61 .0 


0.55 


.019 


.962 


72.5 


A. 25 


.065 


.870 


61.5 


0.60 


.020 


.960 


73.0 


i».70 


.072 


.856 


62.0 


0.65 


.021 


.958 


73.5 


5.25 


.079 


.842 


62.5 


0.70 


.022 


.956 


7A.0 


6.85 


.088 


.824 


63.0 


0.75 


.022 


.956 


74.5 


10.80 


.105 


.790 


63.5 


0.80 


.023 


.954 


75.0 


12,80 


.126 


.748 


6^1.0 


0.85 


.02^ 


.952 


75.5 


14.00 


.149 


.702 


6^.5 


0.90 


.025 


.950 


76.0 


17.30 


.185 


.630 


65.0 


0.95 


.026 


.9^8 


76.5 


19.04 


.232 


.536 


65.5 


1.00 


.027 


.9^*6 


77.0 


29.00 


.336 


.428 


66.0 


1.05 


.028 


.Skh 











TABLE 7.2. Temperature dependence of the anomalous specific heat over the 
lower transition of Dl. The values listed were calculated from 
the data of Clusius and Wolf (1947). 
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Figure 7.2 bhows Lhe order parafiieters derived from the neutron diffraction 
experiments as crosses and those from the calorlmetr ic data as a continuous 
line. Only data recorded with increasing temperature are shown in this 
figure as the specific heat measurements were made with increasing 
temperature. 

A value of 0.015 was chosen for p{T ) , effectively the residual dis- 
order parameter at hSK, as this provided the best overall fit between order 
parameters derived from the calorimetric and the neutron diffraction data. 

The reason for a non-zero residual disorder p(T ) at A5K Is probably a 
combination of several effects; statistical errors in R(T) , temperature 
gradients, uncertainty In the magnitude of the components of the reference 
peak (UTetc.) and possibly also some short-range order over the transition 
region. 

Short-range order implies the formation of groups of molecules whose 
re-orlentatlon Is coupled over short distances. The expression derived for 
(I, equation (132), is then not applicable as In this case the two 
equilibrium orientations are not equally probable. As the combination of 
the effects of temperature gradients, statistical errors, etc. appears to be 
less than ].5%, short-range order is unlikely to be important in the lower 
phase transition of Dl. 



7.2 The upper phase transition 

in the upper transition the crystal structure of Dl changes from two- 
fold disordered orthorhombic {middle phase) to twelve-fold disordered cubic 
(hicjh phase). Specific heat measurements suggest that the transition takes 

plrjcf. over the Leriipcrature range of 100-130K. 
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Figure 7.2 Order parr^meter q(T) derived from neutron diffraction 

experimonta ( + ) and specific heat measurements (sc"" Vd line) 
over the lower phase transition region of Dl , 
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When passing from the middle to the high phase (see figure 3-5), the 

gradual disappearance of the (002) ^, , , . peal< of the middle phase can 
^ ^^ ortnornombic 

be easily monitored since It occurs at a lower angle than the {002} , . 
' cub I c 

peak of the high phase. The (111) peak however hardly changes either in 
position or in Intensity over the upper transition region. Accordingly 
the neutron diffraction study of the upper phase transition of Dl was carried 
out by measuring the Intensity of the disappearing (^^^^oj-t-horhomb i c '^^'^'^ 
as a function of temperature. The Intensity of the nearly constant (111) 
peak was measured as a reference to correct for changes In counter 
efficiency, etc. 

An occupation number, p' (T) at a temperature T, can be defined for the 
high phase as the probability of a deuterium atom being in one of the ten 
equilibrium sites not occupied In the middle phase. It follows from its 
definition that p'(Tj will lie In the range 0_<p'(T) <_ 12. 

It is possible to consider two methods In which the upper phase transi- 
tion occurs; 

a) At low values of p'(T), when the structure Is orthorhombic, the ratio 

of the intensity of the (002) ^, . , . peak at temperature T (when 

orthorhombic 

P'(T) deuterium molecules are in the twelve-fold disordered position) to the 

intensity of the (002) .... (when the deuterium molecules are all in 
' orthorhombic 

the two-fold disordered position) is, 



l.,(002) ^, . ,. ^b, + iib_[l-p'(T)(lO-8cos'4TTx)] 
T orthorhombic , I D -, 

T"(oo2) , : TT" " ^ Ab. + Ab, ^ 



(135) 



orthorhombic 



I 



X is the fractional deuterium co-ordinate and b., b- are tlie coherent neutron 
scattering lengths of the iodine and deuterium atoms. 
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Inserting numerical values gives; 



1^.(002) ,, . .. 2 

T orthorhombic n t o i /-r\ i 

\A002} ,, , ,. ^ f^"7.8p (T)] 
orthorhombic 



(136) 



b) The Dl powder samples used in the studies of the phase transitions 

were not entirely free of temperature gradients (see Appendix two). As the 

temperature and hence the value of p' (T) increases, the temperature 

gradients in the Dl powder sample cause the neutron powder diffraction to 

appear as a mixture of the orthorhombic middle phase and the cubic high 

phase. The intensity of the (002) ^. , , . peak decreases in direct 

orthornombi c 

proportion to the number of unit cells or crystallites which have the cubic 

(high phase) crystal structure. Thus the intensity of the (002) , , . . 

orthorhombic 

peak decreases according to the equation; 



1^.(002) ^, , ,. 
J orthorhomb i c 

orthorhombic 



= l-12p'(T) 



(137) 



Following the treatment of the lower phase transition, it can be shown 
that the number of ways of arranging N molecules so that a fraction p' (T) 
occupy the ten lattice sites (additional to the two already occupied in the 
two-fold disordered orthorhombic middle phase) with equal probability Is 
given by; 

n! 



n'[p'(T)] = 



(^I^i^H;<I^l^l}{p'(T)N!} 



10 



(138) 



where N is the Avogadro constant and < p' (T) < 12 



The two limits to equation (138) are the fully twelve-fold disordered 

high phrr:>o; 



n'C^\2) = 



HI 






(139) 



\Bh. 



and the tv;o-fold disordered low phase; 



fi'(O) 



N! 



W7 



(140) 



which Is in agreement with the expression for the fully tvJO-fold disordered 
middle phase used to describe the lower phase transition equation (132); 



Qil) = 



n! 



(141) 



Using the expressions (130), (133) and (138) it can be shown that for 
the upper phase transition; 



C"(T)dT = 

M 



kN 



{[l-lOp'(T)]JLn [ ^'^°P'^^^ ] + 10p'(T)Jlnp'(T) 



- [l-10p'(T )]£n [JblOEltlol-] - 10p'(Tjllnp'(Tj} (142) 



p'(T ) is the fractional number of deuterium atoms occupying the ten non- 
orthorhombic lattice sites at the temperature T^ and k is the Boltzmann 
constant. 



The ratio of the intensities of the (002)^^^.^^^^^^^^,^ and the (111) 
powder peaks was measured as a function of temperature on the Dido 
'Badger II' neutron powder di f f ractometer. The liquid nitrogen cryostat 
and the vitreous silica sample tube used in these experiments were already 
described in Chapter two. The (002)^^^^^^^^^^^,.^ and (002}^^^^,^ peaks over- 
lapped slightly in the powder diffraction pattern and were resolved by the 
use of the computer program MPR, Johnson (1970). 



The temperature range of ll4-130Kwas covered four times by measuring 
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vnth deci-e.c-ising temperature then increasing, decreasing and increasing c.go'u.. 
The temperature at which measurements v;ere made were interleaved with each 
other to provide a 'closed-loop' of measurements and to check for any change 
in the temperature measuring system; none was found. 

The transition appeared to be complete at 128K when measured with 
increasing temperature. Specific heat measurements, Clusius and Wolf (19^7), 
gave a value of 128. 3K, the difference probably being due to temperature 
gradients and differences of isotopic purity in the sample. In contrast 
to the lower transition, the upper transiton of Dl showed a marked thermal 
hysteresis, the width of the hysteresis loop being 1.75 + 0.5°K. 

The measured variation of ■ ■-p. -j- . -yo r t ho rhomb i c v/ith temperature is 
shown in figure 7-3 and table 7.3- As can be seen, the phase transition 
occurs sharply over 5°K although in reality (because of the temperature 
gradients) it may be still more sharp. This contrasts with the width of 
approximately 30°K shown in the specific heat data of Clusius and Wolf (lS-*7). 

This discrepancy is due to the relative i nsens i t i vl ty of neutron 
diffraction compared to specific heat measurements in revealing the increase 

of molecular disorder occurring in the upper phase transition region. 

, J. 

Table 7.'< lists the variation of the anomalous specific heat C (T) with 

P 

temperature over the upper transition region. To illustrate this difference 
in sensitivity it can be seen in table 7-^ that between the temperatures of 
112K and 12^K the value of C"(T), the anomalous specific heat, doubles but 
the neutron powder peak intensity ratio (which is proportional to l-12p'(T)) 
only changes by 3^- It is only when the specific heat increases by a 
large factor that the intensity ratio decreases by more than ]0%. 



The parameter p'(T) was obtained by graphical integration of the specific 
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T^K 


Di rect ion 


1 [002]^ 


i[ni] 


i[oo2]'^/i [in] 


1 
q'(T) , 


nA.25±0.5 


1 ncreas ing 


39850 


95^00 


0.418±0.02 


1.000±0.05 


n9.25±0.5 


II 


39120 


93^00 


0.419±0.02 


1.000±0.05 


125.25±0.5 


II 


23750 


67500 


0.352±0.02 


0.842±0.05 


126.25±0.5 


ti 


33150 


94310 


0.35I±0.02 


0.839±0.05 


126.?5±0.5 


II 


2^1250 


7V4 50 


0.326±0.02 


0.779+0.05 


127.25±0.5 


II 


19255 


73280 


0.263±0.03 


0.629+0.07 


127.50+0.5 


11 


15^88 


6 80 40 


0,228+0.03 


0.545±0.07 


127.75±0.5 


1 1 


ll^t^S 


70860 


0.162±0.03 


0.387+0.07 


128.00±0.5 


II 


10785 


68980 


0.156+0.03 


0.373+0.07 


128.25+0.5 


II 


14890 


71900 


0.021±0.05 


0.050±0.12 


128.75±0.5 


II 





$k3)0 


0. 00010.0!) 


0.000i0.]2 ' 


126.75+0.5 


Decreasing 


hhB 


70320 


0-006±0.05 




126.50±0.5 


II 


5016 


70484 


0.071±0.05 




126.25±0.5 


M 


IOIA5 


70720 


0.144±0.03 




126. 00 ±0.5 


II 


1^698 


69240 


6.212±0.02 




125.75±0.5 


1 1 


15813 


71080 


0.223±0.02 




125.50+0.5 


II 


I6950 


70900 


0.240±0.02 




125.25±0.5 


II 


20090 


71 800 


0.280±0.02 




123.25±0.5 


II 


s'igos 


92400 


0.378±0.02 





TABLE 7.3. Intensities of neutron powder peaks measured over the upper 

phase transition of DI with increasing and decreasing tempera- 
ture and the ratio and q'(T) values derived from them. 
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T°K 


c;(T) 


P'(T) 


q'{T) 


T°K 


rm 


P ' (T) 


q'(T) 


90.0 


0.00 


.0001 


.9988 


110.0 


1.75 


.0013 


.9a/;4 


90.5 


0.00 


.0001 


.9988 


110.5 


1.85 


.0014 


.9832 


91.0 


0.00 


.0001 


.9988 


111.0 


1.90 


.0015 


.9820 


91.5 


0.10 


.0001 


.9988 


111.5 


1.95 


.0016 


.9808 


92.0 


0.15 


.0001 


.9988 


I112.O 


2.05 


.0017 


.9796 


92.5 


0.20 


.0001 


.9988 


112.5 


2.10 


.0018 


.9784 


93.0 


0.30 


.0001 


.9988 


113.0 


2.20 


.0019 


.9722 


93.5 


0.35 


.0001 


.9988 


113.5 


2.25 


.0020 


.9760 


9^.0 


O.AO 


.0001 


.9988 


114.0 


2,35 


,0021 


.9748 


9'i.5 


0.A5 


.0001 


.9988 


114.5 


2.40 


.0022 


.9736 


95-0 


0.50 


.0002 


.9976 


115.0 


2.50 


.0023 


.9724 


95.5 


0.55 


.0002 


.9976 


115.5 


2.60 


.0024 


.9712 


96.0 


0.60 


.0002 


.9976 


116.0 


2.70 


.0025 


.9700 


96.5 


0.65 


.0002 


.9976 


116.5 


2.80 


.0026 


.9688 


97.0 


0.70 


.0002 


.9976 


117.0 


2.85 


.0027 


.9676 


97.5 


0.80 


.0002 


.9976 


117.5 


2.90 


.0028 


.9644 


98.0 


0.85 


.0002 


.9976 


118.0 


3.00 


.0029 


.9652 


98.5 


0.85 


.0003 


.996^* 


118.5 


3.05 


.0031 


.9628 


99.0 


0.90 


.0003 


.996^ 


119.0 


3.15 


.0032 


.9616 


99.5 


0.90 


.oooA 


.9952 


119.5 


3.20 


.0034 


.9592 


100.0 


1 .00 


.000^ 


.9952 


120.0 


3.25 


.0035 


.9580 


100.5 


1.00 


.OOOA 


.9952 


120.5 


3.35 


.0037 


.9566 


101.0 


1.00 


.0005 


.99^0 


121.0 


3.45 


.0038 


.9544 


101.5 


1.05 


.0005 


.99^0 


121.5 


3.55 


.0039 


.9532 


102.0 


1.05 


.0005 


.99^0 


122.0 


3.60 


.0041 


.9508 


102.5 


1.10 


.0006 


.9928 


122.5 


3.70 


.0042 


.9496 


103.0 


1.10 


.0006 


.9928 


123.0 


3.75 


.0044 


.9472 


103.5 


1.15 


.0006 


.9928 


123.5 


3.80 


.0046 


.9448 


lOi+.O 


1.20 


.0007 


.9916 


124.0 


4.00 


.0048 


.9424 


10'k5 


1.20 


.0007 


.9916 


124.5 


4.20 


.0049 


.9412 


105.0 


1.25 


.0008 


.990^ 


125.0 


' 4.40 


.0051 


.9388 


105.5 


1.30 


.0008 


.9904 


125.5 


4.80 


.0053 


• 9364 


106.0 


1.30 


.0009 


.9892 


126.0 


5.10 


.0056 


.9328 


106.5 


1.35 


.0009 


.9892 


126.5 


48.70 


.0079 


.9052 


107.0 


l.'iO 


.0010 


.9880 


127.0 


92.40 


.0127 


.8476 


107.5 


1.'+5 


.0011 


.9868 


127.5 


136.00 


.0211 


.7468 


108.0 


1.50 


.0011 


.9868 


128.0 


179.70 


.0349 


.5812 


108.5 


1.55 


.0012 


.9856 


128.2 


668.00 


.0484 


.4192 


109-0 


1.65 


.0012 


.9856 


128.3 


51.40 


.0519 


.3772 


109.5 


1.70 


.0013 


.9844 











TABLE 7-4. Te:nperature dependence of the anomalous specific heat over the 
upper phase transition of Dl. The values in this table were 
calculated from the data of Clusius and Wolf (1947). 
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heat data and by nurrierical interpolation of equation (1^(2). Figure 7.^ 
shows two curves; 



a) 



b) 



[1-7.8p'(T)]' 
l-]2p'{T) 



(1^3) 



together with the ratio q'(T) derived from the neutron diffraction data. 
A residual order parameter p' (T ) of value 0.0001 was found to give the best 
overall fit between curve (b) and the order parameters derived from the 
neutron diffraction measurements. Curve (a) is only valid when p'(T) is 
small, i.e. at la>/er temperatures when co-operative effects between neigh- 
bouring molecules are small and short-range order is negligible. 

The poor agreement between these curves and the neutron diffraction 
data can be attributed predominantly to short-range order and partly to 
temperature gradients in the sample. The thermal gradients cause a rounding 
or softening of the order parameter curve as a function of temperature v/hen 
observed by neutron diffraction. 

In the twelve-fold disordered high phase a deuterium atom is assumed 
to lie along one of the directions pointing towards the tv/elve nearest- 
neighbour iodine atoms of its parent iodine atom. On the average each of 
these sites Is occupied with equal probability. This model v^fas used in 
the derivation of the conf I gurat ional entropy change in equation (lA2) and 
it took no account of electrostatic repulsion between deuterium atoms of 
neighbouring molecules. Deuterium atoms in adjacent Dl molecules are 
unlikely to lie along the same bond and this restriction modifies the 
disposition of sites available to each atom. 



The change in conf I gurat ional entropy was calculated by equation (1^2) 
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making the approximation that at a given instant all tv.'elvc sites are 
equally probable. Over a longer period all twelve sites of a given 
molecule are equal ly probable but not as an instantaneous neighbour of 
another deuterium atom (in the next molecule) due to electrostatic repulsion. 

The change in conf igurat ional entropy over the phase transition was 
calculated from the specific heat measurements of Clusius and Wolf (1947) 
and was found to be 2. IAS jmole ^ l<1n 2.9^. The expected change in 
conf igurational entropy from a two-fold to a twelve-fold disordered structure 
is 3.566 jmole" = kin 6. An entropy of 1.419 jmole e kin 2.04l there- 
fore remains due to short range ordering. 

It was shown by Schneider (1955) that, when hydrogen halides form 
chains, the molecular bonds are not col linear with one another but are at 
an angle. This angle is determined by the lone-pair electron orbitals in 
the outer shell of the halogen and for HI/DI is approximately 90°. In the 
low phase of Dl this can be seen as the square rings in the proposed ordered 
model and in the middle phase as the two-fold disordered zig-zag chains. 
Using this information a very simple model can be proposed to describe the 
short range order. 

It is assumed that In the high phase short chains of L molecules joined 
at 90° to each other are formed. Each chain has one molecule at its head 
so that the second molecule can be chosen twelve ways. The next molecule 
must be at 90° and it can be chosen in two ways. Hence the number of ways 
of choosing a chain of L molecules is; 

12 X 2^"2 (145) 

!n the sample there sre H ir-.o'ecules and i!/L chains and the totol nunib-r 
of ways of choosing these chains is; 

[12 X a'--^]""" (iw) 
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The total number of ways of arranging all N molecules in a chain is; 

12 X 2^~^ ^ 2^ (if N is very large) (U7) 

L-2 ^/^ 
. ^ = I12JL1-J (148) 

2' 
The conf igurational entropy per mole due to the chain formation is; 

klnfi = klns'^*- (1^*9) 

and the additional entropy is; 

l.iil9 jmole"' ~ k]n 2.0^1 = kln3'^^ (150) 

so that L - 1.5. 

Therefore the average length of a 90° zig-zag chain is 1-5 molecules. 

7.3 Conclusion 

The lower and upper phase transitions of Dl were studied by neutron 
powder diffraction. Order parameters describing the average configuration 
of disordered molecules were obtained as a function of temperature from 
these neutron powder peak intensities. A similar order parameter was 
derived by integration of the anomalous specific heat over each phase 
transition region. These two parameters agreed well with each other over 
the lower phase transition region, confirming that the lower transition was 
of second order from ordered to 2-fold disordered with an indication of a 
small thermal hysteresis. There was little evidence of short-range order 
in the low temperature region of the middle phase. 

The upper phase transition appeared to be of first-order with a therma 
hysteresis of approximately 1.75°K. The order parameter as a function of 
temperature for this transition obtained by integration of the specific heat 
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did not agree well with that obtained by neutron diffraction. This was 
attributed to some residual short-range order in the cubic high phase. A 
simple model showed that this extended on average over one to two neighbouring 
molecules. 
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CHAPTER EIGHT 

SUMMARY 



8.1 Conclusions of earlier experimental studies conce rning the crystal 
structure of PI 

Experimental investigations carried out prior to this thesis produced 
many interesting results directly relevant to the crystal structure of the 
H!/DI system of which the main results can be summarized as follows. 

!t was the specific heat measurements which first indicated that Dl had 
three solid phases and two well-defined transition temperatures. For 
approximately 20°K below each transition point the specific heats increased 
rapidly then fell back Abruptly thus forming the so-called 'lambda-points'. 
These anomalous changes in the specific heat were Interpreted as due to the 
onset of orientational disorder as either free or hindered rotation of the 
Dl molecules. 

X-ray diffraction experiments indicated that the unit cells were face- 
centred tetragonal (in the low and middle phases) and face-centred cubic 
(in the high phase). In each case the edgelengths of the unit cell were 
approximately 6^ and it contained four Dl molecules. Experiments using 
polarized light indicated that the low and middle phases were anisotropic 
and the high phase was isotropic. These results were in agreem.ent with the 
ideas of two tetragonal phases and one cubic phase from X-ray data. 

Conclusions derived from measurements of the static dielectric constant 
as a function of temperature were similar to those obtained from the calori- 
metric data. An increase in the static dielectric constant on approaching 
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the lower phase transition point was followed by a continuous decrease and 
a sharp drop at the upper transition point. This behaviour suggested an 
increase in the structural disorder at each transition. 

The measurement of the frequency dependence of the dielectric constant 
provided relatively little direct information on the crystal structure. 
Its main contribution was the observation that the complex dielectric con- 
stant of solid Dl was similar to that of several simple polar liquids such 
as glycols. It also showed that the phase transitions were co-operative 
phenomena with re-orientation frequencies higher than those of a typical 
"stiff" molecular vibration. 

Earlyi NMR work indicated that the phase transitions were from ordered 
to disordered structures with several equilibrium orientations rather than 
with freely rotating molecules. Later NMR worl< provided more detailed 
measurements of the different relaxation frequencies which were in general 
agreement with those obtained from the frequency-dependent dielectric 
constant measurements. 

Investigations of the infra-red spectra indicated that the low phase 
of Dl was different from those of HCl and HBr and they also suggested that 
the unit cell was larger at low temperatures. 

Neutron diffraction experiments on solid DBr revealed that the 
deuterium atoms occupied several equilibrium positions In the orientation- 
ally disordered middle and high phases. As the calorimetr ic , dielectric 
constant and NMR data for the compounds DBr and HI/DI were similar this 
provided indirect evidence for the occurrence of orientational ly disordered 
structures in Di. 

By combining all the available experimental evidence one could conclude 
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that the structure of Dl changes from an ordered to successively more dis- 
ordered forms with increasing temperature. The disorder vrauld probably 
be orientational rather than rotational; although the directional hydrogen 
bond would be weal<er in HI/DI than in the other hydrogen halldes. The 
edgelength of the unit cells would be approximately 6^ with its symmetry 
changing from tetragonal to cubic at the upper transition point. 
Co-operative phenomena would give rise to a phase transition extending over 
10-20°K at the lower transition point and a narrower region at the upper 
transition point. 

8 . 2 The cry stal stru c ture of Dl as determined by X-ray and neutron 
diffraction 

The diffraction experiments described in this thesis using X-rays 
and neutrons mostly confirmed the expectations concerning the structures of 
each phase of D!. X-ray experiments on the low phase showed that the 
iodine atoms formed a face-centred orthorhombic unit cell of approximate 
size 6^. However, neutron experiments (vjhich reveal the presence of 
deuterium as well) showed that the true unit cell of the low phase was 
C-centred monoclinic and four times larger in volume than the X-ray unit 
cell. The structure of this phase has not been determined unambiguously. 
A model has been proposed whose calculated diffraction pattern agrees 
reasonably v/el 1 with the observed neutron powder diffraction pattern. The 
iodine atoms lie in face-centred positions of the iodine sublattice while 
the Dl rix)lecules lie in the (OOl) plane. Sets of four molecules exist as 
rectangles which form a continuous stack when viewed along the monoclinic 
C_-axis of unit cell No. 4. 

The middle and high phases were found to be B-centred orthorhombic 
(space group Bbcm) and face-centred cubic (space group Fm3rn) , both structures 
being isomorphic with the corresponding phases of DBr. The edgelengths 
of the unit cells were approximately SK with the middle phase only slightly 



197. 



distorted from a tetragonal cell. The middle phase was two-fold disordered 
with the deuterium atoms having two equilibrium sites at 180° to each 
other. The high phase was twelve-fold disordered with the equilibrium 
sites lying along the twelve nearest-neighbour iodine-iodine bonds. 

The three phases were studied by X-ray and neutron diffraction at 
ii.ZK, SkK and 137K and the structural parameters refined by least-squares 
methods. Anisotropic temperature factors were obtained wherever possible, 
although in the high phase correlation between parameters did not completely 
allow this. 

The lattice parameters were measured at several temperatures by X-ray 
methods and are shown in figure 3.5- The mean linear thermal expansion 
coefficient of Dl was found to be 3-5 + -5 x lo"^ K^~^ As the temperature 
increases the unit cell expands uniformly in the low phase. At the lower 
phase transition point there is no noticeable change in the lattice 
parameters as the Di molecules begin to 'flip' in the (OOl)^ plane. As 
the temperature rises further in the middle phase the £-axis increases 
while the a_ and b_ axes decrease. Near the upper phase transition point 
the lattice parameters converge rapidly to the cubic high phase while the 
volume of the unit cell increases continuously with temperature. In 
general the lattice parameters of Dl are shorter than those of HI, this is 
true also for HBr/DBr and HCl/DCl, because the amplitudes of vibration of 
the deuterium atoms are smaller than those of the hydrogen atoms and con- 
sequently the lattice is packed more tightly. 

8.3 The natur e of the two phase transitions in solid Dl 

At each phase transition neutron powder peaks which were disappearing 

ever the transition region v/ere studied as a function of temperature in 
relation to a stable reference peak. In this way the lower phase transition 
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tji 

was found to extend from 45-80K with a thermal hysteresis whose width was | 

approximately 0.5°K. The upper transition (as observed by neutron powder | 

diffraction) appeared to extend from 122-128Kwith a well-defined thermal |lj 

hysteresis of width 1,5+0.5°K, although the specific heat data indicate I 

that the upper phase transition may be a great deal sharper than observed jj 

1! 
in this neutron study. li 

■ - i| 

Order parameters were obtained from these data as a function of 
temperature. The anomalous specific heat obtained from the calorimetric 
data was integrated to obtain an equivalent order parameter. The order 
parameters derived by these two methods were in good agreement over the 
lower transition. However, the upper transition neutron order parameter : 
showed a more gradual change at the transition point than that calculated 
from the calorimetric data. The difference was possibly due to short-range 
order between the deuterium atoms which are able to form short chains con- 
sisting of a few molecules in the disordered high phase. The overall 
changes in conf igurational entropy between the three phases were evaluated 
as 1.01 H kin 1.74 (lower) and 2.1^(6 = l<ln 2.3^ (upper). The conf i gurat ional 
entropy change associated with the lower phase transition is consistent 
with the proposed structural change from ordered to two-fold disordered. 
Moreover, it suggests that the two equilibrium orientations in the middle 
phase are occupied alnxjst at random, without appreciable short range order. 
On the other hand, the conf i gurat ional entropy change associated with the 
upper transition is considerably smaller than would be expected if the 
transition from two- to twelve-fold disordered were entirely random, with 
no short-range order. Hence the calorimetric data combined with the 
neutron diffraction results suggest the presence of appreciable short range 
order in the high phase of Dl. 
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8.^ Suggested further experiment s on HI/DJ_ 

a) Neutro n pov/der diffraction 

Instruments and facilities in experimental sciences often improve in 
large finite steps. This Is particularly true of "big" science where 
facilities such as reactors and neutron dif f ractometers cost respectively 
tens of millions and tens of thousands of pounds. Over the period of 
three years in which the experimental work for this thesis was carried out, 
the most noticeable change was the considerable improvement in the angular 
resolution of the powder dl ff ractometers at A.E.R.E., Harwell. In 
addition, this Improvement was enhanced by the ability to vary the mono- 

chromator take-off angle on instruments such as Pluto 'Panda' and so to 

. • r- sin6 , 
optimize the angular resolution in the region of — r— where it was most 

useful . 

Increases In angular resolution are always paid for by a reduction in 
the scattered Intensity. High flux beam reactors such as the H.F.R. at 
the Institute Laue-Langevin , Grenoble, France, where the incident neutron 
flux Is at least one order of magnitude higher than at Harwell, will 
with the development of suitable Instruments provide increased intensities 
and higher angular resolution. In addition pulsed reactor or accelerator 
facilities provide opportunities for improved resolution and 'instantaneous' 
data recording. 

Such 'instantaneous' recording can be achieved by the use of position 
sensitive detectors and electronic t ime-of-f 1 Ight instruments. In such 
experiments the Instrument behaves as an electronic powder camera and 
records the whole accessible range of — r— simultaneously. This type of 
Instrument, in addition to reducing the tim.es of experiments, greatly 
simplifies the study of phase transitions as a function of temperature and 
reduces many uncertainties concerning counter efficiencies, etc. 
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The phase transition experiments described in this thesis each took 
approximately three weeks of reactor time and the temperature dependence of 
only a few peaks was studied. Mul t idetector instruments as available at 
I.L.L. and under construction at A.E.R.E. would provide complete diffraction 
scans at improved statistical accuracy, angular resolution or narrov'.'er 
temperature intervals. These data could be analysed using computer pro- 
grams of the type developed by l^ietveld (1967). With an order of magnitude 
increase in reactor flux an experiment previously taking over 2 months 
(which would rarely have been considered in the U.K.) would take one week. 

b) Single-crystal neutron diffraction 

Single crystal diffraction measurements on HI/DI would provide more 
detailed information on the crystal structure and particularly temperature 
factors and the related thermal motion parameters. it appears difficult to 
grow such crystals although some success has been reported for DCl, 
Hoshino (1970), Shimaoka (1970 and Maruyama (1970). High phases of 
molecular crystals can be grown from the liquid by the Bridgmann (1925) 
method. However, on cooling to lower tempeatures through the phase transit- 
Ions they invariably fragment into twins and multiple crystallites- The 
low phase of Dl would most benefit from a single-crystal study out of all 
the hydrogen halides. This is due to its monoclinic symmetry on a pseudo- 
orthorhombic lattice as described in Chapter four. A powder diffraction 
experiment has a built-in degeneracy which combines several different 
peaks with the same d-spaclng. Only a single-crystal experiment would be 
able to separate these reflections. 

c) Inel astic neutron scattering ^ 

When suitable crystal -growing techniques have been developed the 
hydrogen halides will provide a very interesting field of study for 
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inelastic neutron scattering. Improved neutron sources as discussed 
earlier are even more welcome in inelastic scattering where many experiments 
involving energy analysis of the scattered neutrons become possible only 
with these high fluxes. 

In particular, inelastic neutron scattering can provide information 
on the lattice modes of the crystals and enable the study of the 
re-orientation processes occurring in the middle and high phases. Such 
comprehensive data on the frequency of these re-orientations and the 
associated mechanisms and forces will, it is hoped, help the development of 
a satisfactory theory of the phase transitions. 
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APPENDI X ONE 

Computer Progra ms 

Four different computers were used during the course of this work: 

(a) The Science Research Council ICT, Atlas Computer at the Atlas 
Computer Laboratory, Chilton, Oxon. 

(b) The London University ICT, Atlas Computer at the Institue o.f Computer 
Science, Gordon Square, London. 

(c) The London University CDC 66OO computer at the University of London 
Computer Centre, Guilford Street, London 

(d) The IBM 370 computer at A.E.R.E. , Harwell. 

All programs used were written in subsets of FORTRAN, at one stage 
sections of ORFLS were rewritten from FORTRAN V to FORTRAN iV for 
compatbility v^ith all computers. 

Al . 1 Least-s qu ares refinements of cryst al structures 

The least-squares refinements of crystal structures were performed 
using a modified version of ORFLS, Busing (I962). The modifications are 
described elsewhere by Ogunade (I969) and Johnson (1971), and a brief out- 
line only will be given here. 

Because of the relatively poor angular resolution of the neutron 

powder di f f ractometers and the near-equality of the lattice parameters of 

Dl, this ratio being typically 1:1.01:1.09, it was not possible in the low 

and middle phases to obtain a set of intensities for single powder peaks. 

The modified version of ORFLS performs least-squares refinements with 
sets of intensities grouped as multiplets. ORFLS minimizes: 
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NO 2 
E = y W A where NO = Number of observations 



j=i '■'■ 



t h 

W = V/eiqht of £ observation 



and A = F', -sF' , where F' = Observed structure factor 
I obs calc obs 

F' , = Calculated structure factor 
calc 

or A„ :=1f', I^-IsF' , 1^ s = Scale factor 
% * obs' ' calc' 



in the modified program, MP-ORFLS minimizes: 



NO 2 NPK . / p, si 

I W A^ where A, = l^bs " ^ ^Zf.^!^ ^^^D 

ji=1 ^ ^ ^ °^^ ^\ sine sin20 



where NPK = Number of peaks in multiplet 

t h 

j = Multiplicity of m peak 

9 ~ Bragg angle of m component peak 

NPK . C.2 
and I bs == ^ •'"^ ^bs (152) 

° m~l sine sin2e 

mm 



Further modifications were made to this program at a later date and 
included; extended error messages, optimization of the available core 
store and a printing of the calculated Bragg angle beside each peak to 
enable easy checking of the Miller indices. This was useful'when dealing 
with a monoclinic unit cell. 

A 1.2 L east-squ ar es refine ments of cell dimensions 

Least-squares refinements of lattice parameters v;ere performed using 
the program PARAM. This is a part of the overlaid suite of programs 

called X-PiAY 70 (Stewart 1970) which are stored on magnetic disks at all 
three computer establishments. 
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PAR-AM minimizes E where: 



sin L sine ,J-e, 
caic obs 



The program performs 5 cycles of least-squares refinement and, 
depending on the type of lattice specified initially, varies only the 
necessary parameters. 

Al .3 Calculation of powder peak positions 

The program THETA was used to calculate the full set of allowed Bragg 
reflections for a particular lattice. Symmetry conditions of two types 
could be imposed to provide constraints on the allowed reflections, 
e.g. Type 1: A single letter, F,I,A,B,C ... 

Type 2: More general conditions 00£ (£=2n) etc. 

A 1 . 4 Calculation of random orientation structure factors 

A Simpson's rule method of integration was used to evaluate the 
integrals derived in Chapter 3- 

^^•5 LPSG, Low phase symmetry gen erator 

This program used an analytical expression for the structure factor 
of spare group Cc in which four independent molecules were chosen as the 
basis. An inner loop produced the first 13 Bragg reflections and four 
outer loops varied the molecular orientations to generate the separate 
128 structures. 

Each structure was numbered and contained a listing of the deuterium 

atomic co-ordinates. 
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^^•^ Xhe R i et ve 1 d p rog ram 

This program was written for the refinement of neutron diffraction 
powder patterns by Rietveld (1967, 1969). The implementation used here 
is that of Hewat (1973)- The program computes the neutron diffraction 
powder pattern to be expected from a model crystal structure and a 
particular diffractometer. The following parameters are supplied by the 
user; 

(a) Describing "the crystal structure 
(i) Symmetry operations. 
(ii) Atomic co-ordinates, scattering lengths, atom mul tipl iers 
and temperature factors. 
(ill) Scale factor. 
(iv) Unit cell dimensions. 

^^^ Describing the diffractometer 
(i) Wavelength, 
(il) Zero-angle of the 26 scan. 
(iii) The angular resolution as a function of angle. 

The width of the powder peak (FWHM) is described by a function; 



H^(e) = U tan^e + V tane + W (153) 



The program minimizes by successive Iterations the following quantity; 

y. 2 
X^ = I w.[y.(obs)-":!-(calc)] (154) 

.11 c 

The summation Is taken over all the points of the diffraction pattern 
where y.(calc) and y. (obs) are the calculated and observed intensities of 
weight w. and c is the scale factor. 
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Any or all of the parameters listed above may be varied during the [ 
structure refinement. This program has the advantage that all the | 

information contained in the unresolved multiple peal<s of a complicated \ 
powder diffraction pattern can be utilised in refining the structure. The 
complexity of these peaks is due to the poor angular resolution of neutron 
powder dl f f ractometers necessitated by their low Incident flux compared 
with X-ray dl ff ractometers. 

« 

The disadvantage of the program is that convergent refinements only 
occur when all the parameters are close to their correct values, I.e. the 
user must l<now all the instrumental parameters accurately and the structural 
parameters tolerably v^^ell before refinement. 
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APPEtsiPIX TWO 

Temperflture measurement and co ntrol circ uit of the variable temperature 
cryostat 

The control circuit is shown in figure A2.1 and the cryostat system 
in figure 2.10. The sample is held in a vacuum insulated exchange space 
at S. The sample temperature measuring thermocouples, of O.OOSin. diameter 
copper and constantan wire, are fixed at T^ and T2 and pass through neoprene 
vacuum seals, Batchelder (I969) at Si and S2 to a liquid helium reference 
junction, H. The thermocuple leads then pass to 

(a) A chart recorder C2. 

(b) A temperature controller, V/, manufactured by the V/est Division 
of Gulton Industries Ltd., Brighton. This was a West Viscount, 
three term controller, consisting of a preset potentiometer and 
a sensitive balance circuit. The heater output current was 
automatically varied to maintain a constant thermocouple emf. 
The manufacturers specified a control stability of lOuV under 
optimized conditions corresponding approximately to 0.5°K at 8OK. 
However, the recorded stability on the other measuring devices, 
including the centre reading galvanometer in the control 1 er, 

was nearer 2yV. 

Three additional parameters could be varied. 

(1) the percentage of feedback control in the circuit; 

(2) the integration time of the circuit, i.e. the response time of the 
measuring system; 

(3) the differentiation time of the measuring circuit. 
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The feedback control v;as set at the recommended value of 10^. The 
thermal response time of the sample system v;as recorded and found to be 
approximately 5 minutes. The integration time was set to this value and 
the differentiation time to 1/3 of this value as recommended by the 
manufacturers. 

Cent inuous thermocouple leads were used throughout, with no soldered 
joints other than at the measuring junctions. This avoided any errors due 
to thermal gradients at vacuum seals, etc. The leads to the controller 
were screened by an earthed cable, and the junctions to the controller were 
insulated with cotton wool. 

The level of liquid helium in H was monitored by another copper- 
constantan thermocoup^le T3 with a liquid nitrogen reference junction N. 
The emf was recorded on another chart recorder CI which was connected via 
relays to an alarm bell and to a micro-switch in series with the heating 
circuit. The thermocouple T3 was positioned approximately 10cm above the 
Tp T2 reference junctions. The cryostat was filled v/ith liquid helium 
and nitrogen every 8 hours for the duration of the experiment. 

The liquid helium transfer lines were always precooled before filling 
with liquid helium. Due to the configuration of apparatus in the Dido 
shell it was necessary to operate from a gantry some 10ft above the cryostat, 
using multi-sectioned transfer lines. At no time during the experiment 
was the liquid helium level allowed to fall below the level of the 
monitoring thermocouple T3. 

The temperature controller W provided a controlled heater current at 

25V D.C. This v/as passed through the heaters in series and a neoprene 
vacuum seal S3. Reistances Rl and R2 (both I6OQ) were wound on fine 
Eureka wire above and below the sample space S. Resistance R3 (500Q) acted 
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as an antisubl Imation heater. This prevented any sample sublimation to 
cooler parts of the tube while performing experiments at 78K. Because the 
sample space was supported only from above, a thermal gradient existed down 
the tube between T, and Tp to which the heat from R3 also contributed. 
For this reason a trimming potentiometer Rk (50Kfi) was used to balance the 
currents in R1 and R2. Unfortunately the system was only sensitive to 
this adjustment over a period of 2 hours. In order to complete the 
experiment within the allocated period a compromise technique of maintaining 
a minimum gradient at 78K was adopted. 

A graph showing the variation of this gradient with the mean tempera- 
ture is shown in figure A2.2. The gradient measured from the lower 
thermocouple T» increases with lower mean temperature due to the decreased 
effect of Rl and R2. The lower thermocouple T^ was positioned ISmm above 
the centre of R2 and T| was positioned 9mm below the centre of Rl . This 
accounts for the weighted value of AT in fgure A2. 2. 

The thermocouple emf was also recorded on a digital voltmeter D and 
printed on the teleprinter output of the di ff ractometer. 

The radial temperature gradient in the DI sample system was measured in 
the following way. A small diameter glass probe was inserted from the 
top of the cryostat system into the sample space v/ith a thermocouple inside. 
DI was then sublimed into the sample space and a gradient of approximately 
0.5°K was measured between the probe and T, . However, the true gradient may 
be greater than this value because the probe and thermocouple provided a 
substantial input of heat from outside the system which would not be 
present in an actual experiment. 
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Figure A2.2 Temperature gradient across sample, as a function of 
temperature, for the sample tube used for the lower 
phase transition neutron diffraction experiment on DI , 
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